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Abstract 

Chaplygin [9^ proved the integrability by quadratures of a round sphere, rolhng without 
shpping on a horizontal plane, with center of mass at the center of the sphere, but with 
arbitrary moments of inertia. Although the system is integrable in every sense of the word, 
it neither arises as a Hamiltonian system, nor is the integrability an immediate consequence 
of the symmetries. On the other hand, the constants of motion are obtained as a consequence 
of Noether's principle, cf. Section [T] and El The system also turns out to be related to a 
Hamiltonian system, the geodesic flow on the Euclidean motion group for a left invariant 
metric, cf. Subsection 19.31 

In this paper we analyse the global dynamics of Chaplygin's sphere. In the process we 
wiU explain almost all of Chaplygin's results. Readers who are mainly interested in these 
may consult Sections [2 El El El Subsections 19.21 111.21 111.31 111.51 and Section |T21 These can 
be read independently from the rest of the paper with some exceptions, such as Subsection 
19.21 in which Subsection 19.11 has been used. At the end of each section we describe in a 
subsection "Chaplygin" the relation between our text and Chaplygin's. 

We also obtain some new results, such as the proof in Section jH that the level sets of the 
constants of motion in the reduced phase space are two-dimensional tori. In Section |H1 we 
prove that, after a suitable time reparametrization, the rotational motion is quasi-periodic on 
two-dimensional tori. After suitable completion of the level surfaces, this is also true for the 
complexified system, cf. Sectional This shows that the rotational motion is algebraically 
integrable according to the definition of Adler and van Moerbeke. In Subsection 111.41 it is 
explained how this also follows, in a quite different way, from Chaplygin's integration in 
terms of hyperelliptic integrals. 

1 Noether's Principle for Nonholomic Systems 

We use the equations of motion for a system with nonholonomic constraints as given by 
d'Alembert's principle 

{[Lnt), v) = for every v G C^(t), (1.1) 

*A large part of this work has been done during a sabbatical leave in Berkeley, in the fall of 1994, partially 
supported by AFOSR Contract AFO F 49629-92. A more recent version was prepared in July 2000, as a 
chapter in the planned book "The geometry of Nonholonomically Constrained Systems" , together with R.H. 
Cushman and J. Sniaticky. 
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as described in jU Ch. 1, Sec. 2.5]. We may assume that Qq = Q. Let s ^— 7,, be a smooth 
family of smooth curves in Q, with 7 = 70, 6{t) := d'ys{t)/ds\s=o, for which we this time do 
not assume that 6{a) = and 6{b) = 0. Write 



for the 5-component of the momentum, which is a coordinate-invariant quantity. The classical 
variational equation, which is obtained by a partial integration, reads 

^ L(7.(t), m) dt L.0 = -j\[LV{t), S{t)) dt + /(6) - /(a). (1.2) 
If one differentiates ()1.2|1 with respect to b, then one obtains the equivalent form 

y/it) = yL (7.(t), I'M Lo + (W^W, 5(t)). (1.3) 

dt ds 

The second term in the right hand side of (jl.3|) can be viewed as the 5-component of the 
reaction force acting on the system. It is equal to zero if 7 satisfies d'Alembert's principle 
p.l|) and 6{t) G C^(t), which means in words that 6{t) is a virtual displacement. 

The form ()1.3|) of the variational equations is due to Emmy Noether in a version 
where the independent variable t is replaced by a finite number of real variables and L is a 
smooth function on a jet bundle of arbitrary order. In the case that the first term in the 
right hand side of (jl.3p is equal to zero, the formula (|1.3p . is the "momentum equation" of 
Bloch e.a. Thm. 4.5]. 

Suppose that w is a smooth vector field on Q such that ■;^7s(t)|^^o = w{'y{t)). Let w 
denote the lift of w to T Q, the vector field w onT Q such that 

e^'^ = T(e'"'), (1.4) 

if e*"* denotes the flow after time s of the vector field w. In local coordinates w is given by 

w{x, v) = {w{x), Dw{x) ■ v) , (x, v) eTQ, (1.5) 

where Dw{x) denotes the matrix dw^{x) / dxK With this notation, the first term in the 
right hand side of (jl.3|) is equal to the derivative of L at (7(t),7'(t)) in the direction of 
w. This leads to the following version of Noether's principle for variational systems with 
nonholonomic constraints. 

Lemma 1.1 Let L he a smooth function on TQ, of which C is a smooth vector subbundle. 
Let w be a smooth vector field on Q with the following properties 

i) w is a virtual diplacement, which means that w is a section of C . 

ii) At each point of C , the derivative of L in the direction of w is equal to zero, where w 
is the lift of w to T Q as defined by fll.4|) . 
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Then the w-component of the momentum is constant along every solution of 

If there are no constraints, when C = TQ, then condition i) is void and ii) is equivalent to 
the condition that L is invariant under the flow of the vector field w inTQ, which is equal 
to the tangent lift of the flow of w in Q. In this case Lemma ll.ll is due to Emmy Noether 

Question 1.2 Can all the constants of motion in Chaplygin jH] be obtained as applications 
of Lemma 11.11 .'' 

Remark 1.3 Lemma 11.41 below leads to the warning that in the non-integrable case the 
condition ii), under the assumption that i) holds, is not a property of only the restriction of 
L to C, because at the points of C the vector field w need not be tangent to C. 

Lemma 1.4 For each section w : Q ^ C of C the vector field w is tangent to C if and only 
if the subbundle C ofTQ is integrable. 

Proof Let be the flow of w. Then the condition that the lift of w is tangent to C is 
equivalent to the condition that the mappings Tip^ leave C invariant, or that these mappings 
send sections of C to sections of C. This in turn is equivalent to the condition that [w, u] is 
a section of C for every section u of C. That this holds for every section w of C is one of 
the equivalent Frobenius conditions for the integrability of C. □ 

1.1 Chaplygin 

The version of Lemma Fl. II with nonholonomic constraints can be found in Arnol'd ^ p. 82], 
with condition ii) replaced by the somewhat stronger condition that L is w-invariant. Two 
applications have been given in Arnol'd jH p. 83, 84], the first with a reference to Chaplygin 
[Hj and the second with a reference to Chaplygin [H]. In Chaplygin [HI the constants of motion 
have been described as an application of jS]. 

2 Noether's Principle for Chaplygin's Sphere 

The position of a rigid body is given by a pair (A, a), with A E SO (3) and a E R'^. Thus, 
if X G R^ is the position of a material point of the body in its reference position, then 
y = Ax + a is the position of the corresponding point in the moving body. If s is the 
position on the surface 5* of the body in the reference position, such that p = As + a is the 
point of contact of the moving body with the surface P on which the body is rolling, then 
the condition of rolling without slipping means that 

As + a = 0, (2.1) 
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meaning that the at the point of contact the corresponding material point of the body is at 
rest. Correspondingly, {A, a) is a virtual displacement if and only if 

As + a = 0. (2.2) 

If fi denotes the mass distribution of the body in the reference position, which is a finite 
Borel measure on R^, then the kinetic energy of the moving body is given by 

T = I hAx + d, Ax + d) n{dx). (2.3) 

It follows that the {A, a)-component of the momentum is equal to 

= /" {Ax + d, Ax + a) fi{dx). (2.4) 

Let z/ G R'^ be the unique vector such that 

Az = ux (Az), zeR^ (2.5) 

— note that this corresponds to the right trivialization of the tangent bundle. If the condition 
()2.2|) holds, meaning that {A, a) is a virtual displacement, then 

/-^'^)= [ (y, I(x-s))/i(dx) = (j, z/), 

"^r3 

in which 

j := / fi{dx) (y-p) xy (2.6) 

is the moment of momentum about the point of contact p. Here we have used that Ax+d = y, 
A{x - s) = {y - a) - {p- a) = y - p, and {y, u x {y - p)) = {{y -p) xy, u). 

We now turn to the case of Chaplygin's sphere [0] , where the surface S of the body in the 
reference position is a sphere, the center of mass is at the center of S, and the body is rolling 
without slipping on a horizontal plane P. We will take the origin of the reference frame at 
the center of mass = the center of S. If r denotes the radius of 5*, and we take the plane 
P at height — r, then the condition that ^4(5*) + a is lying on top of P corresponds to the 
condition that the third (vertical) component of a is equal to zero. The point of contact then 
is equal to p = a — r if denotes the third standard basis vector, and the corresponding 
point on S, in body coordinates, is equal to 

s = -rA-^e3. (2.7) 

The condition (j2.2|l therefore is equivalent to 

a = ri' X 63, (2.8) 
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where we have also used ()2.5|) . 

From this moment on, we keep the infinitesimal rotation vector v G constant. Then 
(12. 8p implies that a is a constant horizontal vector. If we use the lejt trivialization of the 
tangent bundle of S0(3), corresponding to assigning to A G TaS0(3) the infinitesimal 
rotation given by 

Az = A{u)X z), 2 G (2.9) 

for some vector uj G R^, then the tangent lift of the vector field A such that ()2.5|) with a 
constant v does not effect uj. Because ()2.8p implies that the tangent lift of a does not effect 
a either, the conclusion is that condition ii) of Lemma 11.11 holds if we take L = T. Note 
that for Chaplygin's sphere the center of mass remains at the same height, which means that 
the gravitational potential energy is constant, and therefore can be disregarded. We have 
arrived at the conclusion that 

Proposition 2.1 For Chaplygin's sphere, the moment of momentum about the point of con- 
tact is a constant of motion. 

The kinetic energy of the rigid body is given by 

T=^{lLU,Lu) + ^m{d,d), (2.10) 

where we have used the left trivialization (j2.9|) of the tangent bundle of SO (3). Here / denotes 
the moment of inertia tensor, which is given by a positive definite symmetric matrix, and m 
denotes the total mass of the body. 

In this notation A, given by (j2.5j) . corresponds to 

^5 = ^-^, (2.11) 

because ^ 

Az = u X Az = A [A~^u X z) . 

It follows that 

= {AIuj, u) + mr'^ {Au x e^, u x e^) 
= {AI UJ + mr^es x {Au x e^) , u), 

where we have used that ()2.1|1 . ()2.7|1 and ()2.9j) imply that 

d = r A(^ij X A'^ 63) = r Auj X 63 

and 6 is given by ()2.8j) . With the notation 

u:=A-^e3 (2.12) 

this leads to the formula 

j = A [I UJ + mr'^u X {uj X u)) (2.13) 
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for the moment of momentum about the point of contact. Note that (m, u) = 1 and therefore 

u X {uj X u) = uj — {u, uj) u, (2-14) 

which is the orthogonal projection of u onto the plane which is orthogonal to u. Also note 
that u has the concrete interpretation that — r u is equal to the point of contact on the surface 
of the sphere, in body coordinates, cf. fl2.7|) . 

In order to simplify the notation somewhat, we write 

p := mr^, (2-15) 

and define the symmetric linear mapping : — > by 

Ip,u{i^) ■= {I + p)^^ - p{u, uj)u, (2.16) 

which is equal to / plus p times the orthogonal projection to the plane orthogonal to u. With 
these notations, we have that 

j = AIp,uUJ. (2.17) 

Write 

J:=iI + p)-\ (2.18) 
If Jp^ cij = then uj = J (u + Ou) for some G R, which moreover has to satisfy the equation 

u = i' + 6u — p{u, J {u + 9 u)) u, 

which holds if and only if 

e {l- p{u, Ju)) - p{u, Jv) = 0. 

This leads to the conclusion that the symmetric linear mapping Ip^^ is invertible, with inverse 
given by 

Ip,^-\u) = Ju+ ^ ^^^Z';^ ^ Ju. (2.19) 
1 — p [u, J u) 

If we fix the constant of motion j, then we can use the equation 

cu = J,,„-M-ij (2.20) 

in order to express uj in terms of A, where we note that u is given in terms of A by means 
of ()2.12|) . In this way uj can be eliminated from the equations of motion. 

2.1 Chaplygin 

In §1], Proposition 12. H is stated "■ ■ ■ as a consequence of a generalized theorem of areas", 
with a reference to [S]. 

In jSl §6] the case of [H] appears as the limit when the radius of the big sphere in which 
the small sphere rolls tends to infinity. 
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3 The Equations and Constants of Motion 



3.1 The Equations of Motion 

The equations of motion are 

dp 



r A{iu X u) = r (Alu) x 63, (3.1) 
dt 

dA , , , 

— = Aou;,^, (3.2) 
dt 

— = uxuj, u{t) := A{t)-^ 63, (3.3) 
dt 

d r T 2 ( I ^'^ \ duj\ 



-luj — I uj X iu = mr [ {u, — ) u . (3.4) 

dt \ dt dt / 

Here, as in Section ^ x ^ Ax + a, x, a & R^, A G S0(3) is the rigid motion which is 
apphed to the body in the reference position, with the center of mass at the origin, and 
p := A ?(«) + a = a — r 63 denotes the point of contact between the sphere and the horizontal 
plane. Furthermore Uop denotes the antisymmetric linear mapping u \—>- uj x u : — > R^. 
The equation ()3.2p expresses that G R^ can be viewed as a rotational velocity vector. 

The full system p.ip . ()3.2|) . ()3.4|) is defined in the eight dimensional (p. A, a;)-space 
R2 X S0(3) X R3. The equations (jSH, (El for 

{A, u) e S0(3) X R3 ~ T (S0(3)) 

are the equations for the rotational motion, the system obtained by reduction by the hori- 
zontal translation group. The equations (j3.3|) . (j3.4|) for 

(u, 00) G xR3 

represent the reduction of the system by the left action of the horizontal motion group E(2). 
With the notations (j2.15j) and (j2.16j) . the equation (j3.4|) takes the form 

Ip,u'^ = (luj) X UJ. (3.5) 
dt 

Combining ()2.16|) with ()3.3p . we obtain that 

dlp^u { I dn^ , \ 



and therefore 



hence 



= —p {{u, U X u) U+ {u, u) U X Oj) = —p {u, u) U X UJ, 



— [Ip,uUj) = ——UJ + Ip^u — = -p{uJ, U}UXUJ + [luj) X UJ, 
dt dt dt 

4 (/p,«tu) = {Ip,u^) X UJ, (3.6) 
dt 



because puj x uj = 
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3.2 The Constants of Motion 

In general the total energy is a constant of motion, when the the equations of motions have 
been obtained as a consequence of the principle of d'Alembert, cf. fll.l|) . Because in our case 
the potential energy mgr is a constant, it follows that the total kinetic energy 

T = i {luj, uj) + ^mr^{u x u, u x u) = ^ {Ip,u^, uj), (3.7) 

cf. ()2.1(ij) . is a constant of motion. This can also be verified directly from ()3.5j) and ()3.fi|l . 

because ^ ^ 

{Ip,u^, — ) = (t^, Ip,u — ) = {uJ, (luj) XUJ) =0 

and therefore also 

at dt 

On the other hand, combination of ()3.2|) with ()3.6|1 yields that 
d 

-(A/p,„oj) = —Ip^^uj + A-{Ip^^uj) 

= A {u X (Ip^uuj)) + A {{Ip^^u) X u) = 0. 

In this way we have verified again that the vector AIp uUJ, which according to ()2.17p is equal 
to the moment j of the momentum around the point of contact, is a constant of motion. 

3.3 A Pair of Vectors 

If j is not vertical, then the rotation A is determined by the pair of vectors 

u:=A'^e3 and v:= A-^ j. (3.8) 

More precisely, in this case the mapping A \—>- {u, v) is a diffeomorphism from SO (3) onto 
the smooth algebraic submanifold of R^, which consists of the [u, v) E R'^ x R'^ such that 

{u,u) = l, {u,v)=js, {v,v) = \\jf. (3.9) 

The equations of motion for the rotational motion are given in these coordinates by 

du dv 

— = u X u and — = v x u, (3.10) 
dt dt 

in which u is determined in terms of u and v by the equation 

iO = uj{u, v) = Ip^u~^ V = J V + ^ "^^^ J u, (3.11) 

I — p {u, J u) 

which in view of fl3.8|) is equivalent to ()2.17|) . Here we have used the equation ()2.19|) in order 
to write uo even more explicitly as a function of u and v. 
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In view of fj3.7|) . the kinetic energy can be expressed in terms of uo and v as 



T = ]^{v,u), (3.12) 
which in view of ()3.1H) and ()2.19p can be written in the form 

T = -{v, Jv) + — . 3.13 

2 2 1 — p {u, J u) 

Later it will turn out to be convenient to write the kinetic energy equation in the form 

f{u, v) := Y{u, - X{u) Z{v) = 0, (3.14) 

in which 

X{u) := p-^ - {u, Ju), (3.15) 
Y{u, v) := {u, Jv), and (3.16) 
Z{v) := 2T-{v, Jv). (3.17) 

Note that f{u, v) is a polynomial of degree four, but of degree two in each of the variables 
u and V separately. 

3.4 The Left SO (2) Action 

If i? is a rotation about the vertical axis, then its action from the left sends A and A to RA 
and RA, respectively. It therefore leaves u and u invariant and sends j to Rj. Note that 
the action of the group SO (2) of the rotations in R'^ about the vertical axis is free on the 
set J^' of j G which are not equal to a multiple of 63. In J'', the SO(2)-orbits are equal 
to the level curves of the functions F{j) = {j, 63) = js and G{j) = {j, j) = where 
j G Jj' corresponds to the condition that F{j)^ < G{j). Substituting (|2.17p we obtain the 
constants of motion 

J3 = (j, 63) = {Ip,uUJ, u) = {luj, u) (3.18) 

and 

\\3\? = {3,3) = {Ip,uu:,Ip,uu:) (3.19) 

for the left E(2)-reduced system for (m, oj) G S^ xR^. 
Let 

TT : {A, uj) y-^ {u, u) = (^A~^ e^, 

denote the projection from the phase space T (SO (3)) ~ SO (3) x R^ of the rotational motion 
onto the phase space S^ xR^ of the left E(2)-reduced system, which maps each left S0(2)- 
-orbit to a point. The fact that the action of S0(2) on J'' is free implies that vr is a 
diffeomorphism from the submanifold of SO (3) x R^ determined by the equation A Ip^uJ = j 
onto the the submanifold of S"^ x R^ determined by the equations (j3.18p and (j3.19p . where 
each of these submanifolds is invariant under motion of the system. The left S0(2)-invariance 
of the system means that vr intertwines the rotational motion with the flow of the E(2)- 
reduced system. 
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3.5 Chaplygin 

In the left column of the the following table we list the variables and some formulas which 
appear in Chaplygin [HJ §2], with our corresponding notations in the right column. It is 
assumed that the moment of inertia tensor / is in diagonal form, in accordance to the 
"principal axes of inertia attached to the sphere" of Chaplygin [HI §2]. 

Chaplygin |21 §2] our notation 

(p, g, r) uj 
{u, V, w) 

iP,Q,R) A-'j=v 

(7, y, 7") es = u 

m m 

p r 

D = m p^ p = mr"^ 

{L, M, N) diagonal of / 
{A, B, C) diagonal of / + p = J"^ 

(1) (EH) 

(2) dSHD 

u in (3) —{^, u) 

X X{u) in (ITTHI) 

Y Y{u, v) in (pnH|) 

(6) {I + p)uj = v + ^u, p{u,u) = ^ 

(7) (Umi) and (USD 

n \\jf = {j^ j) = {v,v) 

h J3 = (j, 63) = { v, u) 

I 2T = {v, u), cf. (I3"T2|) 
Z Z{v) in dSHI) 

(10) dsn 

The only comment of Chaplygin to his formulas (1) and (2) consists of the preceding sentence 
"We easily find 

No explicit notation has been introduced in Chaplygin \^ §2] for the rotation A. However, 
when Chaplygin said "(7) are the equations of motions of the sphere", it is clear that he 
meant our Subsection 13 .31 

The constants of motion of the E(2)-reduced system are Chaplygin's n, h and /, which 
correspond to our Hjp, js and 2T, respectively. 

4 The Level Surfaces of the Constants of Motion 
4.1 Fixing the Moment 

The system of equations ()3.2|) . ()3.4|) in the (A, ci;)-space SO (3) x describes the rotational 
motion of Chaplygin's ball. It is equal to the system which is obtained by ignoring the 
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equation ()3.1|) for the motion of the point of contact (or the center of gravity), which is the 
same as the R^-reduced system, obtained by working modulo the symmetry group of the 
horizontal translations {A, uj, a) i— > {A, u, a + b), where b G is viewed as a horizontal 
vector in R^. The constants of motion, viewed as functions of [A, u, a) in the phase space 
S0(3) X R'^ X R^, do not depend on the horizontal translations a, and therefore will be 
considered as functions of {A, uj) in the phase space SO (3) x R^ of the rotational motion. 

As observed at after (j2.2(Jj) . the constant of motion j (= the moment of the momentum 
about the point of contact) can be used in order to eliminate u) from the equations of motion. 
In other words, j : SO (3) x R^ R^ is an analytic (rational) fibration, of which each fiber 
is equal to the graph of an analytic (rational) function {u = ujj{A) a function of A), such 
that the projection {A, u) A is an analytic (rational) diffeomorphism from the level set 
of j onto SO (3). 

We note that j is not invariant under the full symmetry group E(2) (= the horizontal 
motion group) of the system. If i? is a rotation around the vertical axis then it acts on the 
phase space by sending {A, u, a) to {RA, u, a). It leaves u = A~^ invariant and we read 
off from ()2.17|) that it sends j to Rj. Therefore the level set is E(2)-invariant if and only if 
the level j is vertical. If j is not vertical and i? is a non-trivial rotation around the vertical 
axis, then the action of R on the phase space sends the level set at the level j to the different 
(disjoint) level set at the different level Rj. 

The level set of all the constants of motion j and T together is diffeomorphic to the level 
set in S0(3) of the function Tj on S0(3), defined by Tj{A) = T{A, u) when AIp^u^^ = j- It 
follows from ()3.7|1 that 

T,(A) = i {Ip^^uj, c.) = i (j, Auj) = i (j, AI,,^-'A-'j). (4.1) 

If j = then, a; = 0, = 0, Tj = 0, and Chaplygin's sphere is at rest. We will exclude 
this rather trivial case in the remainder of our discussions. The vector field of the motion in 
SO (3) is given by (j3.2p . where u = ujj{A) is given by (j2.17p . If j 7^ then ujj{A) 7^ and 
the vector field on S0(3) has no zeros. It may also be observed that replacing j by cj with 
a constant c leads to replacing u by cu, multiplying the vector field on S0(3) by c, whereas 
Tcj = (? Tj. The solutions of the equations of motion are changed only by a rescaling of time 
by the constant factor c. 

If T is a regular value of Tj, then the level set {A G SO (3) | Tj{A) = T} is a smooth 
(algebraic) closed two-dimensional submanifold of S0(3), compact because S0(3) is compact. 
It is oriented by the area form Q/ dTj, where f2 is a volume form on SO (3). Let C be a 
connected component of a regular level set. We conclude that C is a compact connected 
oriented two-dimensional smooth (algebraic) manifold which carries a tangent vector field 
without zeros, which implies that the Euler characteristic of C is equal to zero. According to 
the classification of compact oriented surfaces, this in turn implies that C is diffeomorphic 
to the two-dimensional torus R^/Z^. The considerations below will lead to a much more 
detailed description of the level sets, from which the conclusion that the regular ones consist 
of tori can be obtained without using the just mentioned facts from differential topology. 
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4.2 The Critical Points of the Energy 

For any vector v G let R^(A) G S0(3) be the tangent vector of S0(3) at A G S0(3) 
which is given by 

R,(v4) := A o z/op. (4.2) 

In order to determine the derivative Ri/tUj of the vector- valued function iOj on SO (3) in 
the direction of the vector field R,y, we begin with the observation that ()2.17p implies that 

= Rj, {Alp^u^^) = A (Z/ X Ip^u^j + ^ylp^u^j) . 

Furthermore, 

Ri. M = YLy 63 — v X A^^ 63 = —u X u, 
and therefore it follows from ()2.16|) that 

Rylp^u^j = P {{^j, u X u)u + {Uj, u) u X u) + Ip^u Ru^j- 

Therefore K^, ujj is determined by the equation 

= z/ X ((/ + p) ojj) + p {uj, u X u) u + Ip^u Ri/CUj- (4.3) 

Because Ip^u is symmetric, we have that 

The derivative of Tj, cf. (13. 7|) . in the direction of R^/ therefore is equal to 
Ri^Tj = - ((Ry/p^u) Uj, ujj) + (/p^„ Ri/CUj, Uj) 
= -{{ujj, V X u) U+ {cUj, u) V X u 

-{u X ((/ + p) Uj) + p {Uj, U X U)U, LUj), 

from which we obtain that 

RuTj = {ujj X (/ + p)ujj, u) = {ujj X lujj, u). (4.4) 

Let Sj denote the set of critical points of Tj. It follows from ()4.4|) that A G Sj if and 
only if = ujj{A) satisfies oj x I oj = Q. In view of ()3.4|1 this is equivalent to = 0. Because 
the function Tj is invariant under the motion, the Ri^^-fiow on S0(3), S^- is invariant under 
the motion. Therefore, if ^(0) G Sj then we have for every t that A{t) G S^, which implies 

that ^^ii^ = for every t and therefore uo = uJj{A{t)) is a constant. It follows then from 
()3.2|) that A{t) = A{0) oe*'^ describes a circle, which we will call a critical circle. 

We have u; x Jcj = if and only if I u = lu, which means that u is an eigenvector of the 
moment of inertia tensor /, with eigenvalue l equal to one of the principal inertial moments 
/i, I2 or /s. If this is the case, it then follows from ()2.17|) that 

j = {l + p) Auj - p{Auj, 63) 63, (4.5) 
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Taking the inner product with 63 we obtain that 



]^ = i{Au,e^), (4.6) 
which can be inserted into ()4.5|) in order to yield that 



Au=^ fj + ^ea). (4.7) 



i + P \ L 



From ()4.7|1 we obtain that 



l^ll' 




^+2^ + f^n, (4.^ 



and combining ()4.7p with ()4.ip we obtain that the critical level is equal to 



1 1 (w -112 , pis 



Tcnt = :.^ I|jr + ^ • (4.9) 

Note that the right hand side of ()4.9|1 is a monotonously decreasing function of 6, which 
implies that if the principal inertial moments are taken in increasing order, then the corre- 
sponding critical levels of the kinetic energy appear in decreasing order. 

It follows from ()4.7|) and ()3.H) that the point of contact p{t) moves along a straight line, 
with constant velocity equal to 

| = ^jxe3. (4.10) 
dt i + p 

We note that for the uniformly rolling sphere the axis of rotation HAuj need not be horizontal 
as one might expect. It follows from (j4.6j) that it is horizontal if and only if the vector j is 
horizontal. 

The point of contact is at rest if and only if the moment of momentum is vertical (j is 
equal to a multiple of 63), which according to ()4.7|) corresponds to the case that the vector 
Au is vertical. In this case the sphere is spinning around the vertical axis, which then 
coincides with an inertial axis. 

Let Ej^^ denote the set of G ker(J — t) such that ()4.8|1 holds, and let Sj^^ denote the 
set of A G S0(3) such that ()4.7j) holds for some u G Ej^^. Then is a smooth algebraic 
circle bundle over the smooth algebraic manifold Ej^^ and therefore is an smooth algebraic 
submanifold of SO (3). For different values of the eigenvalue l of /, the sets ^ are disjoint 
and is equal to the union of the t, where l runs over the principal inertial moments. 

If 6 is a simple eigenvalue of / then Ej, consist of two opposite {u ^ —u) eigenvectors of 
I for the eigenvalue l, and Sj ^ consists of two disjoint critical circles. Note that the function 
Tj (and also the equation of motion) is invariant under a transformation A AR, where 
R G SO (3) commutes with I. There exists such R which maps u to —uj and for each such R 
the mapping A^ AR interchanges the two critical circles in In the generic case that 
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all the principal inertial moments Ji, I2, I3 are different, we obtain six critical circles, two 
for every choice of t = Ii, I2, h- 

The case of two equal principal moments of inertia Ii = I2 ^ I3 is that of a body of 
revolution with surface equal to a sphere and center of mass at the center of the sphere. This 
is the example of Routh's sphere with center of mass at the center of the sphere, or Bobylev's 
sphere, which will be discussed in some more detail in Sectional If l = Ii = I2 I3 then 
Ej^^ is a circle in ker(/ — i) and ^ is a two-dimensional torus in S0(3). As discussed before, 

/3 consists of two critical circles. 

If all the principal moments of inertia are equal, l = Ii = I2 = I3, or equivalently / is 
equal to t times the identity, then = S0(3), the function Tj is constant, and all solutions 
of the equations of motion are of the form A{t) = A{0) o e*'^ with a constant vector u. 

We will now verify that each Sj^^ is a is a nondegenerate critical manifold of Tj in the 
sense of Bott j7j, which means that for each A G the null space of the Hessian Tj'{A) of Tj 
at A is equal to the tangent space T^ E^^^ of Ej^^ at A. Because always T^ Ej^^ C ker Tj'(y4), 
we only need to verify that the dimension of the null space of the Hessian is at most equal 
to the dimension of the critical submanifold. This follows from 

Lemma 4.1 Let Iujj{A) = Lujj{A) and K^{A) G kerTj'(yl). If l is a simple eigenvalue of 
I, then u is a multiple ofujj{A). If l is a double eigenvalue of I and a is a nonzero vector 
which is orthogonal to ker(/ — l), then u is a linear combination of u)j{A) and a. 

Proof It follows from (jOj) that R^{A) G kerTj'(A) if and only if, at the point A G SO (3), 

= Rj, (/i^j X Uj) = I ujj X ujj + 1 ujj X Rj, Uj = (/ R.^ Uj — l Ki, ujj) x uj, 

or (/ — l) Ry ojj is equal to a multiple of Uj. This implies that {I — l)'^ R^, = 0, which in 
turn implies that (/ — l) Rj, ujj = 0, because J — i is a diagonal matrix. 
Substituting luj = LUj and / R^, Uj = t R^, Uj in ()4.3p . we obtain 

= u X ({l + p) ujj) + p {uj, u X u) u+ {t + p) RuUJj — p {Ru Uj, u x u) u 
= {1 + p)9 - p{9, u) u, 

in which 9 := v x uoj + R^ uOj. It follows that ^ = 0. 

If i is a simple eigenvalue of I then / R^ ojj = lR^ uj implies that R^, Uj is equal to a 
multiple of 00 j and it follows from ^ = that vxujj =0 and R^ cu^ = 0, which in turn implies 
that z/ is a multiple of Uj. 

If ker(/ — l) is two-dimensional, then ^ = implies that = (z/j, a) = (z/, Uj x a). 
Because ujj G ker(/ — i), we have that a is orthogonal to Uj and it follows that z/ is a linear 
combination of ujj{A) and a. □ 

Because Tj is a continuous function on the compact set SO (3), it attains its maximum 
and its minimum. The variational principle says that the points A where Tj attains its 
maximum (minimum) are critical points for Tj. Therefore, if Ii < I2 ^ h, then the maximum 
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(minimum) value of Tj is equal to the right hand side of ()4.9p . with t = Ii {l = J3). If Ii < I2 
{I2 < I3), then Ej consists of two critical circles and on a transversal two-dimensional 
manifold the function Tj has a nondegenerate maximum (minimum), the nearby level sets 
of which are small loops around the critical points. It follows that the level sets of Tj near 
consist of narrow tubes around the critical circles. This implies that the critical circles 
corresponding to i = Ji < /2 < /s or to i = Is > /2 > h are stable periodic solutions of the 
system. 

Now assume that Ii < I2 < Is- A Morse theoretic argument then yields that Tj has 
cannot have a local maximum or minimum at Sj./j- Note that the index of Tj', the number 
of negative eigenvalues of Tj', is constant along each of the two critical circles in /j. 
If -R = diag(— 1, —1, 1) or i? = diag(l, —1, —1), then A AR leaves Tj invariant and 
interchanges the two critical circles in Sj which implies that the index of Tj' is constant 
along Suppose that it is equal to two, which means that Tj has a local maximum at 

Ej_/2- Let gradTj denote the gradient vector field of Tj with respect to a given Riemannian 
structure on SO (3). Let 7* denote the flow of gradT,-. Define 5*1 and S2 as the set of 
A G S0(3) such that, when t ^ 00, 7*(A) converges to Sj^/^ and T,jj^, respectively. Si and 
5*2 are nonvoid disjoint open subsets of SO (3), with union equal to SO (3) \ 7-3. The set 
SO (3) \ /3 is connected, because S^- 7-3 is a codimension two submanifold of the connected 
manifold SO (3). This leads to a contradiction. In a similar way the assumption that Tj has 
a local maximum at ^j,i2 leads to a contradiction and the conclusion is that the index of Tj' 
is equal to one along Hj^ , which means that transversally Tj has a saddle point behaviour. 

Let P denote the linearization of the return map to a transversal plane (the Poincare map) 
of the flow along the critical circles in Sj^/j. It follows from Corollary 17.21 that det P = 1. 
This implies that | trace P| < 2 if and only if P is conjugate to a nontrivial rotation, whereas 
I trace P| > 2 if and only if P is a hyperbolic map with real eigenvalues A, 1/A such that 
A 7^ ±1, in which case the critical circles in Sj^/j are linearly unstable. The saddle point 
behaviour of Tj near (when Ii < I2 < I3) excludes that P is conjugate to a nontrivial 
rotation. The number trace P depends in a real analytic fashion on j G \ {0}. As we will 
see at the end of Subsection^lbelow, if j is vertical then the critical circles in are linearly 
unstable, which implies that | trace P| > 2 when j is vertical. It follows that the set N of 
j G \ {0} such that trace P = ±2 is a proper closed analytic subvariety of j G R^ \ {0}. 
Using the invariance of the equations of motion under the action {A, u) 1— {RA, uj), j ^-^ Rj 
of the rotations R around the vertical axis and the homogeneity {A, u) 1-^ {A, cu), j 1-^ cj, 
it follows that N, if not empty, is equal to the union of finitely many cones in R^ which 
are invariant under the rotations around the vertical axis. For all j in the complement of 
A^, which is an open and dense subset of R^ \ {0}, the critical circles in T^jj^ are linearly 
unstable. This in turn implies that every critical circle with I uj = I2UJ is unstable with repect 
to the flow in the full phase space. 

Conversely, at any critical circle which is linearly unstable the invariant function Tj must 
have transversal saddle point behaviour. Because the linearly unstable critical circles are 
dense, this leads to a proof that the index of Tj is equal to one along S^^/g without using 
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Morse theory. 

Question 4.2 When Ii < I2 < I3, is every critical circle with I u = I2UJ linearly unsta- 



4.3 The Moment Mapping 

It is also instructive to consider (j, T) as a mapping from the phase space SO (3) x R'^ of the 
rotational motion to x R. The set of singular points of (j, T) , the set of points where 
the rank of the tangent map is less than four, is equal to the union of S0(3) x {0} and the 
set of all {A, cu) where j G R^ \ {0}, A e T.j and to = ujj{A). If Ii < I2 < I3 then the latter 
set is a smooth conic closed submanifold of codimension two in S0(3) x (R^ \ {0}). 

Because of the scaling {A, uj) 1— > {A, cur), which maps (j, T) to (cj, c^T), and because 
the case T = 0, when a; = and everything is at rest, is not very interesting, we restrict 
ourselves to an energy hypersurface where T is equal to a positive constant. (An isotopy 
argument as below shows that the energy hypersurface is diffeomorphic to the Cartesian 
product of S0(3) with a two-dimensional sphere.) The moment of the momentum around 
the point of contact then defines a mapping jx from the energy hypersurface to R^. The 
singular points of Jt are the above singular points on the energy hypersurface. 

According to ()4.9p . the set of singular values of jr consists of the points j G R'^ such that 



in which t = Ii, l = I2 or t = I3. If Ii < I2 < I3, then this set is equal to the union of 
three disjoint ellipsoids with center at the origin and which are invariant under the rotations 
around the vertical axis. The inner and the outer one correspond to t = Ii and 1 = 1^, 
respectively. The points in the phase space of the rotational motion which by jx are mapped 
to these ellipsoids correspond to the extremal critical circles (the stable ones) of the functions 
Tj, and it follows that the inner and outer ellipsoids together form the boundary of the image 
of Jt- Therefore, the image of Jt is equal to the set of all j e R^ such that ()4.11|) holds for 
some L E [Ji, 

The unstable critical circles are mapped to the intermediate (interior) ellipsoid described 
by (j4.11|) with l = I2. For the singular values j in this interior ellipsoid, the level sets are 
two-dimensional, with a singularity of normal crossing type along the two unstable critical 
circles in the level set. 

For the regular values of Jt, the points j such that ()4.11|1 holds for some l such that 
Ii < L < I2 or I2 < L < I3, the level sets are smooth (alegbraic) two-dimensional compact 
oriented submanifolds of the phase space for the rotational motion. 

4.4 Isotopy of the Fibration 

Following Cushman p. 412], a smooth function on a compact manifold, for which the 
critical set consists of nondegenerate critical manifolds (possibly with varying dimensions). 



ble? 








(4.11) 
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will be called a Bott-Morse function. We will use an isotopy lemma for families of Bott-Morse 
functions, which should be well-known. However, because we did not find a reference in the 
literature, we include a proof. 

Lemma 4.3 Let M be a compact smooth manifold and f^ a familie of smooth functions on 
M, depending smoothly on a real parameter e. Furthermore assume that the set of critical 
points of fe consists of finitely many disjoint compact connected smooth submanifolds Ce,i, 
1 < i < N , depending smoothly on e and such that C^^i is a nondegenerate critical manifold of 
fe,i- On Ce,i the function f^ is constant, let F^^i he the value of f^ on C^^i- We finally assume 
that the ordering of the real numbers F^ i, I < i < N , does not change with varying e. Under 
these assumptions there exist smooth diffeomorphisms ip^ and $e of R and M respectively, 
such that tp^o f^o $^ does not depend on e. The ip^ can be chosen to be order-preserving. 

Proof The assumption that the ordering of the critical values does not change implies that 
there exists a family of order-preserving smooth diffeomorphisms ip^ of R such that the real 
numbers Gi := ip^ 1 < i < do not depend on e. It is also quite easy to prove that 

there exists a smooth family of diffeomorphisms of M, depending smoothly on e, such 
that the manifolds Di := \E'7^ {Ce,i) do not depend on e. The functions '■= ipe o fe ° 
have the same properties as f^, but now with the constant nondegenerate critical manifolds 
Di on which g^ has the constant critical values Gi. 

We now follow the idea of the proof of Moser [21]. The condition for a smooth family of 
diffeomorphisms of M, depending smoothly on e, that g^oE^ = E* does not depend on 
e, is equivalent to the condition that = -^E* g^ = E* (X^gi, + ^) , or 



If the equation ()4.12|1 for can be solved locally near every point of M, then a global 
solution can be obtained by means of a smooth partition of unity. 

Near a noncritical point of we can use g^ as one of the local coordinates and ()4.12j) 
then amounts to prescribing the corresponding component of the vector field X^. 

If G Di, then one can introduce a local coordinate system near x^^^ in which x^^^ = 
and X ^ Di corresponds to Xj = for 1 < j < c, if c denotes the codimension of Cj in M. 
Writing x = {y, z) with y G R'^, z G R'^, we view geiy-, z) as a family of functions of y, with 
e and z as parameters. A second order Taylor expansion with respect to ?/ at ?/ = 0, in which 
the remainder term in intergral form is absorbed into the second order term, yields that 



where Q = Qeiv, z) is a nondegenerate symmetric matrix, depending smoothly on all the 
variables. (This is also one of the steps in the proof of the Morse lemma with parameters 




(4.12) 



(4.13) 



g,{y, z) =Gi + - {Qe{y, z) y, y), 
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of Hormander [201 Lemma 3.2.3].) If we take = {Y^, 0) with G R'^, then the equation 
(j4.12p for is equivalent to the equation 

my, Y., ,) + 1 ( Y.) y. ,) + i y. y) = 

for y^, which is satisfied if 

For sufficiently small y the latter equation has a unique solution Y^^ which depends smoothly 
on y, z, and e. 

Piecing together the local solutions by means of a smooth partition of unity, we obtain 
a smooth vector field X^ on M, depending smoothly on e, such that (j4.12p holds. Define 
as the solution of the e-dependent ordinary differential equation ()4.13|) . with initial 
condition So(x) = x. Using the compactness of M we obtain that the are globally defined 
smooth diffeomorphisms of M, depending smoothly on e. Reading the paragraph preceding 
fl4.13|l backwards, we obtain that o does not depend on e. This proves the lemma with 

In order to emphasize the dependence on p of the kinetic energy function on SO (3), we 
now write Tpj instead of Tj. Applying Lemma [4.31 to = T^j, we obtain that there exists 
an order-preserving smooth diffeomorphism ipofH and a smooth diffeomorphism $ of SO (3) 
such that ip oTpj o ^ = Tqj. 

We can make ip and $ to depend smoothly on j when j varies over the unit sphere in 
R'^. Extending the transformations by homogeneity for the scaling {A, uj) i— > [A, cuj), one 
obtains a diffeomorphism $ of S0(3) x (R^ \ {0}), and a diffeomorphism \1/ of (R'^ \ {0}) x R 
of the form (j, T) i— > (j, ip{j, T)), such that (cj, c^T) = c^ip{i, T) for every c > 0, such 
that \l/ o (jp, Tp) o $ = (jq. To). This implies that the smooth diffeomorphism $ maps the 
whole fibration, together with its singularities, of the phase space defined by the constants 
of motion for p = to the one for our given value of p. 

If p = 0, then Ip u = Io,u = I- The equation of motion (|3.4|) then turns into Euler's 
equation of motion 

d 

—Iu = Iujxuj (4.14) 
dt 

for the Euler top, and the constants of motion j of (|2.17|) and T of (|4.H) are given by the 
familiar formulas 

j = AIuj (4.15) 

and 

= 1(1^,^) = I r'A-'j) (4.16) 

for the moment of momentum around the center of mass and the kinetic energy of the Euler 
top, respectively. For details about the Euler top, we refer to Cushman and Bates Ch. 
III]. 
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It follows from ()4.16p that the kinetic energy Tqj of the Euler top is invariant under the 
circle action A t-^ e*-'°p oA, the orbits of which are the fibers of the mapping vj : A A^^ j. 
Note that Vj : S0(3) — *• S||j|| is a smooth fibration of S0(3) over the Euler sphere the 
sphere in with center at the origin and radius equal to ||j||. On the Euler sphere, the 
kinetic energy is equal to the restriction to S||j|| of the quadratic form v ^{v, I~^v), 
defined by the positive definite symmetric matrix 

The critical levels of Tqj are equal to | where l = Ii, I2 or I3, and the regular values 
are the numbers in between the critical levels. On the Euler sphere each regular level set 
has two connected components, opposite to each other, each of which is a smooth closed 
curve, diffeomorphic to a circle. The preimages of these under the mapping vj are circle 
bundles over these circles and therefore each regular level set in SO (3) has two connected 
components, each of which is diffeomorphic to the two-dimensional torus. 

For each extremal level the level sets consists of two critical circles surrounded by narrow 
tubes. For the intermediate critical level, the level set on S||j|| consists of two opposite 
critical points. The complement of these in the level set has four connected component, 
each of which is a smooth curve running from one of the critical points to the opposite one. 
It follows that the level set in SO (3) of the intermediate critical level contains two critical 
circles, the complement of which has four connected components each of which is a smooth 
cylinder running from one of the critical circles to the other. 

For more details about the fibration in SO (3) x for the Euler top we refer to Cushman 
and Bates fl2\ Ch. Ill, Sec. 5]. The point of the isotopy lemma is that there exists a 
diffeomorphism $ which sends the whole fibration with singularities for the Euler top to 
the one for Chaplygin's sphere for an arbitrary value of p. In particular all the qualitative 
statements about the level sets remain true. The regular level sets have two connected 
components, each of which is diffeomorphic to the two-dimensional torus. The level set of 
an intermediate critical level contains two critical circles, the complement of which has four 
connected components each of which is a smooth cylinder running from one of the critical 
circles to the other. 

4.5 Chaplygin 

In the beginning of jHl §6] Chaplygin gave a short description of the critical circles, but 
without relating these solutions to the points where the derivatives of the constants of motion 
are linearly dependent. He also stated that the ones corresponding to the extremal moments 
of inertia are stable and the ones corresponding to the intermediate moment of inertia are 
unstable, but without any proof. 

The question of the smoothness of the level surface of the constant of motion, which 
is related to the question of the linear independence of their derivatives, does not occur in 
Chaplygin 
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5 When the Moment is Vertical 



Next to the critical circles, the solutions with vertical moment j of the momentum around 
the point of contact, j = j^e^, form another interesting special family. 

Note that the condition that j is vertical defines a smooth codimension two algebraic 
submanifold of S0(3) x (R^ \ {0}) x R^. If Ii < I2 < I3 then also the set of singular points 
of the constants of motion, corresponding to the critical circles, is a smooth codimension two 
algebraic submanifold of S0(3) x (R^ \ {0}) x R^. The intersection of these submanifolds of 
special motions consists of the rotations of the sphere around a vertical axis which is equal 
to an axis of inertia, during which the point of contact is at rest. These motions define a 
submanifold of codimension four in S0(3) x (R'^ \ {0}) x R^. 

5.1 Invariance under Rotations about the Vertical Axis 

If j is vertical, then j is invariant under the group of rotations around the vertical axis, 
which is the action of the E(2)-symmetry group, the horizontal motion group, on j. This 
implies that the level set of j is invariant under the transformations {A, u) ^ {RA, u) with 
R G SO (2). The orbits of this action are equal to the fibers of the projection {A, u) {u, u), 
u = A~^ 63, where the space S xR^ of the (m, tu) G R'^ x R'^ such that \u\ = 1 is viewed as 
the phase space for the E(2)-reduced system. Because T is E(2)-invariant, it follows that the 
function Tj on S0(3), which represents T on the j-level set, is S0(2)-invariant, a fact which 
can also be deduced directly from 1)4.11) . In particular the two-dimensional regular level sets of 
Tj are S0(2)-invariant, and therefore are mapped by the projection A 1— >• (m, a;), m = A~^ 63, 
io = ujj{A), to smooth compact one-dimensional algebraic submanifolds of S xR^. Each 
connected component of the regular level set in the {u, C(j)-space therefore will be a closed 
curve, which implies that the motion in the E(2)-reduced phase space is periodic. 

As a consequence, we can apply the reconstruction technique in Hermans Sec. 3.2] in 
order to obtain information about the flow in the full {A, u, a)-phase space S0(3) x R^ x R^ 
of the rolling body. There it is assumed that the symmetry group is compact, but the only 
thing which is needed is that the centralizer in the group of the shift element is a torus. 
Now the centralizer of the element {B, b) G E(2) ~ SO (2) x R^ is a circle subgroup of E(2) 
when B I, or B is a nontrivial rotation around the vertical axis, whereas it is equal to the 
translation subgroup R^ if i? = /. Because the i?(2)-reduced phase space S^ xR^ is equal 
to the S0(2)-reduced phase space of the R^-reduced phase space S0(3) x R'^, we obtain the 
following conclusion. 

Proposition 5.1 Let j be vertical. Each solution on a regular level set in the {u, uj)-space, 
the E{2)-reduced phase space, is periodic. The corresponding rotational motion in the {A, u)- 
-space is quasiperiodic on an analytic two-dimensional torus, depending analytically on the 
parameters j andT. 

For each periodic solution in the {u, u) -space such that the corresponding solution in the 
{A, u) -space is not periodic with the same period, the motion in the full phase space SO (3) x 
R^ X R^ is quasiperiodic on an analytic two-dimensional torus, depending analytically on j 
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and T. In particular the point of contact p{t) remains in a bounded subset of the plane in 
this case. 

If the rotational motion, the motion in the {A, u) -space, is periodic with the same minimal 
period as the motion in the {u, u) -space, then the translational motion, the motion of the 
point of contact pit), is equal to the superposition of a straight line motion with constant 
speed and a periodic motion with the same period as the motion in the [u, uj) -space. 

Some more information about the periodic solutions mentioned in Proposition 15.11 is given 
in Proposition 15.21 below. Proposition 15. HI implies that for at least half of these periodic 
rotational motions the motion of the point of contact p{t) is actually periodic, with the same 
period. 

For more details on the behaviour of the sphere when the moment is vertical, see Kilin 
[221 Sec. 3.3], which also contains computer pictures of orbits of the point of contact. 



5.2 Fourier Series 

In order to appreciate the statements in Proposition 15.11 we recall what it means that 
the motion in a suitable open subset M of the phase space is quasiperiodic on analytic r- 
-dimensional tori, with analytic dependence on s parameters, varying in some open subset 
E of R*. It means that there exists an analytic diffeomorphism $ from (R^/Z*^) x to M 
and an analytic function i/ : E ^ W , such that the pull-back w = ^* v of the velocity field 
f on M is of the form w{x, e) = (z^(e), 0), x G H'" /Z^, e G -E. This implies that the solution 
curves 7 in M are of the form 7(t) = $ (xq + 1 z^(e), e), t G R. In our case r = 2. 

If we apply this to the rotational motion, then it follows that the right hand side r {Auj)x 
63 in ()3.1|) is of the form p [xq + t z^(e), e), where p is an analytic mapping from (R^/Z^) x E 
to the horizontal plane R^. Using Fourier expansion it follows that we can write 

p{xo + tu{e), e)=Y^ Cfc(e) e^-^^o+t^^U)^ (51) 

in which the Fourier coefficients Cfc(e) depend analytically on the parameters e. The ana- 
lyticity of the function p implies that the Fourier coefficients Cfc(e) are rapidly decreasing as 
\\k\\ 00. Formal termwise integration of ()5.H) would lead to 

p{t) = p{0)+t J2 Cfc(e) e2"^<^°''=> (5.2) 

I ^fc(^) 27Ti(xo+tu(e),k) /gox 

^ ^ 2ni{u(e), k) ' ^ ' ' 

The coefficient of the linear term in t, the secular term, contains, apart from the term 
Co(e) (which is equal to the average of the function p over the torus), terms for nonzero k if 
and only if the components z^i(e) and z/2(e) have a rational ratio. If we assume that 1^2 (e) 7^ 
and z/i(e)/z/2(e) = ni/n2, with Ui G Z, ^2 G Z>o and gcd (ni, ^2) = 1, then 

t (z/(e), k) =t ^^^^ (ni ki + n2 ^2) , 
n2 
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which is equal to an integer if t is equal to an integral multiple of T := n2/'^2(e)- This in 
turn implies that the function t p{xo + tu{e), e) is periodic with period equal to T. For 
generic analytic functions z^i(e), z>'2(e) the ratio z/i(e)/z/2(e) is rational for a dense subset of 
parameter values, where the denominator n2 and therefore also the period T is unbounded 
in every nonvoid open subset of the parameter space. The coefficients of the secular term in 
general would have a correspondingly wild behaviour as a function of the parameters. 
In the case that the fraction z/i(e)/z/2(e) is irrational, when the function 

p{xo + tv{e), e) 

is not periodic, then the sum in (j5.3p over the k & 7? such that (z/(e), fc) 7^ need not 
converge, due to the possibility that the denominators (z^(e), k) may become arbitrarily 
small. This problem can arise, despite the rapid decrease of the Fourier coefficients Cfc(e) 
when ||A;|| 00. If the fraction Vi{t)/v2{^) satisfies suitable diophantine inequalities, then 
the sum is convergent and defines a quasiperiodic function of t on a two-dimensional torus. 

Proposition 15. II implies that for the motion of Chaplygin's sphere when j is vertical none 
of the above complications occur: no wild behaviour of the coefficients of the secular term 
and no problem with convergence of the Fourier series in ()5.3|) . Note that the conclusions of 
Proposition 15. II have not been obtained by means of an analysis with Fourier series, but by 
using the reconstruction technique in Hermans fHl Sec. 3.2] instead. 



5.3 Euler's Equations 

In this subsection we give some explicit formulas, which among other show that if j is vertical, 
then the rotational motion is determined by Euler's equations. An extensive discussion of 
the Euler top can be found in Cushman and Bates ^21 Ch. III]. 
If j = is 63 , then it follows from (j4.ip that 

2T={Au:,3)= js {A cu, eg) = J3 {to, u) (5.4) 

Because (j3.3|) implies that (cj, %) = (cj, u x cu) = 0, we obtain from (j5.4j) that 



( — , u) = —(io, u) = 0. 



But then fj3.4p implies that 



{I + p)uj X u = I u X uj = —{I + p)lj, (5.5) 

dt 

which are Euler's equations with I replaced by / + p. 

The equation 63 = j = A Ip^^^ yields in view of (|5.4|) that 

2T p 

j-iU= {I + p)uj - p{u, u)u= {I + p)uj ; — M, 

J3 
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which leads to the formula 

u = .2 i\rr. {I + P) (5-6) 

which expresses u in terms of u. The differential equation ^ = m x u; then leads to the 
differential equation 

du j^ + 2Tp 

— = u X J u, (5.7) 

for u only, which also is an equation of Euler type. Here J = (/ + p)~^, cf. fj2.18p . 
From ()5.6|) and ()5.4j) we obtain that 

2T 
J3 +2Tp 

which shows that the solutions u{t) run on the intersection of the sphere {u, u) = 1 with 
an ellipsoid defined by the symmetric matrix J. Using ()5.6p it follows that also (/ + p)uj{t) 
runs over the intersection of a sphere with an ellipsoid, explicitly given by 

{{I + p)uj,{I + p)uj) = (^^^-J^) , (5.9) 
{uj,{I + p)uj) = 2T {l + 2Tp/js'). (5.10) 
The kinetic energy T is given in terms of u by 

2T= f^^/;\ , (5.11) 

1 — p{u, J u) 

cf. ()5.8|) . The critical points of T on the unit w-sphere are the unit eigenvectors u of J, with 
the eigenvalues l/(i + p), with l = Ii, I2 or I3. The corresponding critical value is equal to 

Tent = (5.12) 

which is in accordance with ()4.9p when j = J3 63. It is well-known that the eigenvectors for 
the two extremal eigenvalues of J + p are stable equilibrium points of the Euler equation, 
whereas the eigenvectors for the intermediate eigenvalue are linearly unstable equilibrium 
points of Euler's equation (|5.7|) . See for instance Cushman and Bates jT21 p. 117]. The 
latter implies that, when j is vertical, the critical circles of Tj for the intermediate critical 
values are linearly unstable. 



5.4 The Translational Motion 

Let u = u{t) be the solution of (j5.7|) on the intersection curve of the ellipsoid (j5.8|) with 
the unit sphere U. The rotational velocity vector uj = uj(t) is determined in terms of u by 
means of ()5.6|) . The projection A \^ u = 63, cf. fl2.12|) . exhibits the Tj-level set in S0(3) 
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as a principal S0(2)-bundle over the aforementioned curve on U, where SO (2) is the group 
of rotations about the vertical axis, acting on SO (3) by means of left multiplications. The 
projection intertwines the flow on the T^-level set in SO (3) with the motion on the curve on 
U determined by ()5.7|) . where the flow on the Tj-level set in SO (3) is given by the differential 
equation ()3.2|) in which u = u{t) is determined by ()5.6|) . 

If Tj is close to the intermediate critical level, then the solution u = u{t) of ()5.7|) will 
stay for a long time near one of the unstable equilibria Ucnt of (jS.Tj) before it moves on to 
the other one. During this time uj = uj{t) will stay close to the nonzero vector ujcnt which 
is determined by ()5.6p with u replaced by Merit- It follows that A = A{t) will make many 
rotations during that time. This leads to the following conclusion. 

Proposition 5.2 LetTper denote the levels ofTj such that the motion on the corresponding 
torus in SO (3) is periodic with the same period the motion in the [u, uj) -space, cf. Proposition 
\5.1l Then T^er is an infinite subset o/R with the intermediate critical level T = j^^ /2I2 as 
its only accumulation point. This accumulation point is approached by T^er both from above 
and from below. 

The translational motion, the motion of the point of contact p, is obtained by integrating 
the right hand side r [Auj) x 63 of ()3.1|) . Here the vector 9 := Auj., the rotational velocity 
vector in space coordinates, the herpolhode in Poinsot's description of the Euler top, lies in 
view of (j5.4j) in the fixed horizontal plane 63) = 2T / j^. We will use the discussion in 
Cushman and Bates ^1 II. 7. 2] of the construction of Poinsot. (Additional information can 
be found in Routh Art. 151, p. 98 and pp. 471-473]. An interesting fact is for instance 
that the herpolhode is always concave towards the interior, without inflexion points. In his 
Theorie nouvelle de la rotation des corps, 1834, Poinsot drew the herpolhode like a snake (= 
herpes in Greek), which therefore is misleading.) 

To begin with, the image of the T,-level set in S0(3) under the projection A ^ 9 = Auj 
is invariant under the rotations around the vertical axis, and is therefore known if we know 
how (6*, 9) = {uj, uj) varies as uj is coupled to u by means of ()5.4|) and u E U runs over the 
curve determined by (j5.8|) . 

We will restrict ourselves to one of the two opposite connected components of the intersec- 
tion of U with the ellipsoid ()5.8p . which is the orbit of the motion on the w-sphere. On it, the 
function {9, 9) has four critical points which all are nondegenerate. Two of these correspond 
to the maximal value -R^ax two to the minimal value -Rmin; where -Rmax > -Rmin > 2T/j3. 
To be more precise, if Tcrit,i = jlli denote the critical values of Tj, cf. (|5.12|) . then we 
have the following two cases. 

i) Tcrit,3 < T < TCTit,2- We have M3 7^ on the intersection of U with the ellipsoid ()5.8|) . 
and therefore M3 has one sign on the w-orbit. The critical points corresponding to the 
maximal value are equal to the two intersection points (±Mi, 0, ^3) of the orbit with 
the coordinate plane U2 = 0. The critical points corresponding to the minimal value 
are equal to the two intersection points (0, ±M2, M3) of the orbit with the coordinate 
plane ui = 0. 
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ii) Tci.it,2 <T < Tcrit, 1- We have Mi 7^ on the intersection of U with the elhpsoid ()5.8|) . 
and therefore Ui has one sign on the w-orbit. The critical points corresponding to the 
maximal value are equal to the two intersection points {ui, 0, ±^3) of the orbit with 
the coordinate plane U2 = 0. The critical points corresponding to the minimal value 
are equal to the two intersection points {ui, ±U2, 0) of the orbit with the coordinate 
plane U3 = 0. 

It follows that the image of the Tj-level set in SO (3) under the projection A\—>-p = rAujxes 
is equal to a circular annulus in the plane, with center at the origin. When u{t) runs around 
the orbit in U then it subsequently passes an intersection point with the coordinate plane 
U2 = 0, then an intersection point with a second coordinate plane, then the other intersection 
point with the coordinate plane U2 = 0, and finally the other intersection point with the 
second coordinate plane, before it closes. The corresponding point p in the annulus will then 
reach the outer circle, with a second order contact, then touch the inner circle, return to 
the outer circle at a point which is rotated over an angle a as compared to the first contact 
point with the outer circle, and then touch the inner circle for the second time before the 
curve in U closes. 

If the rotational motion A{t) is periodic with the same period as the motion u{t) on U, 
then the third point of contact with the outer circle is equal to the first one, which means 
that 2a is equal to an integral multiple of 27r. There are two cases. 

a) a itself is not an integral multiple of 2tt, which correspond to the case that the second 
point of contact with the outer circle lies opposite to the first one. In this case the 
p-orbit is symmetric about the origin, its time average is equal to zero and the motion 
of the point of contact p{t) with the horizontal plane is periodic. In other words, the 
speed of the straight line motion in Proposition 15.11 is equal to zero. 

b) a is equal to an integral multiple of 27r, which correspond to the case that the second 
point of contact with the outer circle is equal to the first one. Equivalently, q{t) is 
periodic with a period equal to half the period of the motion u{t) on U. In this case 
the p-orbit is not symmetric about the origin. 

The argument preceding Proposition 15.21 yields that a is given by a smooth real-valued 
function of T which tends to ±00 as T — > Tcrit,2- 

Moreover, if T converges to the intermediate critical level T = ^3^/2/2, then the inner 
boundary circle of the p-annulus shrinks to the origin. Because u{t) stays for a long time 
near the critical point ±62, to which the intersection point of the w-orbit with the coordinate 
plane ui = in case i) and M3 = in case ii) is close, the conclusion is that p{t) stays for a 
long time close to the inner boundary circle, running many times around it in the process. It 
follows that the time average of p{t), which is equal to the speed of the straight line motion 
in Proposition 15.11 converges to zero when T converges to the intermediate critical level 

T = J3V2/2. 

We therefore arrive at the following conclusions. 
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Proposition 5.3 Let T^^^ denote the levels of Tj such that the motion of 

q(t) = rA{t) uj{t) X 63 

is periodic with half the period of the motion in the {u, u) -space. Then T^^^ is an infinite 
subset of H with the intermediate critical level T = /2I2 as its only accumulation point. 
This accumulation point is approached by T^^^ both from above and from below. If the speed of 
the straight line motion in Proposition \5.1\ is nonzero, then necessarily T G T^cv & T^^^ 
converges to the intermediate critical level T = j3^/2/2, then the speed of the corresponding 
straight line motion, cf. Proposition \5. 1\ converges to zero. 

We conjecture that for most values of T G T^^^ the speed of the straight hne motion in 
Proposition 15.11 is not equal to zero. Here the word "most" can mean all except finitely 
many, or for generic values of /i, h- 

If T ^ T^gj, approaches a value Tq G T^'^j. for which speed of the straight line motion in 
Prop osit ion 15 ■ 1 1 is not equal to zero, then the bounded area in which the quasiperiodic motion 
of p{t) takes place "opens up to infinity" , and closes again to a bounded subset when T has 
passed the value Tq. If the above conjecture holds, then this scenario takes place infinitely 
often when T approaches the intermediate critical level T = j3^/2/2, but with the average 
speed of the point of contact converging to zero. 

5.5 Chaplygin 

The remainder of Chaplygin jHl §6], starting with the sentence "In addition, there is an 
exceptional case, ■ ■ ■" , consists of a discussion of the case that j is vertical. This discussion 
contains several interesting observations, but it does not give the qualitative information 
about the motion as in Proposition 15.11 

6 Bobylev's Sphere 

A rigid body is called a solid of revolution if it is dynamically symmetric with respect to all 
rotations R about a given axis, for which we can take the vertical axis. This means that 
both the surface S and the inertial tensor I are invariant under such rotations R. Because 
for Chaplygin's sphere the surface S already is invariant, it is a solid of revolution if and 
only if Ji = I2. In this case Chaplygin's sphere is equal to Routh's sphere with the center of 
mass at the center of the sphere. This case has been studied by Bobylev [S], where in our 
case there is no gyroscope as mentioned in the title of 0. 

The assumption of having a solid of revolution implies that the equations of motion are 
invariant under the right SO (2) -action 

{A, Lu, a) ^^ {AR~^, Ru, Ra) 

of rotations R about the vertical axis. The right SO(2)-action also leaves the kinetic energy 
T and the moment j invariant. The quotients of the connected components of the regular 
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(j, T)-level surfaces by the right SO(2)-action are diffeomorphic to circles, which imphes that 
the solutions of the right S0(2)-reduced system, for the regular levels of (j, T), are periodic. 
The reconstruction method in Sec. 3.2] then leads to the following conclusions. Note 
that the vector Au is invariant under the right SO(2)-action, because {AR~^) {Ru) = Auj. 

Proposition 6.1 Suppose that two of the moments of inertia are equal to each other. Then 
the right SO (2) -reduced rotational motion on the regular levels of {j, T) is periodic, with a 
period which depends analytically on j and T. In particular the vector Auj, and therefore 
also ^ = r (A uj) X 63 performs a periodic motion. 

It follows that the rotational motion on the regular levels of (j, T) is quasiperiodic on 
two-dimensional analytic tori, depending analytically on j and T. Also, the motion of the 
point of contact p is equal to the superposition of a straight line motion with constant speeed 
and a periodic motion with the same period as that of the right S0(2)-reduced motion. 

As in the case when j is vertical, cf. Proposition 15.11 there are no problems with secular 
terms or with the convergence of Fourier series for the motion of the point of contact p. In 
contrast with Proposition 15. !( we obtain here that ^ is always periodic. 

It follows from (HTT^ that 

F = Iix + hyz, X = U1UJ1 + U2UJ2-, y = u^s, z = U3. 

This is known as Jellet's integral, cf. Routh Art. 243]. 

In the same notation, the constant of motion G of ()3.19|) takes the form 

G= (Ji+p)' w + ih + pf y^-2p{x + yz){{h + p) x + ih + p) y z) + p" {x + yzf. 

Together with the kinetic energy we thus obtain three constants of motion in the four- 
dimensional E(2) X S0(2)-reduced phase space, the "fully reduced" phase space. The regular 
level sets of all the constants of motion are algebraic curves, the periodic motion on which 
can be obtained by means of quadratures. We do not go into further details about this here. 

6.1 Chaplygin 

In his Introduction, Chaplygin ^ referred to the papers of Bobylev j3] and Zhukovsky for 
the case that two moments of inertia are equal. In the beginning of P| §6], Chaplygin wrote 
"We will not treat the case when two or all three principal moments of inertia are equal, 
because the motion of such a sphere has already been investigated (see the Introduction)." 
Apparently Chaplygin did not feel that the articles of Bobylev and Zhukovsky, which I 
have not seen, needed further comments. In particular I wonder whether Bobylev and/or 
Zhukovsky used the moment j of the momentum around the point of contact as a constant 
of motion. In the paper to which Chaplygin referred for the fact that j is a constant of 
motion, there is no reference to Bobylev or Zhukovsky. 

The description of Bobylev's paper j3] in the Fortschritte der Mathematik says: "After 
the proposed integration, which can be performed with the help of the elliptic functions of 
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Weierstrass, the author reaches the conclusion that the center of the sphere describes a curve 
which is enclosed between two parallel straight lines and has a periodic character, where it 
successively reaches the one and the other straight line with constant distances between the 
successive contact points on each of the straight lines." This corresponds to the description 
of the point of contact p in Proposition I6.ll 

As observed before, Chaplygin's sphere with two equal moments of inertia is equal to 
Routh's sphere with the center of mass at the center of the sphere. In [HH Art. 243], where 
Routh's sphere is treated, no special attention is paid to the case that the center of mass is 
at the center of the sphere. 

7 An Invariant Volume Form 

We return to the general case of Chaplygin's sphere, with arbitrary moments of inertia, 
arbitrary total mass and radius of the sphere and arbitrary moment of the momentum 
around the point of contact. In the following lemma we use the notation of Subsection 13.31 
Note that X{u) in ()3.15|) is strictly positive when = 1, because the eigenvalues of the 
symmetric matrix J = (/ + p)^^ are equal to 1/ (Jj + p), i = 1, 2, 3, and therefore strictly 
smaller than 1/p. 

Lemma 7.1 Consider the equations of motion ^ = uxLU,^ = vxujin the (m, v)-space 
X R'^, in which uj is determined in terms of u and v by (j3.1H) . Let fl be the volume form 
in which is equal to X(u)~^^'^ times the Euclidean volume form, in which X(u) is defined 
by (j3.15p . Then Q is invariant under the flow in . 

Proof The velocity field of the flow is equal to the vector field 



in which u = uj{u, v) is determined by (|3.11|) . The divergence of this vector field is equal 
to the trace of the derivative, and therefore equal to the sum of the trace of the derivative 
Du of u u X uj{u, v) and the trace of the derivative of v v x uj{u, v). Because 
the traces of the linear mappings 6u ^ 6u x uj{u, v) and 6v ^-^ 6v x uj{u, v) are equal to 
zero, the divergence is equal to trace + trace Z)^, in which : 5u u x ^^^q^ Su and 
D'^ : 6v V X '"^ 6v. From ()3.11|) we obtain that 



in which ipi'^, "v) '■= y{u, v)/X{u) and Y{u, v), X{u) are given by ()3.16|) . ()3.15p . respectively. 
We now use that the trace of Su u x J 6u is equal to 



R^ : {u, v) {u X u, V X u) , 



uj{u, v) = J V + ip{u, v) Ju, 



(7.1) 



3 3 



u X a, Ci) = 0. 



i=l i=l 
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It follows that trace is equal to the trace of the rank one mapping 



and therefore equal to 



because 



dibiii V) 

6u (-^ ( — 6u ] u X J u, 

ou 



dtPiu, v) ^ J = {ux Ju, Jv)lX{u), (7.2) 
ou 

dX{u) 



[u X Ju) = —2{u X Ju, J u) = 0. 



du 

Similarly the trace of 6v ^ v x J 6v is equal to zero, and therefore the trace of is equal 
to the trace of the rank one mapping 

/ d^{u, v) \ 

\ — d / ' 

which is equal to zero because 

{u, J {u X J u)) = {J U, U X J u) = 0. 

The conclusion is therefore that the divergence of the vector field is equal to ()7.2|1 . On 
the other hand the derivative Ro, X of the function X in the direction of the vector field R^^ 
is equal to 

—2{u X uj, Ju) = —2{u X Jv, J u) = 2{u x Ju, Jv), 

and therefore 



div R^ = ^X-i R^ X. 



It follows that 



div j^^j _ Rc.^ + X-^/^ div R^ = 0, 

which completes the proof of the lemma. □ 

Let M be a smooth manifold of dimension m, Vt a. smooth volume form on M and / a 
smooth function on M such that d/ 7^ at every point of the level set := {x G M | 
f{x) = c}. Then Mc is a smooth (m — l)-dimensional submanifold of M, and there is a 
unique volume form u on Mc such that 

Q:r{vi, Vm-1, Vm) = UJ^ {vi, Vm-l) <^fx {Vm) (7.3) 

whenever x G Mc, Vi G T^. Mc every 1 < i < m — 1, and Vm G T^. M. The volume form uo on 
Mc is smooth and nonzero at every x G M where Vt^ 7^ 0. It is called the relative quotient 
of Q and df and denoted hj u = Q/ df. 
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Jacobi observed in [211 10-14. Vorlesung], that if f is a smooth vector field on M such 
that V f = and its divergence C^Q/fl with respect to Q is equal to zero, then the flow of 
V leaves Mc, df and Q invariant, and therefore uj as well. In other words, if Vc denotes the 
restriction of v to Mc, which is tangent to Mc, then the divergence Cy^u/uj of Vc with respect 
to u is equal to zero. 

Applying this principle succesively to the functions in the left hand sides in fl3.9p . which 
are all invariant under the action {u, v) t— > {Ru, Rv) of arbitrary rotations R, and using 
that the set determined by (j3.9|) can be identified with SO (3) if j is not vertical, we arrive 
at the following corollary, where for vertical j we can apply a continuity argument. 

Corollary 7.2 Let dA be a Haar volume form on S0(3), a volume form on S0(3) which 
is invariant under right (or left) multiplications with elements o/SO(3). Then X^^/^ dA is 
invariant under the Ki^^-fiow on SO (3). 

On the regular level surfaces for the kinetic energy function Tj, the area form 

dA/dT,- 

is invariant under the H^.-flow on SO (3). 

If M is a two-dimensional smooth manifold, v a nowhere vanishing smooth vector field 
on M, and a is a nowhere vanishing smooth area form on M which is f -invariant, then the 
fact that the three-form da is equal to zero on M implies that 

= £^ a = d (i^ a) + if, (da) = d (i^ a) , (7.4) 

or that the nowhere vanishing one-form 

/? := if, a 

is closed. Let g{x) be the function which is obtained by integrating (3 along a curve in 
M starting at some base point and ending up at x. (This is an allowed procedure in the 
"integration by quadratures" philosophy.) Then dg = (3, hence v g = and the orbits of the 
f-solution curves correspond to the level sets of g. Note that dg = (3 is nowhere vanishing, 
which implies that the connected components of the level sets are smooth curves. Also note 
that the function g is always globally defined on the universal covering space of M, but that 
on M in general it will be a multi- valued function, with the indeterminacy that g{x) has to 
be replaced by g{x) + ([7], [u]) if the curve ending up at x is followed by a loop 7 which 
starts and ends at x. Here [7] G Hi(M) and [u] G H^(M) denote the homology and (de 
Rham) cohomology class of 7 and uj, respectively. 

If Qq is a nonzero smooth volume form on a manifold M, then every smooth volume 
form n on M is of the form Q = fiQo, for a unique smooth function fi. Therefore the 
search for an invariant volume form is a matter of finding the right factor (multiplier) /i. If 
one has such a multiplier on an n-dimensional manifold and one also has n — 2 independent 
constants of motion /i, . . . , /„_2, then taking successively the relative quotient volume forms 
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on the level manifolds /, = Cj, one obtains multipliers on the (n — j)-dimensional level sets 
fi = Ci, . . . , fj = Cj. For j = n — 2 one finally obtains the "last multiplier" on the two- 
dimensional level surface of all the constants of motion, to which one then can apply the above 
integration by quadratures. This is method for integration of vector fields by quadratures, 
which has been introduced in j^H 10-14. Vorlesung], is called Jacobi's last multiplier method. 

If the two-dimensional M is compact and connected then the fact that v has no zeros im- 
plies that M is diffeomorphic to a torus, as we have observed before at the end of Subsection 
14.11 Siegel [HUl Lemma 3 and 4] proved that there exists a smooth closed loop C in M such 
that V is everywhere transversal to C and that for every such C and every f-solution curve 
7 there exists a t > such that 7(t) G C. The transversality of f to C implies that if a; G C, 
and t I— > 7(t, x) denotes the f-solution curve with 7(0, x) = x, and T(x) is the smallest t > 
such that 7(t, x) G C, then T depends smoothly on x G C and we obtain a smooth Poincare 
map P : X ^ 7(T(x), x) : C ^ C . The restriction f3c of the above one-form /3 = i„ a to C 
is a smooth one-form on C without zeros, and there an angle coordinate ^ on C such that 
d^ = (5c 1 which is unique up to an additive constant. The f-invariance of the area form 
a implies that (3 is f- invariant. In turn this implies that I3c is invariant under P and we 
conclude that P{9) = + c, where c is a constant. In other words, the return map is a 
rotation. 

By modifying the speed of the solution curves before they arrive at C, we can arrange 
that the return time T[x) is a constant. In other words, there exists a strictly positive 
smooth function f on M (which we can choose to be non-constant only in a thin strip at 
one side of C), such that T{x) is equal to a constant if we replace v hj f v. Let wc be the 
unique tangent vector field of C such that i^c «c = 1- Because wc is invariant under P, we 
can carry wc around with the f-fiow and obtain an extension w of wc which is a smooth 
vector field on M and commutes with /f by construction. Is is also clear that w and fv 
are everywhere linearly independent. 

We now recall the argument of Arnol'd and Avez [Bj Appendix 26] that the /f-fiow is 
quasiperiodic. It follows that if e*" denotes the fiow after time t of the vector field v, then 

(t, s)^e'^^oe''" (7.5) 

defines an action of on M. Because of the linearly independence of f v and w, the orbits 
are open subsets of M. because the orbits form a partition of M and M is connected, there 
is only one orbit, equal to M. In other words, the action is transitive. If e*-^" oe*"'(a;) = x 
for some x G M then e*-^" o e*"'(x) = x for every x G M. The period lattice 

n := {{t, s) G R2 I e*^''oe*"' = 1} 

is a discrete additive subgroup of R^, and because for each x G M the mapping {t, s) t— > 
g^/f oe*"'(x) induces a diffeomorphism from R^/II onto the compact manifold M, the con- 
clusion is that the lattice 11 is two-dimensional. In the coordinates with respect to a Z-basis 
of n, R^/n is equal the standard torus R^/Z^. Because in these coordinates the vector 
field f V is constant, we arrive at the conclusion that the /f-fiow is quasiperiodic on a 
two-dimensional torus. 
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It is clear from its introduction that the function / is far from unique. Actually, Kol- 
mogorov [221 proved that if the rotation number of the Poincare map P satisfies suitable 
diophantine inequalities, ensuring that it cannot be approximated too rapidly by means 
of rational numbers, then also the f-flow, without the time-reparametrizing factor /, is 
quasiper iodic. He also showed that for this conclusion the diophantine inequalities for the 
rotation number are essential, in the sense that in general it is not sufficient to assume that 
the rotation number of P is irrational. 

In Corollary 18.41 we will obtain that for the Chaplygin sphere the rotational motion is 
quasiperiodic on two-dimensional tori, if the reparametrization of the time corresponds to 
multiplication of the vector field by the specific function / = where X(u) is given 

by (nrTH|i . 

7.1 Chaplygin 

The last part of Chaplygin §2] , starting with "To solve the problem completely, . . . " , 
contains the proof of Lemma (7.11 followed by the conclusion, in one line, that Jacobi's last 
multiplier method can be applied in order to solve the equations of motion by quadratures. 
Apparently at the time of this method was so well-known, that no further explanations 
or references were needed. 

8 Two Commuting Vector Fields 
8.1 The Second Vector Field 

It follows from ()4.4|) that the tangent spaces of the level surfaces in SO (3) of the function 
Tj are spanned by the vector fields R^^^. and R(i+p)^.. It is therefore natural to investigate 
the divergence of the vector field R(/+p)a;j with respect to the area form dA/ dTj, in 

analogy with Corollary O Note that I + p = J-^, d. (IZTHj) . 

Lemma 8.1 Consider the vector field in the {u, v)-space x R^, defined by Rj^ u : = 
u X u , K,^ V := V X u , where v = v{u, v) := (/ + p) uj{u, v) and uo = uj{u, v) is determined in 
terms of u and v by (|3.1H) . Let X{u) be defined by (|3.15|) and let Q be the volume form in 
R^ which is equal to X{u)^^^'^ times the Euclidean volume form. Then the divergence ofK^ 
with repect to Q is equal to zero. 

It follows that if dA is a Haar volume form on SO (3), then the divergence of the vector 
field R(/+p)tjj on SO (3) with respect to X"^^"^ dA is equal to zero. Also, the divergence is 
equal to zero of the vector field R{/+p)(^j on any regular level surfaces ofTj, with respect to 
the area form X"^^"^ dA/ dTj. 

Proof The proof follows the same lines as the proof of Lemma 17. ![ The calculations are 
actually somewhat easier this time, because the expression 

//(■U, !>) = (/ + p) Uj{u, V) = V + 1p{u, v) U 
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for V is simpler than the formula ()7.H) for io. 

The divergence of R,^ is equal to trace D^^ + trace D^, in which and is equal to the 
trace of the derivative of m x z/(m, = m x w and v x v) = ip{u, v) ■ v x u with respect 
to u and v, respectively. It follows that the divergence is equal to the trace of the rank one 
linear mapping 



( — Sv ] V X u, 

ov 



and therefore equal to 



dv 

On the other hand 



[v X u) = {u, J{v X u))/X{u). 



KyX = —2{u X u, J u) = —2{u X V, Ju) = 2X div R^, 
which in the same way as at the end of the proof of Lemma 17.11 implies that 

div j^^^ _ -^^'^^^ RuX + X-^l^ div R^ = 0, 

or that the divergence of Rj, with respect to is equal to zero. 

The statements about the volume form on SO (3) and the area form on the level sets of 
Tj follow in the same way as CoroUarv 17.21 follows from Lemma f7. II □ 

Lemma 8.2 Let dA he the Euclidean volume form on S0(3). Then 

(dA/dT,)(R^^.,R(,+p)^J = L 



Proof It follows from 14.41 that if i/ = Uj x (I + p) coj, then 

R,^ Tj = {ujj X (I + p) Uj, ujj X (I + p) ujj) = dA {uj, (/ + p) uj, ujj x (I + p) ojj) , 

which implies the statement of the lemma in view of the defining equation (j7.3j) of the relative 
quotient of a volume form and the total derivative of a function. □ 

Proposition 8.3 Let the factor X{u) he defined hy ()3.15|1 . Define the vector fields ^ and 
r/ on SO (3) hy ^ := X(uY^'^ R^^. and i] := X{uy^'^ R(/+p)wj; respectively. Then the vector 
fields ^ and rj commute. 

Define the area form aj and the one-forms (3 and 7 on the regular level surfaces ofTj hy 
aj := Xlu)"^ dA/ dTj, (3 := i^otj and 7 := ijjttj, where dA is the Euclidean volume form 
on S0(3). Then 77) = 1, = /? A 7 and the one-forms (3 and 7 are closed. The area 
form Uj and the one-forms (3 and 7 are invariant under the flow of hoth vector fields ^ and 
V- 
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Proof It follows from ()7.4|) and the fact that the divergence of v with respect to a is equal 
to zero if and only if the one-form i„ a is closed. Because for any functions / on M we have 
that ga = f g iyO, it follows that for every nowhere vanishing smooth function / we 
have that the divergence of v with respect to a is equal to zero if and only if the one-form 
(/ ^ a) is closed. 
If we apply this with v = R^., 

a = dA/dTj = X^/^aj, 

and / = X^/^, then it follows from CoroUarv 17.21 that /? = i^ aj is closed, or equivalently 
that C^aj = 0. In a similar manner it follows from Lemma l8.ll that 7 = 1,, is closed, or 
equivalently C^uj = 0. 

On the other hand Lemma 18.21 implies that t]) = 1, from which it follows in turn 

that 

= C^aj{^,r]) = aj{^, r]]) , (8.1) 
= Cr,aj{^, ri) = aj{[ri, ri) = aj{r], r]]) . (8.2) 

Here we have used in (|8.ip that C^aj = 0, = [^, ^] = 0, and C^r] = rj] = 0, whereas 
in in (18. 2p we have used that = 0, = [r], C,] = — r]], £^77 = [r], rj] = 0, and the 
antisymmetry of aj. It follows from ()8.1|) that [^, t]] is a multiple of ^ and from ()8.1|) that 
[^, 77] is a multiple of rj. Because ^ and r] are everywhere linearly independent on the regular 
level sets of Tj, it follows that [^, rj] = there. Because the regular level sets are dense, it 
follows by continuity that the vector fields C, and rj commute on all of SO (3). □ 



Corollary 8.4 // the rotational motion is parametrized by a time variable t which is related 
to the time i by ^ = X(u{t))~^^'^ , then the rotational motion on the regular level sets is 
quasi-periodic on two-dimensional analytic tori, depending analytically on the parameters j 
and T. 



Proof We have that 

dr dt ' ^ 

The conclusion of the corollary follows from the discussion of the action ()7.5j) with M, v and 
w replaced by a regular level surface, ^ and rj, respectively. □ 



Question 8.5 Is there a proof of Corollary E21 Lemma f8. II and Lemma f8.2[ and therefore 
also of Proposition IHSl which is based on a general principle, in the same way as Proposition 
12. II follows from Noether's principle for nonholonomic systems in Lemma fl.lL '' 
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8.2 A Zero Average 

According to (jH.lll . the time derivative of the j-component (p, j) of the point of contact p 
of the sphere with the horizontal plane is equal to r times the quantity We begin with 

{{Auj) X eg, j) = {Auj, 63 X j) = {uj, {A-^es) x A'' j) 
= {u, u X v) = {u, u X I u) = — det {u, u, I oj) . 

Here we have used ()2.12|) and ()3.8p in the third equation, and in the fourth equation we have 
used fl3.11|) and ()2.16|) . together with the facts that u x u = Q and uj is orthogonal to m x a;. 
We therefore obtain in view of ()3.1|) that 

^{P{r), j) = -rX(n)i/2 det^ loo). (8.3) 
dr 

Lemma 8.6 Let M be a connected component of a regular level set of Tj in S0(3). Let 
det {u, Uj, I ujj) be viewed as a function on M. Then, for any continuous function f on U, 
the integral of f{u) det (m, Uj, I ojj) over M with respect to the area form aj, cf. Proposition 
\8.cl[ is equal to zero. 

Proof We have 

X rju = X{u) (u X ujj) X {u X {I + p) ujj) 

= —X{u) {ujj, u X (I + p) Uj) u = X{u) det {u, Uj, I Uj) u, 

because {u, u x (I + p)lj) = and therefore the Uj-term drops out. Because rj) = 1, it 
follows that aj is equal to the pull-back of 1/X{u) det {u, Uj, I Uj) times the standard area 
form d2 M on f/ by means of the projection n : A \—>- u = A~'^ 63 from the level set M to U. 
Now we have, for any smooth mapping tt from a compact oriented manifold M to a compact 
oriented manifold U, and any volume form Q on U, the formula 

[ 7^*^] = deg(7r) • [ Q, (8.4) 
Jm Ju 

see for instance Guillemin and Pollack p. 188]. In our case 

f{u) det (m, ujj, I Uj) aj = tt* [ {}^\ 



Xiu) 

Moreover, the degree of vr is equal to zero, because n{M) ^ U, see for instance the description 
at the end of Remark 111.21 Therefore the conclusion of the lemma is obtained by applying 
(lOl to = ^ daM. □ 

It follows from Lemma (8.61 that the average of the right and side of ()8.3|) over the level 
surface, with respect to the area form aj, is equal to zero. In view of Corollarv 18.41 we can 
apply ()5.3j) with t replaced by r and p(t) replaced by (p(r), j), in the case that the rotational 
motion is not periodic. In this case the coefficient of the secular term is equal to co(e), which 
is equal to the average of the right hand side of ()8.3p over the level set, with respect to the 
area form aj. This leads to the following conclusion. 
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Corollary 8.7 Assume that the rotational motion on the regular level set is nonperiodic 
and that the series in i5.!^) . with t replaced by r and pit) replaced by {pij), j) , is uniformly 
convergent. Then the function t i— {p{t), j) is quasiperiodic on a two-dimensional torus. 
In particular, {p{t), j) remains bounded in this case. 

Note that the series mentioned in Corollary 18 . 71 converges uniformly when the irrational ratio 
i/i(e)/^'2(e) mentioned after ()5.3|) is sufficiently slowly approximated by rational numbers. 
The set of irrational numbers for which this happens has full Lebesque measure on the real 
axis. 

Question 8.8 Do the complications with secular terms when the rotational motion is 
periodic, and convergence of Fourier series for nonperiodic rotational motions, as discussed 
after (j5.3p . really occur? 

For the critical circles we have the simplification that is constant, cf. ()4.10j) . When 
j is vertical, and when two of the moments of inertia are equal, we can reconstruct the 
full motion by means of Lie group techniques from a periodic motion, which implies that 
the aforementioned complications do not arise. See Proposition 15.11 and Proposition 16. !( 
respectively. 

In the general case I have neither been able to find any analogous special features which 
would eliminate the complications with the primitives of the Fourier series, nor did I find a 
proof that the complications really occur. See also Question I11.121 

8.3 Chap ly gin 

The themes of Section |H1 do not occur in Chaplygin jnj. The only exception may be that the 
existence of two closed one-forms f3 and 7 such that (3 = and i^7 = 1, cf. Proposition 
18. 3| is implicitly contained in the formulas in Chaplygin [9, §3, (29)], cf. Remark 111.31 As 
discussed after R emark 1 1 1 . 31 b elow . the vector field ^ corresponds to an explicitly determined 
constant vector field on the Jacobi variety J{C) of a hyperelliptic curve C. In Remark (jll.4j) 
we determine the constant vector field on J{C) to which ri corresponds. 

9 Some Simplifications 
9.1 A Polynomial System 

The square root in the factor X(-u)^/^ in front of the commuting vector fields ^ and t] in 
Proposition 18.31 becomes double- valued if we extend the vector field holomorphically into 
the complex domain. One can make it single-valued by passing to the space of {z, u, v) e 
C X X C'^ where [u, v) still satisfy ()3.9p and ()3.13|) and the new variable z is coupled to 
u by means of the equation 

z^ = X{u) = p-^ - {u, Ju). (9.1) 
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This means that we pass to the branched (=ramified) double covering of the complexification 
of the level surface, which branches along the complex one- dimensional submanifold (complex 
curve) defined by the equation X{u) = 0. As observed in front of Lemma l7.ll we have 
X{u) > when u is real, which implies that the curve defined by the equation X{u) = has 
no real points. 

We have that C, = Rzlu, V = '^z{i+p)u), where 

zu = z J V + z"-^ Y{u, v) J u, and (9.2) 
z{I + p)uj = zv + z^'^Y{u, v)u, (9.3) 

in which we have used the abbreviation ()3.16|) . We recall that K,^ u = uxu and K^v = vxu. 
The actions on z are given hj 2z ^ z = z"^ = —2z {u x Jv, J u) and rj z = —{uxv,Ju), 
because {u x J u, J u) = and {u x u, J u) = 0. It follows that the vector fields ^ and i] 
are rational. The functions ^ z and r] z are regular, but ^ u and ^ v seem to have poles along 
z = 0. 

Combining the kinetic energy equation in the form (j3.14j) with (|9.1|) . we see that Y{u, v) = 
± z Z(vy^'^, and it follows that zcu = z J v ± Z{v)^l'^ J u and z{I + p)uj = zv ±. Z(t>)^/^ u. 
Note that the sign choice in ± Z(t;)^/^ is coupled to the choice of the sign of z by means of 

dsini). 

However, the vector fields are still double- valued at 2; = 0, where also the manifold of the 
solutions ((m, z), v) of the equations ()3.9j) . ()3.14|1 and ()9.1|1 is singular. These singularities 
can be resolved by introducing one more variable C, which is coupled to v by means of the 
equation 

e = Z{v) = 2T-{v, Jv). (9.4) 

The kinetic energy equation Y{u, f)^ = X{u) Z{v) = z"^ ("^ leads to Y{u, v) = ±z(. With 
the choice of the minus sign, 

Yiu,v)=-zC, (9.5) 
the vector fields ^ = Rzuj and r] = Rz{i+p)u} are given by 

= zuxJv — (uxJu, (9.6) 

= zvxJv — (vxJu, (9.7) 

^ z = —{u X J V, J u) = det('U, J u, J v), (9.8) 

= {v X J u, J v) = det{v, J u, J v), (9.9) 



and 



rju = zu X V, (9.10) 

r]v = —(v X u, (9.11) 

rjz = —{u X V, J u) = det(f, u, Ju), (9.12) 

r/C = X J v) = det(f, u, Jv), (9.13) 



respectively. 
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The equations ()9.6p — ()9.13p define polynomial vector fields ^ and 77 in C^, homogeneous 
of degree three. The vector fields ^ and 77 in commute with each other. Finally both 
vector fields ^ and rj are divergence- free, and have the six functions 



as constants of motion. (We have chosen the minus sign in ()9.5|) in order to get a pkus sign in 
/s.) All these statements are true when J, which appears as the parameter in ()9.6p — ()9.13|) . 
is a symmetric 3 x 3-matrix. 

For the generic values of the fi, the equations (j9.14p define a smooth two-dimensional 
affine algebraic variety M in C^, on which ^ and rj are commuting vector fields. Note that 
Chaplygin's sphere corresponds to the case that 



In contrast with the real case, the complex surface is not compact and therefore we cannot 
conclude that the flows of C, and rj with complex times lead to an identification of M with 
a complex torus. Also, the flows of ^ and 77 on M with complex times will not be complete 
in the sense that these are not defined for all complex times and therefore do not define 
an action of on M. In Proposition 110.81 we will obtain a completion M of M which is 
isomorphic to a complex torus on which the vector fields ^ and t] are constant (and linearly 
independent). A very different construction, based on Chaplygin's integration of the system 
in terms of hyperelliptic integrals, will be described in Subsection 111.41 

The system (|9.6p — (j9.9p and even more so the integrals (j9.14|) resemble the system (2) and 
the integrals (X, X), (3), (5) and (6) in the article of Adler and van Moerbeke 2j. However, 
there are also differences: we have the intersection of six quadrics in an eight-dimensional 
space, whereas in j2] one has the intersection of "only" four quadrics in a six-dimensional 
space. In Subsection 19.31 we will show that the system ()9.6|1 — ()9.9|1 can be mapped to the 
geodesic flow on the Eulidean motion group for a left invariant metric, cf. Subsection 19.31 
The latter system resembles the ones in the article of Adler and van Moerbeke j2j even more 
closely, only with the six- dimensional Lie algebra of SO (4) replaced by the six-dimensional 
Lie algebra of the Euclidean motion group in the three-dimensional space. Like in |2], the 
vector field in this Lie algebra is homogeneous of degree two and has four quadratic constants 
of motion. 

9.2 Reduction to Horizontal Moment 

In this subsection we assume that the moment j of the momentum around the point of 
contact is not vertical. We will investigate how the vector field ^ defined by (j9.6|) — (j9.9j) 
changes if we apply a linear substitution of variables of the form 



/i := {u, u), /a := {u, v), /g := {v, v), 

fi:= {u, Ju) + z'^, f5 := {u, Jv) + zC, f^^ ■= {v, J v) + C 



(9.14) 



/i = 1, /2 = J3, h = lljf , /4 = 1/p, h = 0, /e = 2T. 



(9.15) 



u = au + bv and v = cu + dv. 



(9.16) 
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The equations ()9.16p are equivalent to du — bv = Du, —cu + av = Dv, in which 
D := ad — be. Substituting (|9.6|) and (j9.7p in D = d^u — b^v, D = —c^u + a^v =, 
we obtain with a straightforward calculation that 

D = D {z c — ( a) u X Ju + D{zd — (b)ux Jv, 

and in a similar fashion that 

^v = {z c — ( a)v X Ju + {zd — (b)vx J v. 

These equations are of the form ()9.fi|l . ()9.7|1 . with u, v, z and ( replaced by u, v, z and (, 
respectively, if we take ^ 

z:=zd — (b, —(■.= zc — (a. (9-17) 

Applying ^ to ()9.17|) and substituting ()9.8|) . ()9.9|) and ()9.16|) . we obtain with a straightforward 
calculation that 

^z = D^ det ^M, Ju, Jv^ , —^C = D^ det Jv, Ju^ . 

It follows that the vector field is invariant if we arrange that 

{ad-bcf = = 1. (9.18) 
The constants of motion ()9.14j) are related to the corresponding ones 



/i := {u, u), ^ /2 := {u, v), ^ _ h '■= i^' ^ 
/4 := {u, Ju) +z^, /5 := {u, Jv) + f& := ((v, Jv) + C 



(9.19) 



with tildes over all the variables, by means of the formulas 

/i = A + 2a 6/2 + 62/3, ^ (9.20) 

/2 = ach + {ad + bc)f2 + bdh, (9.21) 

h = ^h + 2cd% + d''h, (9.22) 

/4 = U + 2abf, + b^ fe, (9.23) 

/5 = ac74 + (ad + 6c)/5 + 6ci76, (9.24) 

/e = c^U + 2cd% + d^U (9.25) 

Here we have used in (jHIlSl), (EllD, that (PTTTI) imply that az + b( = Dzand 

cz + d( = D (, whereas ()9.18|) implies that D~'^ = 1. 



Consider the matrices M= , ,F= „ „ ,6*= „ „ .IfFis 

\ c d J \ f2 h J \ h h / 

equal to the matrix F with tildes over the coefficients, then the equations ()9.2()j) . ()9.21|1 . 
()9.22|) are equivalent to the matrix equation F = M F M*. Here M* denotes the transposed 
of M. Similarly, if G is equal to the matrix G with tildes over the coefficients, then the 
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equations ()9.23j) . ()9.24p . ()9.25|) are equivalent to the matrix equation G = MGM*. An 
obvious consequence is that, for any A G C, 

p(A) := det(G -XF) = det(G - A F). (9.26) 

Note that p{X) = a A^ — /5 A + 7, in which 

a = detF = /l/3-/2^ /3 = /i/6 + /3/4- 2/2/5, and 7 = det G = /4 /e - /s'- 

Therefore the invariance of the polynomial p is equivalent to the three equations det F = 
detF, detG = det G and /i /e + hU - 2/2/5 = 7i /e + fsU - 2/2/5. In the case of 
Chaplygin's sphere, we have p{X) = aA^ — /3A + 7 with 

«=l|jir-j3'=Ji'+j2', P = 2T+\\jf/p, and 7 = 2T/p, 

cf. fFTTHll . 

Suppose that all the coefficients are real and that F is positive definite. Then, first 
diagonalizing F by means of an orthogonal transformation, one susbsequently can obtain a 
real 2 x 2-matrix A such that F = A A*. It follows that (det A)'^ = det F > 0, which implies 
that A is invertible. There exists an orthogonal transformation O such that D := G' O = 
G' {0~^)* is diagonal, where G' is equal to the symmetric matrix A^^ G (A^^)*. In other 
words, G = AG' A* = AO D (AO)*. For an arbitrary invertible diagonal matrix B we have 
now arranged that G = MG M*, with M = AO B^^ and G := B D B* diagonal. In order to 
arrange that (detM)^ = 1, it is sufficient to take 1 = {det A)^ {det B)-^ = det F {det B)-^ , 
or {det By = detF. We now have F = M F M* with F = BB*. If we choose B = 
diag (1, (detF)^/^), then we have arrived at the situation that /i = 1, /2 = 0, /a = detF, 

/5 = 0. Note also that in O we still have the freedom to precede it by the matrix which 
switches the two basis vectors, which means that we still can switch and /e. 

Remark 9.1 A diagonal matrix remains unchanged if it is multiplied from the left and 
the right by the matrix diag(l, —1), which has determinant equal to —1. Therefore we can 
arrive at the same diagonal matrices F, G with the help of a matrix M which satisfies the 
stronger condition det M = 1 instead of det M = ±1 as required in ()9.18p . 

In other words, SL(2, C) acts as a symmetry group for the vector fields ^ and rj. The 
functions fi are not invariant under the action of SL(2, C) and actually the SL(2, C)-action 
can be used to change to levels fi such that /i = 1, /2 = 0, /s = det F, /s = 0. On the other 
hand, the coefficients of the polynomial p in ()9.26|) are invariant, which fact can be used in 
order to determine /s, /4 and /e in terms of the fi. 

In Subsection 19.31 we will give a description of the quotient space under the action of 
SL(2, C) 

In the case of Chaplygin's sphere, when we have ()9.15p . we have that det F = js^ = 

ji^ + > because j is not vertical. Therefore F is positive definite and we have a 
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reduction to the situation that /i = 1, /2 = and /s = 0. It follows then from ()9.26p that 
fs = det F = ii' + j2^, whereas /4, /e satisfy the equations /e = det G = 2T/ p and 

2T + \\j\\yp = u + {j,'+j,') u 

These fi are again the levels of a Chaplygin's sphere, but with j, p, T replaced by j, p, T, 
respectively. If js = 0, which means that the new moment j is horizontal, the length of j 
is equal to the length of the horizontal projection of j, and p and T satisfy the equations 
T/p = T/p and 

2f+{n'+j,')/p = 2T+\\jr/p. 

The solution which depends continuously on the parameters and satisfies T = T, p = p when 
js = 0, is given by 

p = pf/T and 2f = T+\\jf/2p + 

Because T ^ T when j is not already horizontal, we have p ^ p- Because we do not change 
the matrix J = {I + p)~^ , this means that we have to change the moment of inertia tensor 
I to the new one I = I + p — p. 

If we allow complex coefficients, then we can arrive at any F, G such that ()9.26p holds, 
provided that the polynomial p is of degree two and has two distinct zeros, cf. Hodge 
and Pedoe Vol. II, p. 278]. These conditions are equivalent to the conditions that 
/i /s — 7^ and the discriminant 

A = (/i /e + hh- 2/2 hf - 4 (/i /3 - h') [h h - h') (9.28) 

of p is not equal to zero. Again we can arrange that /i = 1, f-^ = fi — f2- The new values 
/4 and /e are determined from the condition that and /e/Zs are equal to the zeros of p, 
which are unique up to their ordering. 

Remark 9.2 In the case of Chaplygin's sphere, when we have (j9.15p . and similar equations 
with tildes over all the symbols, the equations (j9.2U|) . (j9.2ip . (|9.22|) are equivalent to the 
statement that there exists a rotation C G S0(3) such that C 63 = aes + hj and C j = 
€63 + dj. These equations are equivalent to ()9.16|) if u = 63, v = j, u = 63 and 
V = A~^j, in which A = C A, A & S0(3). The invariance of the vector field ^ under the 
substitutions ()9.1fij) then means that if A{t) denotes the rotational motion as a function of 
the reparametrized time r, then A{t) = G A{t) satisfies a differential equation of the same 
form as A{t), but with j, p, I and T replaced by j, p, I and T, respectively. 

Remark 9.3 If j is not vertical, then it is impossible to make j vertical. Actually the 
rotational system with nonvertical j can not be transformed in any algebraic fashion to the 
system with vertical moment of Section El 

Indeed, if j is vertical, then the completion of the complexification of the phase space of 
the S0(2)-reduced system is equal to the elliptic curve defined by ()5.9|) . (jS.lOj) . On the other 



(T 



'/2py 



+ '2j3'T/p 



(9.27) 
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hand the complexification of the SO(2)-action is a free action of C^, where denotes 
the multiphcative group of the nonzero complex numbers. (The mapping t ^ is an 
isomorphism from the additive group C/27riZ onto the multiphcative group C^.) In this 
way the completion of the complexification of the phase space of the system for vertical j is 
a C^- bundle over an elliptic curve, which is not isomorphic to the Jacobi variety Jac(C) of 
the hyperelliptic curve C which we obtain when the moment is not vertical. In particular 
the C^-bundle is not compact because is not compact, whereas Jac(C) is compact. 

The passage from nonvertical to vertical j is an example of "a degenerate limit of an 
abelian variety, as an extension of a power of by an abelian variety", mentioned by 
Mumford ^ p. 3.53]. 

For the motion of the point of contact p, we observe that our substitutions imply that, 
with uj as in ()3.1H) . zu = z J v — (Ju = zJv — (Ju = zu. Therefore p.ip yields that 

' ~ z r {{Alu) X 63, j) = z r {uj X u, v) = z r D {lu X u, v) , 



dr 

which is equal to a constant times the same function for the new system with horizontal 
moment, and therefore equal to a constant times the rational function on the double covering 
of Jac(C), described after pi.74|) . Here C denotes the hyperelhptic curve corresponding to 
the new system with the horizontal moment. 

A similar calculations yields for the j x es-component of p that 



d(p, j X 63) 



dr 



DI [a2T + hC/p^ . 



The function i^has a similar description as a rational function on a double covering of Jac(C), 
with 7 replaced by — 1/p. However, the function zC, = —{u, Jv) does not seem to have an 
equally straightforward description in terms of Jac(C). 

9.3 Geodesic Flow on the Euclidean Motion Group 

The form of the equations ()9.6|) — ()9.13p suggests to introduce the vectors 

q := u X V and r:=zv — C^u. (9.29) 

The vector (g, r) G x = represents the exterior product of the vectors (m, z) 
and (f , C) in C^, and therefore the mapping ((m, z), {v, Q) t— > (g, r) from to will be 
denoted by A. 

For Chaplygin's sphere, it follows from ()3.8|) that 

q = uxv = A-^ (eg) x A'^ j = A'^ (eg x j) , 

which means that q is equal to the vector eg x j in body coordinates. Furthermore, it follows 
from ()9.2j) and ()9.5p that Jr = zu, which is equal to the rotational velocity with respect to 
the time variable r which is related to t by dt/ dr = X^/^ = z, cf. Corollarv 18.41 and ()9.H) . 
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We have 

= (^"w) X f + M X (^f ) = {uxJr)xv + ux{vxJr) = {uxv)xJr = qxJr, 

in which the third identity follows from the Jacobi identity (uxr) xv + {r xv) xu + {v xu) xr 
in 00(3). Similarly 

= {u, J u X Jv)v—{v, J u X Jv) V + zv X J r — ( u X J r 
= {u X v) X {Ju X Jv) + r X Jr. 

For any vector w we have that 

{Ju X Jv, J w) = det{Ju, J V, J w) = det J det(-u, v, w) = det J {u x v, w), 

which in view of the symmetry of J implies that J {J u x J v) = det Juxv,oyJuxJv = 
det J J~^{u X v) when J is invertible. It follows that 

Jux Jv = J'^^iuxv), in which = diag ( J2 ^3, ^3 ^1, ^1 ^2) • (9.30) 

Note that J™ = (det J) J~^ when J is invertible, but ()9.30p also holds for noninvertible J. 
This means that the mapping A intertwines the vector field ^ in with the vector field ^ 
in defined by 

iq = q X Jr and = q x J^° q + r x Jr. (9.31) 

Similarly, we have 

rjq = (jju) X V + u X {rjv) = zqxv + ux(q = qx {zv — (u) = q x r, 

and 

rjr = (rjz) v — {ri() u + z (rjv) — ( {rju) = — {Ju, q)v+{Jv,q)u + zvxr — (uxr 
= —{u, J q) V + {v, Jq)u + rxr = — (m x v) x J q = —qx J q. 

Therefore the mapping A intertwines the vector field rj in with the vector field rj in 
defined by 

rjq = q X r and rjr = —q x J q. (9.32) 

If q and r is interpreted as a rotational and translational velocity vector, then the (g, r)- 
-space can be identified with the (complexified) Lie algebra e(3) of the Euclidean motion 
group E(3) in the three dimensional Euclidean space, with the Lie brackets defined by 

[{q, f)^ {q', f')] = {q X q' , q X r' + r X q'). (9.33) 

Therefore the vector field ^ defined by ()9.3H) has the form of a Lax pair dX/ dr = [X, L^{X)], 
where is the linear transformation in e(3) defined by 

L^iq,r) = iJr, J^° q) . (9.34) 
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Similarly the vector field t] defined by ()9.32j) has the form of a Lax pair dX/ dr = [X, L^(X)], 
where L^j is the linear transformation in e(3) defined by 

L,{q,r) = {r,-Jq). (9.35) 

The Lax pair form of ^ and rj implies that the flows of both vector fields leave the conjugacy 
classes in e(3) invariant. Because the functions 

hi{q, r) := (g, q) and /i2(g, r) := (g, r) (9.36) 

are constant on the conjugacy classes, it follows that both functions hi and /i2 are constants 
of motion for ^ and rj. 

Let G be a Lie group and / a function on the cotangent bundle T* G of G which is 
invariant under all left multiplications by elements of G. The quotient space of T* G by 
means of the left action of G on T* G is naturally identified with the dual space g* of the 
Lie algebra q of G, and we denote the restriction of / to g* with the same letter. The 
canonical Possion structure on T* G induces a Poisson structure on g* in such a way that 
the Hamiltonian vector field of / on g* is given by 

-(X, H;(/)) = ([X, d/(0],/), XGg, /Gfl*. (9.37) 

The vector field is tangent to the coadjoint orbits in g*, on which the Poisson structure 
is given by a symplectic structure. This means that on each coadjoint orbit the vector field 
Hj is Hamiltonian with respect to this symplectic structure. This construction has been 
introduced already by Lie in [25, Kap. 19] under the name "Die dualistische der adjungierte 
Gruppe". It has been rediscovered independently by Kostant, Kirillov, and Souriau. The 
coadjoint orbits with their symplectic structure are the Marsden-Weinstein reduced phase 
spaces of T* G for the left action of G on T* G. See also Abraham and Marsden ^1 Sections 
4.3, 4.4]. 

When g = e(3), cf. (j9.33l then the right hand side in (j9.37p takes the form 



Therefore, if we identify the point (a, h) G g* with the point (g, r) G g with q = h and r = a, 
then we recoginze from ()9.31|) that ^ is equal to the Hamiltonian vector field of the function 
/ = /i3/2, in which 

hiq,r):={q,r°q) + {r,Jr). (9.38) 

With the same identifications the vector field r] is equal to the Hamiltonian vector field 
defined by the function 

h4{q,r):= -{q, Jq) + {r,r). (9.39) 
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Remark 9.4 The Lie algebra e(3) has not been identified with its dual e(3)* by means of 
the Killing form — 4(g, g), which is degenerate, but by means of the nondegenerate invariant 
quadratic form /i2 = ((?, r). 

The function defines a left invariant metric on T* E(3), which can be used to identify 
T* E(3) with the tangent bundle of E(3). Under this identification the fiow of the Hamilto- 
nian system of corresponds to the geodesic fiow on TE(3) defined by the dual metric 
on TE(3). 

We therefore obtain the following, somewhat roundabout correspondence between this 
geodesic fiow and our vector field ^. First pass from the geodesic fiow of the left invariant 
metric on the tangent bundle to the Hamiltonian system of the function /13/2 on the cotangent 
bundle, using the metric in order to identify the tangent bundle with the cotangent bundle. 
Then pass to the reduced system on e(3)* by means of the left action of E(3) on T* E(3). In 
the next step, pass to the vector field ^ on e(3) using the identification of e(3) with e(3)* by 
means of the quadratic form (g, r). Finally the mapping A intertwines the vector field ^ in 
the ((u, z), (f , C))-space with the vector field ^ in the (g, r)-space e(3). 

The Poisson brackets {/, g} := H/(7 define a Lie algebra structure on the space of 
functions. In particular it is antisymmetric, which implies that Hj/ = and H/^f = 
if and only if / = 0. It follows from Hj / = that and is a constant of motion for 
^ and ?7, respectively. Moreover, 

77/13/2 = (r^g, J™g) + (r/r, J r) = {q x r, J""" q) — (g x Jg, Jr) = 0, 

because 

(g X Jg, Jr) = (g, Jg x Jr) = (g, J""" (g x r)) = (J''°g, g x r), 

cf. fl9.3m) . This implies that and are constants of motion for both vector fields ^ and 
rj. The Jacobi identity of the Poisson structure implies that [H/, H^] = H{j g}. Because we 
just proved that {h^, h^} = 0, it follows also that the vector fields ^ and rj commute. 

In particular the system in e(3) ~ defined by the vector field ^ is completely integrable, 
as a Hamiltonian system on the four- dimensional coadjoint orbits determined by fixing the 
values of the functions hi and defined in (I9.36|) . The Hamiltonian function is the function 
/i3/2 with /13 defined in fl9.38p and the function /14 defined in fj9.39p is the additional constant 
of motion which yields the complete integrability of the Hamiltonian system. 

Remark 9.5 The system in the (g, r)-space resembles the one in the paper of Adler and 
van Moerbeke very closely: it is defined on a six-dimensional (dual of a) Lie algebra 
(e(3) instead of the Lie algebra so(4) of [2]) and has four quadratic constants of motion. 
Furthermore it is Hamiltonian on coadjoint orbits, and we have two polynomial vector fields 
which are homogeneous of degree two. In this respect the vector fields are simpler than the 
vector fields ^ and rj defined in fj9.6|) — ()9.13p . 

Also the behaviour at infinity the level surface of the hi (with /i2 = 0) in the complex 
projective space is very similar to the behaviour of the level surface in so(4) as described in 
Mumford's appendix to ^ . See Subsection I1U.51 
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On the other hand it turns out that this behaviour is more singular than that of the 
projective closure of the level surface M of the functions (j9.14j) . see Subsection I1U.51 For 
this reason we start Section^] with a discussion of the system in C^, because this seemingly 
more complicated system has a simpler behaviour at infinity. 

We now turn to a closer examination of the mapping A from to defined by fl9.29|) . 
To begin with, if (g, r) belongs to the image of A, then (g, r) = 0, which means that A is 
a mapping from to the hypersurface /t2 = in C^, cf. ()9.36|) . Conversely, any g G 
can be written as m x t> for some u, v E C^. If g 7^ 0, then u and v are linearly independent 
and span the orthogonal complement of g, which implies that for every r G C'^ such that 
(g, r) = there exist z, ( E C such that r = z v — ( u. If g = 0, then we have ()9.36|) for 
z = 1, V = r, ( = 0, u = 0. It follows that A is surjective from onto the hypersurface 
/i2 = in 

Remark 9.6 If g 7^ 0, then (g, r) = if and only if there exists a vector xq such that 
g X a;o + r = 0, and every vector x such that gxa; + r = Oisof the form x = Xo + cq for some 
scalar c. Therefore the condition that (g, r) means that the infinitesimal motion (g, r) G e(3) 
is either equal to an infinitesimal translation (g = 0), or to an infinitesimal rotation about 
some axis in the three-dimensional case: "no spiralling". 

If we interpret (g, r) as an element of /\^ C^, then the condition (g, r) = means that the 
rank of (g, r) is smaller than four. If (g, r) 7^ (0, 0), then its null space is two-dimensional 
and is spanned by {u, z) and (f, Q if (g, r) = (m, z) /\{v, Q = {u x v, zv — (u). Note 
that (m, z) and (f , () are linearly independent in this case. If (g, r) = («', z') A {v', ('), 
then [u', z') and {v', (') are also contained in the null space of (g, r), and therefore there are 
unique a, b, c, d E C such that 

[u', z') = a{u, z) +b{v, C) and [v' , (') = c{u, z) + d {v, (), (9.40) 

and for such vectors {u', z') and (f', C,') we have that (g, r) = (u', z') A (f', C,') if and only 
ii ad — he = 1. In other words, if (g, r) = and (g, r) 7^ (0, 0), then the fiber of (g, r) for 
the mapping A is equal to the orbit in the ((«, z), (f , C))-space of the action of SL(2, C) 
defined by (ICTIl . 

Let U denote the set of the ((-u, 2;), (t>, Q) G such that the vectors {u, z) and {z, Q 
in are linearly independent, and let V be the set of (g, r) G such that (g, r) = 
and (g, r) 7^ (0, 0). Then U is an open subset of and V is an open subset, equal to the 
smooth part, of the 5-dimensional hypersurface (g, r) = in C^. The action of SL(2, C) on 
U is free and the mapping A identifies V with the orbit space of the SL(2, C)-action on U . 

In Subsection 19.21 we had observed that the action of SL(2, C) leaves the vector fields 
^ and 7] invariant, cf. Remark 19.11 Therefore the projection /\ : U ^ V intertwines the 
vector fields ^ and rj on U with uniquely defined vector fields on V, which we denoted with 
the same letters. The fact that the vector fields ^ and t] in U commute implies that their 
push-forwards under A, the vector fields ^ and rj in V, commute as well. In the beginning of 
this subsection we showed that the vector fields ^ and rj in defined by (j9.31|) and (j9.32j) 
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extend the vector fields ^ and t] in V. The fact that the vector fields ^ and r] in commute 
is stronger than the fact that their restrictions to V commute. 

As observed in Remark l9. 11 the functions fi are not invariant under the SL(2, C)-action. 
However, the coefficients of the polynomial p in ()9.26p are invariant under the SL(2, C)- 
action, which means that these coefficents can be written as functions of (g, r) G V. Actually, 
we have: 

h,oA = hh-h\ (9.41) 
/120A = 0, (9.42) 
/i3oA = hh-h\ (9.43) 
(traceJ)/iioA + /i4oA = hh + hf^-V^h. (9-44) 

from which the compositions /ijO A of the functions hi with the mapping A can be determined 
in terms of the functions fi defined in (j9.14p . Here the functions hi are defined by (|9.36p . 
(FHIHll and (PM . 

Proof The equation 1)9.411) follows from 

det F = /i /s - = (n, u) ■ {v, v) - {u, v)^ = {ux v, ux v) = {q, q) = hi. 

For ()9.43|) we write 

detG = Uh-h''={{u,Ju) + z^) {{v,Jv)+e)-{{u,Jv) + zCf 

= {u, J u) ■ {v, J v) — {u, Jt>)^ + z'^ {v, J v) + (m, J u) — IzC, (m, J v) 

= {u, {v, J v) J u — {v, J u) J v) + {z V — ( u, J {zv — (u)) 

= {u, V X {Ju X J v)) + (r, J r) = (g, J""" q) + (r, Jr) = hs, 

where we used fl9.30p in the fifth identity. 

Finally, we have that fi /e + fs fi — 2/2 /s is equal to 

{u, u) [{v, Jv) + C^) + (1;, v) ((m, Ju) + z^) - 2(m, v) ((n, J v) + zC) , 

for the computation of which we write 

(m, u) ■ {v, J v) + (t>, v) ■ {u, J u) — 2{u, v) ■ {u, J v) = {u, a), 

in which 

a = {v, J v) u + {v, v) J u — 2{u, J v) V 

= {v, J v) u — {u, J v) V + {v, v) Ju — (f , J u) V 
= Jv X {u X v) + V X {Ju X v). 

It follows that {u, a) = {u X v, b), in which b = uxJv + Juxv. Now we have, for any 
vector c, 

(6, c) = det(u, J V, c) + det(JM, v, c) = (trace J) det(u, v, c) — det(M, v J c) 
= (trace J) (m x v, c) — {u x v, J c) = ((trace J — J){u x t>), c). 
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which imphes that b = (trace J — J){u x v). Collecting all the results, we arrive at 



/i /e + fsh- 2/2 /s = {q, (trace J - J){q)) + {r, r) = (trace J) hi + h. 



'4, 



from which ()9.44|) follows. 



□ 



In the case of Chaplygin's sphere, we can insert the values ()9.15|1 of the functions fi, 
which leads to the values 



for the constants of motion hi of the system in the six-dimensional (g, r)-space. 

The action of the matrix M G SL(2, C) with the coefficients a, b, c, d leaves a given level 
surface of the functions fi invariant if and only if, in the notation of Subsection 19.21 we have 
that M F M* = F and M G M* = G. Assume that detF 7^ and that F and G are not 
proportional. Then we obtained in Subsection 19.21 that there exists an A G SL(2, C) such 
that the matrices F = AF A* and G = AG A* are diagonal, and it follws that F and G 
axe not^ proportional. With the notation M = AMA~^, we now have M F M* = F and 
M G M* = G. A straighforward calculation, in which we u^ that the diagonal matrices F 
and G are not proportional, leads to the conclusion that M = ±1, which in turn implies 
that M = ±1. If M = —1 then it acts on the {{u, z), {v, C))-space as the antipodal map 
X ^— s> —X. 

It follows that the restriction of the mapping A to the level surface M of the function 
fi, where we assume that detF 7^ and F and G are not proportional, defines a twofold 
unbranched covering from M onto the level surface of the functions hi, for the levels given by 
fl9.41|) — ()9.44p . The fibers of A\m are pairs of antipodal points, and therefore the mapping 
A leads to an identification of the level surface A(M) on /i2 = of the functions hi, /13, /14 
with the quotient M/±l of the surface M by means of the antipodal mapping. 

Because the commutation of vector fields is a local property, the fact that the vector 
fields ^ and t] on A(M) commute implies, together with the fact that A : Af — > A(M) is a 
covering, that the vector fields ^ and t] on M commute. This leads to a proof of Corollary 
18.41 which is based on the facts that ^ is a Hamiltonian vector field on a coadjoint orbit and 
has the functions fi as constants of motion. 

9.4 Symmetric Matrices 

The equations of motion ()3.2|) . in which u = ujj{A) is given in terms of A by ()2.20|) . ()2.19|) . 
(|2.12p . can be entirely expressed in terms of the (positive definite) symmetric matrix 



hi = ]i^+]2\ h2 = 0, h3 = 2T/p, (trace J) hi + h, = 2T + \\jf/p 



(9.45) 



B:=A{I + p) 



'A 



-1 



(9.46) 



We have 



dB 



[e(i?)op, B] :=e(5)opoB-Boe(B) 



(9.47) 



dt 



op- 
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Here ^(S)op denotes the antisymmetric linear mapping v ^ ^(-B) x v and 



i{B) := Au.iA) = Bj + " Be,. (9.48) 

1- p{Be3, 63) 

The velocity of the point of contact p also is a function of B: 

^ = re(5)xe3. (9.49) 
at 

The equation ()9.47p is a Lax system, cf. [21], and therefore the eigenvalues of B are 
constants of motion. More explicitly, it follows directly from ()9.46|) that we have the the 
following three constants of motion in the six-dimensional vector space of the symmetric 
3 X 3-matrices: 

traced* = trace (/ + p)-\ i = 1, 2, 3. (9.50) 

These equations are homogeneous of degree one, two and three, respectively. The kinetic 
energy equation ()3.14|1 is a polynomial equation of degree two in B: 

{Bj, es)^ - (p-i - {Bes, e,)) {2T - {Bj, j)) = 0. (9.51) 

When the moment of inertia tensor J, assumed to be diagonal, has three different eigen- 
values, then the mapping A B = A {I + p)^^ A^^ is a fourfold covering from SO (3) onto 
the manifold determined by ()9.50|) . where A and A are mapped to the same symmetric ma- 
trix if and only A! = AoRm which R is one of the four diagonal rotations. Apparently the 
reduction of this symmetry group leads to quite a reduction of the degrees of the constants of 
motion. After passing to a covering on which the vector field X is single valued and regular. 
Proposition II 1.61 would lead to a mapping to a fourfold covering of the Jacobi variety of a 
hyperelliptic curve of genus two, on which X corresponds to a constant vector field. 

The system fl9.47p has a strong resemblance to the equations of van Moerbeke ^Hl formula 
(17)] in the lowest dimensional case = 3. In his case the kinetic energy equation 1)9.5111 is 
replaced by the condition that the "modulus" , the product of the upper triangular elements 
of B is kept constant. His equations of motion also lead to a constant vector field on the 
Jacobi variety of a hyperelliptic curve of genus two. 

9.5 Chaplygin 

Subsection 19.21 reflects our understanding of Chaplygin's [01 §4]. His A, A', /i, p! in (35) 
correspond to our a, 6, c, d in ()9.16|) . The equations in his (37) correspond to our ()9.20|) . 
()9.21|) . ()9.22p in which the /j and /j are given by ()9.15|) and the same formulas with tildes over 
all the variables. The equation ID X X' + p p' = after Chaplygin's formula (40) corresponds 
to our ([^O^ with h = h = 0. 

The sentence "The sphere rolls in a direction perpendicular ..." in front of [HI (47)] has 
not been formulated very accurately. The velocity of the point of contact is neither exactly 
orthogonal to the moment, nor is it a periodic function of (the reparametrized) time. 
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The formula (47)] corresponds to our description of d{p, j) / dr. The variable ^ in P 
§4] corresponds to our {p, j x 63), but our equations for it differ from the equations which 
Chaplygin obtained for it at the end of P| §4] . In the very last formula in P, §4] , we believe 
that the right hand side has to be replaced by its primitive with respect to the time t, an 
expression which is not much more transparent than our constant times the primitive with 
respect to r ofz(. 

10 Complexification and Completion 
10.1 A Smooth Complex Surface 

Let M = M(c) denote the set of the solutions ((u, z), {v, ()) G of the equations fi = Ci, 
in which fi are the quadratic forms defined in ()9.14|) . In this subsection we will assume that 
J is a diagonal matrix with three different eigenvalues Ji, J2, J3 on the diagonal. We will 
furthermore assume that the constants q satisfy 

ci C3 7^ 0, C2 = 0, C5 = 0, ci C6 7^ C3 C4, (c4 - Ji ci) (c6 - Ji C3) 7^ for z = 1, 2 and 3. 

(10.1) 

As we have seen in Subsection 19.21 we can arrive at the first three conditions ci C3 7^ 0, 
C2 = 0, C5 = and ci ce 7^ C3 C4, if and only if, in the original system, det -F = /i /3 — 7^ 
and the discriminant A of the polynomial p : A 1— > det(G' — A F), given by ()9.28p . is not equal 
to zero. Because the polynomial p is invariant under the transformations in Subsection 19.21 
the last condition means that none of the Jj is a zero of p. Summarizing, the conditions 
mean for the original system that the polynomial p is of second order and has two distinct 
zeros, none of these equal to one of the Jj's. In Subsection 19.21 we actually arranged also 
that Ci = 1, as one always has for Chaplygin's sphere. 

In the case of Chaplygin's sphere, where we have (j9.15|) and Ji = 1/ (li + p), these 
assumptions mean the following. 

i) The Ji are different: there are three different principal moments of inertia /j. 

ii) /i /3 — = IIj P — js^ = ji^ + 32^ 7^ 0: the moment vector j is not vertical. 

iii) A = (2T — II + Sjs^ T / p ^ Q: this is automatically true when J3 7^ 0, or j is not 
horizontal, because this implies that j 7^ 0, which in turn implies that T > 0. When 
j is horizontal, then the critical energy levels are equal to Tcnt,i = ||i||^/2 (/« + p) < 
||j|p/2p, cf. fl4.9p . Because T is less than or equal to the largest critical energy level, 
we have that 2T < ||j|p/2p when js = 0, and therefore we have always that A > 0. 

iv) For Chaplygin's sphere, we have that 

P(A) = (lljf - Js') X' - (2T + IIjIIVp) a + 2T/p. 

The equation ()4.9p for t = Ii, Ji = 1/ {Ii + p), turns out to be equivalent to p (Ji) = 
when T = Tcrit,j. Therefore the last condition in ()10.H) follows from the condition that 
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T is not equal to one of the critical energy levels, or that the real part of the level set 
is a smooth two-dimensional manifold. 

After the reduction to the situation that the moment j is a nonzero horizontal vector, we 
also have that C2 = 0, whereas C5 = always holds for Chaplygin's sphere. We conclude that 
for Chaplygin's sphere the conditions are satisfied if the moment vector j is not vertical and 
T is in between the critical energy energy levels. 

Proposition 10.1 The derivatives dfi of the functions fi defined in ()9.14j) are linearly in- 
dependent at each point of M, and therefore M is a smooth complex two-dimensional affine 
algebraic variety. 

Proof We have to prove that if {{u, z), {v, ()) G M and 1 < i < 6 are constants such 
that Yl^=i '^i ^fi ~ 0' then all the are equal to zero. The equations for the amount to 



2ai M + 0:2 f + 2a4 Jm + as Jt; = 0, (10.2) 

2^4 -2 + as C = 0, (10.3) 

a2 u -\- 2a3 V -\- J u -\- 2ae J V = and (10.4) 

a5 z + 2aeC = 0, (10.5) 



corresponding to the derivatives with respect to u, z, v and (, respectively. 
If we take the inner product of (jl0.2j) with v, we obtain that 

= a2 C3 - 2a4 2: C + "5 (c6 - C^) = "2 C3 + as ce, 

where we have used {v, u) = 0, {v, v) = C3, {v, J v) = —z(, {v, J v) = Cq — C^, and (|10.3p . 
If we take the inner product of p0.4|) with u, we obtain that 

= a2 ci + as (c4 — z^) — 2a6 2 C = «2 ci + as C4, 

where we have used {u, u) = 1, {u, v) = 0, {u, J v) = —z(, {u, J u) = C4 — 2;^, and (jl0.5p . 
These two equations for a2 and as lead in combination with C3 C4 7^ ci ce to the conclusion 
that a2 = as = 0. 

If we substitute this in fjl0.2p then it follows that, unless ai = a4 = 0, the vectors Ju 
and u are linearly dependent, which in turn implies that J u = JiU ioi some i = 1, 2, 3. It 
follows that {u, J u) = Ji {u, u) = JjCi, and therefore 

^2 = C4 - JiCi ^ 0. (10.6) 

From (jl0.3p with as = we obtain that a4 z = 0, which in view of (jl0.6|) implies that a4 = 0. 
Now ()10.2p is equivalent to ai m = 0, which implies that ai = because {u, m) = ci 7^ 
implies that u ^ 0. 

Unless a3 = ae = 0, it follows from ()10.4|) . in which a2 = as = 0, that Jv and v 
are linearly dependent, which implies that Jv = JiV for some i = 1, 2, 3. It follows that 
{v, J v) = Ji {v, v) = Ji C3, and therefore 

C'=C6- JiC3^0. (10.7) 
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From ()10.5|) with ^5 = we obtain that a^^ = 0, which in view of ()10.7|) imphes that ag = 0. 
Now (|lU.4p is equivalent to a^v = 0, which imphes that 03 = because {v, f ) = C3 7^ 
imphes that v ^ 0. □ 



Proposition 10.2 The polynomial vector fields ^ and rj are linearly independent at every 
point of M. 

Proof If M X Ju = 0, which means that u = UiCi for z = 1, 2 or 3. This imphes p0.6|) 
and therefore z ^ 0. Because {u, v) = 0, it foUows that the vectors C,u = zu x J v and 
rju = zu X V can only be linearly dependent if v and Jv are linearly dependent, which 
implies that v = VhCh for some = 1, 2 or 3. This implies (jlU.Tp with i replaced by h and 
therefore C 7^ 0. On the other hand we have that Ui^ = {u, m) = ci 7^ 0, = {u, v) = UiVi 
hence f j = and therefore —zC, = {u, J v) = Ui JiVi = 0, which leads to a contradiction. In 
a similar way we obtain that ^ and r] are linearly independent when v x Jv = 0. 

In the sequel of the proof we therefore may assume that u and J u are linearly independent 
and that v and Jv are linearly independent. Assume that = 0, which in turn implies 
that (u, J v) = —zC = Q, and therefore {u x J v) x u = Ci J v, whereas {u x v) x v = CiV 
because {u, v) = 0. It follows that u x Jv and u x v are linearly independent. If z 7^ 
then C = and we obtain that = zu x Jv and rju = zu x v are linearly independent. 
If C 7^ then 2 = and we obtain that E,v = —(v x Ju and rjv = —(v x u are linearly 
independent. 

li z = ( = 0, then it follows from = —z( = {u, Jv) = {Ju, v) and {u,v) = that 
there are nonzero a, b E C such that u = av x J v and v = bu x J u. Inserting this in ()9.8|) . 
(El), we obtain that 

b^Z = {v, Jv) = Cq - C"^ = C6, 

— a^C = {u, J u) = C4 — z ( = Ci, 

brjz = {v, v) = C3, and 
—ari( = {u, u) = Ci. 

It follows that 

-a b [(^ z) {ri C) - C) iv z)] = Ci Ce - C3 C4 7^ 0, 

cf. fllO.ip . which implies that ^ and rj are linearly independent. 

Finally, suppose that uxJuy^O,vxJ^O,z^O, (^0, and a ^ + /5 r/ = for some 
a, P E C. It follows that there are 7, 5 G C such that 

a{zJv — (Ju) + Pzv = ■JU and (10.8) 
a {zJv - CJu) - PCu = 5v. (10.9) 

If we take the inner product of (jl0.8j) with v, then we obtain that 

= a [z [ce - C^] + C ^ C) + /5 -2 C3 = 2 (a ce + /? C3) , 
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where we have used that (f, u) = 0, {v, J v) = cq — C^, {v, J u) = —z(, and {v, v) = C3. 
Similarly the inner product of (|1(J.9|) with u yields 

= a {-z^ C - C [ci- z^]) - pCci = -C (a C4 + /? ci) . 

Because cq ci 7^ C4 C3, 2 7^ and C 7^ 0, these two equations for a and /5 imply that a = P = 0. 
□ 

Let L denote the manifold of the solutions {u, v) e x of (jT^ and The 
projection {{u, z), [v, ()) ^ {u, v) exhibits M as a two-fold covering of L and intertwines 
the vector fields 1] and rj defined by ()9.6|) — ()9.13p with the previously defined unoriented 
vector fields ^ = X(uY^'^ and t] = X(m)^/^ R.(/+p)aj on Mq, where the word "unoriented" 
refers to the fact that the latter vector fields are only determined up to their signs. 

10.2 At Infinity 

The (homogeneous) polynomial nature of the vector fields ^ and rj, given by ()9.6p — ()9.9p 
and (j9.1(J|) — (j9.13|) . respectively, and of the constants of motion (j9.14p makes it natural to 
investigate the system in the eight-dimensional complex projective space CP*^. The complex 
projective space is obtained by adding one more variable, which we denote by e, and then 
taking the quotient of \ {0} with respect to the actions x ^ cx of the multiplicative 
group of the nonzero complex numbers c. The standard coordinate charts correspond to 
the quotients of the sets of x for which one of the coordinates, say Xi, is nonzero, and then 
the coordinates for this chart are obtained by putting Xi = 1 and using the the Xj with j i 
as the coordinates. The changes of coordinates are obtained by using the identification of x 
with ex. Although it would be clearer to do so, we will not introduce different notations for 
the coordinates in the various charts, in order to avoid heavy notations. If we put e = 1, then 
we obtain the affine space as a subset, equal to one of the standard coordnate charts, 
of CP^. The complement CP^ := CP*^ \ corresponds to taking e = in the other 
coordinate charts, in which one of the coordinates of ((m, z), {v, ()) is taken equal to 1. In 
this way CP^ \ is identified with CP^. We will refer to C CP^ and CP^ - CP^ as 
the affine (or finite) part of CP^ and the projective space at infinity, respectively. 

For any choice of the constants Cj, we denote by M(c) the set of the solutions in C® of 
the equations /i = Cj, the level set of the constants of motion. As a subset of CP^, the set 
M(c) is obtained by homogenizing the equations. In view of the homogeneity of the fi of 
degree two, this corresponds to replacing the equations fi = Ci by 

g,{{u, z), {v, 0, e) := M{n, z), {v, ()) - Q = 0. (10.10) 

Let N{c) denote the set of solutions of (|l().10jl in CP*. Note that M(c) = N{c) n C* = 
iV(c) \ CP^ is equal to the affine part of N{c). 

Our goal in this section is to study the closure M(c) of M(c) in CP*, especially in the 
case that M(c) = M with c and M as in Subsection 110. II Here the closure is taken with 
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respect to the ordinary topology, but it is known that M(c) is equal to a projective algebraic 
variety, and therefore also closed in the Zariski topology. See Lojasiewicz j2ni P- 383]. 

The solutions of the equations (jlU.lU|) at CP^ are obtained by putting e = in (jlU.lUp . 
in which case we obtain the equations 

{u,u) = 0, {u,v) = 0, {v,v) = 0, nnn^ 

{u, Ju) + z^ = 0, {u, Jv)+zC = 0, {v, Jv) + C = 0. ^ ' 

The solutions of ()10.11|) in form the conic affine algebraic variety M(0), and the cor- 
responding projective variety in CP^, which we denote by M(0)oo, is equal to the set of 
the solutions in CP^ of the equations (ri().l()|) . Note that M(0)oo = N{c) n CP^ does not 
depend on the choice of the Cj. 

Lemma 10.3 The set M(0) is a three-dimensional conic subvariety ofC^, consisting of the 
{{u, z), {v, 0) such that {u, u) = 0, {u, Ju) + z'^ = and the vectors {u, z) and (f, Q in 
are linearly dependent. Ifu = and z = 0, then we have to add the conditions that {v, v) = 
and {v, J v) +(^ = 0. The corresponding projective variety M(0)oo = ^(c) flCP^ at infinity 
is a smooth two-dimensional subvariety o/CP^. 



Proof If M 7^ 0, then {u, v) = implies that v = u x w for some w G C^. Using {u, u) = 0, 
we obtain that 

= {v, v) = {u X w, u X w = {u, u) ■ {w, w) — {u, = — (n, w)"^, 

or {u, w) = 0, which in turn implies that w = u x a for some a G C^. But then 

V = u X (u X a) = {u, a) u — {u, u) a = {u, a) u, 

which shows that v = Xu for some A G C. 

Now assume conversely that u, z are solutions of (m, u) = 0, {u, J u) + z"^ = and 
that V = Xu for some A G C. Then we have automatically that {u, v) = X {u, u) = and 
{v, v) = X^ {u, u) = 0, whereas the equations 

= {u, Jv) + zC = X{u, Ju) + z( = -Xz^ + z( = z{( - Xz), 

= {v,jv)+e = x'{u,ju)+e=-x'z'+e = {c+xz){c~xz) 

are equivalent to( = XzoTz = (-\-Xz = 0. In the second case z = ( = 0, and therefore the 
conclusion is that the equations (jlU.llj) hold if and only if ( = Xz. 

In a similar way we obtain that if v, ( are solutions of {v, v) = 0, {v, Jf ) + = and 
M = /if, then the equations p0.11|) hold if and only ii z = For yU 7^ this corresponds 
to the solutions in the previous paragraph with A = whereas for /i = we obtain the 
missing solutions with u = 0, which implies that 2; = in view of = {u, u) + z"^ = z^. □ 

Suppose that M(c) = M with c and M as in Subsection llO.il Then M is a smooth two- 
dimensional affine algebraic variety, and its closure M in CP® with respect to the ordinary 
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topology is a projective algebraic variety, cf. Lojasiewicz [211 P- 383]. It follows that the 
intersection Moo '■= M n CP^ of M with the projective space at infinity, the set of the 
limit points of M at infinity, is an algebraic variety in CP^. It is known that in general 

dim Moo = dimM — 1, cf. Lojasiewicz [26, p. 388], and therefore Moo is aan algebraic curve 
in the projective space at infinity. (Actually, the explicit computations below lead to an 
independent verification of this, see Proposition 110.71 ) 

It follows from Lemma ITTOl that N{c) = M U M(0)oo, which implies that M C N{c), 
or Moo = M n M(0)oo- N{c) is not irreducible, because it has the two-dimensional varieties 
M and M(0)oo, which intersect along the curve Moo, as proper components. As we will 
see below, the curve Moo = M{c)oo depends on the choice of the constants Cj, and actually 
the surface M(0)oo is equal to the union of the curves M(c)oo for the various c's such that 

dim M(c) = 2. 

The fact that M(0)oo is higher-dimensional than M(c)oo is surprising, because for generic 
polynomials Qi the codimension of M(0)oo in the projective space at infinity is equal to the 
number of the equations. Because the codimension of M(c)oo C M(0)oo cannot be larger, 
it follows that, for generic Qi, M(c)oo is equal to the union of some of the components of 
M(0)oo- Compared to this, our set of polynomials gi is quite degenerate. We still have to 
determine, in the case that dim M(c) = 2 < dimM(O), which curve in the projective surface 
M(0)oo is equal to the limit curve M(c)oo of M at infinity. 

Remark 10.4 At the subset M(0), both vector fields ^ and t] are equal to zero. Actually, 
the set where both ^ and t] are equal to zero is much larger. One component consists of 
the ((m, z), (f, C)) for which the vectors (m, z) and (f, () in are linearly dependent, this 
component is five-dimensional. At this component the values Cj of the functions have the 
property that the matrices 



both are singular and one is a multiple of the other. 

Further components consist, for each z = 1, 2 or 3, of the ((m, z), {v, ()) such that u and 
V are multiples of Cj and z and ( are arbitrary. These components are four-dimensional. At 
this component the matrices F and G — JiF are singular. 

It follows that the condition that detF ^ 0, which allowed us in Subection 19.21 to make 
a reduction to the case that C2 = and ci C3 7^ 0, implies that we avoid the variety where 
both vector fields vanish. 

Note that if u, z, t>, ^ are real, then the equations (u, m) = and (f , f ) = imply that 
u = and f = 0. Subsequently the equations (m, Jm) + = and (f, Jf ) + = imply 
that z = and C = 0- Therefore M(0)oo has no real points, which implies that Moo has no 
real points either. This corresponds to the fact that the real points of M form a compact 
subset of R^. 

Suppose that we are at a point p = {{u, z), (t>, (), e) of M(0)oo, where e = 0, {u, u) = 0, 
{u, J u) + = 0, (m, z) 7^ (0, 0), and, for some A G C, f = Am and C, = \z. Recall that in 




and 




(10.12) 
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the standard charts we have to put one of the coordinates of m, v, ( identically equal to 
1. The case u = 0, z = is covered by interchanging the role of (m, z) and {v, (). 

At such a point p the equation Yl^=i ^9iip) = for the constants 1 < i < 6, 
amounts to the equations 

= (2ai + A 02) u + (2a4 + A a^) J u, = (2«4 + A a^) z, 
= (a2 + 2A tts) u + (0:5 + 2A a^) J u, = (0:5 + 2A oq) z, 

cf. fjl0.3|) — ()lU.5p . If M 7^ 0, then Jm is linearly independent of u, because otherwise u would 
be a nonzero multiple of one of the basis vectors Cj, in contradiction with {u, u) = 0. In that 
case the equations are equivalent to the four equations 

2ai + A aa = 0, 2^4 + A = 0, 0^2 + 2A = 0, ag + 2A ag = 0. (10.13) 

On the other hand, if u = then z ^ and the same conclusion holds. It follows that at 
all points of Moo the rank of the matrix of the df^j's is equal to four, instead of the expected 
five. 

The equations (jl(J.13|) are equivalent to 0^2 = — 2Aa;3, ai = A^ 0^3 and = — Aag, 
0:4 = A^ ag, in which and ag are free. In other words, at the aforementioned points of 
^(0)00 we have that the derivatives at p of A^ gi — 2\g2 + gs and A^ g^ — 2Xg^ + g^ are equal 
to zero. This implies that dg^i^p) and dgei^p) are equal to linear combinations of the dgj{p) 
with j = 1, 2, 4, 5. 

Let B denote the common zeroset in CP® of the gj with j = 1, 2, 4, 5. Because the 
dgj{p) are linearly independent, we have that near p the set i? is a smooth four- dimensional 
complex projective subvariety of CP®. The tangent space of i? at p is equal to the common 
null space of the dgj{p), j = 1, 2, 3, 4, which in turn is equal to the common null space of 
the dgi{p), 1 < i < 6. Note that in the following lemma one of the coordinates of p is kept 
equal to zero, corresponding to the projective coordinate chart in which we are working 

Lemma 10.5 The tangent space of B at p is equal to the common null space of the dgi{p), 
1 < i < 6. It consists of the vectors p = {'z, u, v, C, e); such that {u, u) = 0, {u, v) = 0, 
{J u, u) + zz' = and {J u, v) + z ( = 0. 

Proof The first equation is equivalent to dgi{p)p = 0. Assuming v = Xu and the first 
equation, the second equation is equivalent to dg2{p)p = 0. The equation dg3{p)p = 
follows from the combination of v = Xu, the first and the second equation. 

The third equation is equivalent to dg4{p)p = 0. Assuming v = Xu and the third 
equation, the fourth equation is equivalent to dg5{p)p = 0. The equation dgQ{p)p = 
follows from the combination of v = Xu, the third and the fourth equation. □ 

Assume that -2 7^ 0, which means that we can work in the projective coordinate system 
for which z = 1. In this case ( = X z = X. Therefore, if we define the polynomials g and h 
by 

h := C'^7i-2C^72 + ^73 = (C^-^,C'«-^)-(ciC'-2c2C + C3) and (10.14) 
k := C'^74-2C^75 + ^?6 = (Cw-^^, ^(C^-^^))- (C4C'-2C5C + C6) e^, (10.15) 
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respectively, then dh{p) = and dk{p) = for all p G M(0)oo- Moreover, the set N{c), the 
common zeroset in CP® of all the Qi, 1 < i < Q, is equal to the zeroset in B of two functions 
h and k. Note that M(0)oo C N{c) C B. Also recall that M{0)^ = N{c) f] CP^, that 
M = M{c) = N{c) \ CP^, and that Moo = M n CP^ = the set of limit points for e of 
solutions of (IIO.IOII with e 7^ 0. 

Because h and k vanish up to second order at p G M(0)oo, their second order Taylor 
expansions at p are canonically defined quadratic forms on Tp B, given by 

h^^\p) = (u?, «;)-(ciC'-2c2C + C3) e^ (10.16) 
k^^\p) = {id,Jiv)-{c,C-2cU + ce)e\ (10.17) 

in which w := (u + C,u — v. These formulas are obtained by replacing (u — v and e in p0.14j) 
and ()10.15p by their first order approximations w and e". 
The equations {u, u) = and {u, J u) + 1 = imply that 

det(M, Ju, u X Ju) = {u X Ju, u x Ju) = — (m, Jm)^ = —1, (10.18) 

and therefore the vectors u, Ju and u x Ju form a basis of C'^. Furthermore z = 1 implies 
that z = 0, and the equations in Lemma 110.51 for p & TpB imply that {u, w) = and 
{Ju,w) = 0, which in turn imply that w = 6u x Ju for some 5 G C. Also note that the 
condition that p ETp M(0)oo is equivalent to w = and 0, or 5 = ■£ = 0. 

Let a, b, S be functions of z, u, v, (, which together with e form a regular system of 
coordinates for B near p, in such a way that a = b = 6 = e = corresponds to the point p 
and, near p, the equations 6 = e = define M(0)oo- Then the tangent vector d/d6 at the 
origin corresponds to a tangent vector p oi N{c) at p, such that ^=0 and p is not tangent 
to M(0)oo- We can arrange this such that w = u x Ju. At 5 = e = the functions h and 
k and their first order derivatives with respect to 6 and e are equal to zero. Their Taylor 
expansions with respect to 6 and e start with quadratic terms of the following special form 

h^^'^ = hi{a, b) 6^ - h2{a, b) e\ k^^'^ = ki{a, b) 6^ - k2{a, b) e^. (10.19) 

The structure of the common zeroset of h and k in B will now be clarified in the following 
lemma. 

Lemma 10.6 Let h = h{a, b, S, e) and k = k{a, b, S, e) be two holomorphic functions de- 
fined in an open neighborhood of the origin in C^. Assume that their Taylor expansion 
at 6 = e = with respect to 6 and e start with quadratic terms as in ()10.19p . Write 
A(a, b) := hi{a, b) k2{a, b) — h2{a, b) ki{a, b), so that the equation A = means that the 
quadratic forms h^'^^ and A;'-^^ are proportional. If not both hi(0, 0) and /i2(0, 0) are equal 
to zero and not both ki{0, 0) and k2{0, 0) are equal to zero, then the origin can only be ap- 
proached by points in the common zeroset of h and k for which (5, e) 7^ (0, 0) ^/ A(0, 0) = 0. 

// conversely A(0, 0) = 0, hi{0, 0) 7^ 0, /i2(0, 0) 7^ and the derivative at (0, 0) of 
A(a, b) with respect to (a, b) is not equal to zero, then near (0, 0) the common zeroset of h 
and k is equal to the union of two smooth complex analytic surfaces which intersect cleanly 
along the smooth curve through the origin which is determined by the equations 6 = e = 0, 
A(a, b) = 0. 
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Proof Suppose that hi{0, 0) 7^ and /i2(0, 0) 7^ 0, which conditions are equivalent to 
the condition that h^'^^ is a nondegenerate quadratic form in S and e. (If fci(0, 0) 7^ 
and ^2(0, 0) 7^ 0, then we can interchange the roles of h and k.) Let 6 = 6{a, b) be a 
square root of h2{a, b)/hi{a, b) which depends holomorphically on (a, b) in a neighborhood 
of (0, 0). The Morse lemma with parameters, cf. Hormander fl^ Lemma 3.2.3], yields 
that there is a holomorphic change of the coordinates {6, e) to coordinates (x, y), depending 
holomorphically on (a, b), such that, near the origin, h = hi{a, b)xy. We can moreover 
arrange that in first order approximation at 5 = e = we have that x = 6 + 6{a, b) e and 
y = S — 9{a, b) e. 

The Taylor expansion of x ^— k{a, b, x, 0) at x = now starts with a quadratic term, 
which implies that we can write k{a, b, x, 0) = K{a, b, x) x^, in which K{a, b, x) is a holo- 
morphic function of (a, b, x) near the origin, and 

K{a, b, 0) := ki{a, h)/A - ^2(0, 6)/4^(a, bf = (/i2 ki - hi k2) /4/i2. 

For X 7^ the equation k{a, b, x) = is equivalent to the equation K{a, b, x) = 0, and we 
conclude that the point p cannot be a limit point of M when K{0, 0, 0) 7^ 0, or A(0, 0) 7^ 0. 

Assume conversely that K{0, 0, 0) = 0, which means that k2/ki = /i2/^i at (0, 0), and 
that the derivative at (0, 0) of (a, b) K{a, b, 0) is not equal to zero, which is equivalent 
to the condition that the derivative at (0, 0) of A is not equal to zero. For instance, assume 
that dK{0, b, 0)/db when 6 = 0. Then the implicit function theorem yields that there 
exists a holomorphic function B{a, x) of (a, x) near (0, 0) with B(0, 0) = 0, such that, 
for (a, b, x) near (0, 0, 0) the equation K{a, b, x) = is equivalent to 6 = B{a, x). This 
describes a smooth complex analytic surface, and we obtain the description of the common 
zeroset of h and k near the origin as in the lemma. 

The only case which we have not discussed yet is that A(0, 0) 7^ but not /ii(0, 0) 7^ 
and /i2(0, 0) 7^ and not A;i(0, 0) 7^ and ^2(0, 0) 7^ 0, for instance when hi{0, 0) 7^ 0, 
^2(0, 0) = 0, A;i(0, 0) = 0, and /c2(0, 0) 7^ 0. However, in this case we obtain, for an 
arbitrarily small positive constant c, that the points (a, b, 6, e) near the origin in the zeroset 
of h satisfy an estimate of the form \6\ < c\e\ and those in the zeroset of k satisfy |e| < c\S\, 
and the conclusion is that S = e = for the points (a, b, S, e) near the origin in the common 
zeroset of h and k. □ 

In our case the coefficients in (I1U.19|) are given by 

hi{a,b) = {uxJu,uxJu)=-l, /i2(a, 6) = Ci - 2c2 C + C3, flO 20) 

ki{a, b) = {u X J u, J {u X J u)), A;2(a, 6) = C4 — 2c5 C + ce. ' 

Until now we did not really use that M(c) = M with c and M as in Subsection 110.11 but 
from now on this assumption will be essential. Then /i2 = ci + C3 and k2 = c^^'^ + cq. If 
/i2 = then /c2 7^ because of the assumption that ci cq 7^ C3C4. It follows from Lemma H 0.61 
that the points of M(0)oo where h2 = ki = do not belong to M. 

Inserting /12 = Ci(^ + C3 and ^2 = 04^^ + cq in p().2()|l we obtain that 

A := h2 ki -hik2 = (ci ki + C4) + £3^1 + Cg, (10.21) 
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with ki = {u X Ju, J {u X J u)). Here the vectors u run over the eUiptic curve E given by 

(m, m) = 0, (m, Jm) + 1 = 0. (10.22) 

The u E E together with the free C, are parametrizing M(0)oo- The equation A = deter- 
mines a curve in the (u, C)-space. We have dA = at a zero of A if and only if A = 0, 
(ci fci + C4) C = and (ci C + C3) dfci = 0. 

If C 7^ 0, then Ci /ci + C4 = and A = yields that C3 ki + = 0. This leads to a 
contradiction with the assumption that Ci Cg 7^ C3 C4. 

If C = then dki = because ci + C3 = C3 7^ and A = yields that C3 fci + ce = 0. 
The tangent space of E is spanned by the vector u x J u, on which dki = if and only if 

= ((^i X J u) X Ju, J [u X Ju)) + {u X J {u X Ju), J [ux Ju)) 

= {u, J u) ■ {J u, J {u X J u)) — {Ju, J u) ■ {u, J {u X J u)) = —{Ju, J {u X Ju)), 

where in the first equality we used that J is symmetric, in the second that {u x J u) x J u = 
{u, Ju) Ju — {J u, Ju)u, and in the last that {u, J u) = —1 and once more that J is 
symmetric. A straightforward calculation shows that 

{Ju, J{ux J u)) = ( Ji - J2) ( J2 - -^3) (Js - Jl) Ui U2 Us, 

which, in view of the assumption that the Ji are different from each other, is equal to zero if 
and only if = for some i = 1, 2 or 3. Writing the indices of the coordinates of u modulo 
3, the condition that u E E now amounts to Ui+i^ + Ui+2' = 0, Jj+i Mj+i^ + Jj+2 + 1 = 0, 
or Mj+i^ = 1/ (Ji+2 - Ji+i), = 1/ (Ji+i - Ji+2). This implies that 

ki = {u X J U, J [u X J u)) = Ji (Ji+2 - Ji+l)'^ Uj+i^Mi+2^ = -Ju 

and we obtain a contradiction with the equation C3 /ci + Cg = 0, in view of the assumption 
that ce — C3 Jj 7^ 0. 

Applying Lemma 110.61 we obtain the conclusion that in the domain where ^ 7^ the 
curve Moo coincides with the subset of M(0)oo determined by the equation A = 0. It is 
smooth and near it M is equal to the union of two smooth complex analytic surfaces which 
intersect cleanly along M^o- 

If z = and u 7^ 0, then we work in a chart where, for some z = 1, 2 or 3, = 1. 
Then Vi = \ ui = \ and it becomes expedient to replace the functions h and k in ()10.14j) and 
(jlO.lSp near such a point p by 

h := Vi^ gi-2vig2 + 92. = {w, w) - {ciVi^ + cs) and flO 23) 

k := Vi^ gA-2vig^ + gQ = {w, J w) + {viZ - C)'^ - {dVi^ + Cq) e^, 

respectively, where w := ViU — v. The manifold M(0)oo near p is now parametrized with the 
curve of the {u, z) with Ui = 1, {u, u) = 0, {u, J u) + z'^ = 0, and the coordinate fj. 

At 2; = we have that {u, u) = and {u, Ju) = 0, which imply that {uxJu, uxJu) = 0. 
On the other hand it follows from Lemma 110.51 that {u, u) = 0, {u, v) = 0, {Ju, u) = 0, 
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{Ju, v) = 0, whereas ? and ( are free. It follows that u, u, v and therefore also w are 
multiples of u x J u, and we conclude that {w, w) = 0. If we let the vector d/d6 in the 
paragraph preceding (|1(J.19|) correspond this time to the vector p such that S = 0, u = 0, 
V = 0, ( = 1 and ^ = 0, then we obtain (|1U.19|) with hi = 0, /12 = Ciu,^ + C3, ki = 1, 
k2 = c^Vi^ + Cg. It follows from Lemma llU.61 that, even when k2 = 0, the point p can only 
be approached hj M if h2 = 0, which implies that ^2 7^ because of the assumption that 
Ci Cq 7^ C3 C4. Moreover, because Ci 7^ 0, the derivative of A = /i2 fci = Ci fj^ + C3 with respect 
to Vi is nonzero when A = 0. Again we can apply Lemma f 1 . 61 and conclude that the curve 
Moo is smooth at p, and that M near p is equal to the union of two smooth complex analytic 
surfaces which intersect cleanly along M^o- 

The case that z = and m = is treated by interchanging the role of the vectors {u, z) 
and (f , Q. In this case we use the assumption that C4 — ci Jj 7^ for every i. Again we 
can apply Lemma 110.61 and conclude that the curve M^o is smooth at p, and that M near p 
is equal to the union of two smooth complex analytic surfaces which intersect cleanly along 

The projection (u, Q ^ u exhibits M^c, as a branched covering over the elliptic curve 
where E is defined by ()10.22|) and C ^ C U {00} are the solutions of A = 0, with A as 
in ()10.21|) . Here = 00 corresponds to 2; = 0, in which case we interchange the role of the 
vectors (m, z) and (f , Q. The branching occurs when C3 /ci + Ce = or ci ki + = 0, and 
all these branch points are simple. A straightforward calculation shows that the equations 
u & E and fci + c = are equivalent to 

Ui^ = {J^ - c) / - Ji) {Ji - Ji+i) , I e Z/3Z. (10.24) 

It follows that there are 2 ■ 2'^ = 16 branch points, all of which are simple, because Jj — c 7^ 
for c = C4/C1 and for c = Cq/cs. 

The Riemann-Hurwitz formula says that if one has an n-fold branched covering from a 
curve r onto a curve C and B is the set of branch points in F, then 

genus(r) - 1 = n (genus(C) -l) + Yl o^der(6)/2, (10.25) 

beB 

cf. Farkas and Kra, ilS, p. 18]. Here the order of the branch point b is equal to m if the first 
m derivatives of the mapping at b are equal to zero. Because the genus of an elliptic curve 
is equal to one, it follows that the genus of M^c minus one is equal to 16/2 = 8, or the genus 
of Moo is equal to 9. We have proved: 

Proposition 10.7 Suppose that M(c) = M, with c and M as in Subsection 110. IL Then 
Moo := M n CP^ is determined by the condition that the quadratic forms h^'^^ and k^"^^ are 
proportional. Moo is a smooth closed algebraic curve in CP^ of genus equal to 9. Near Moo, 
the variety M is equal to the union of two smooth complex analytic surfaces which intersect 
cleanly along Moo- 
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The singularities of M can be resolved by considering the bundle G over CP^, of which 
the fiber Gp at p G CP^ consists of the space of all two-dimensional linear subspaces of the 
tangent space at p of CP®. Let Gm denote the restriction of G to M and let tm be the 
section of Gm which is obtained by assigning to p e M the tangent space TpM of M at 
p, which is regarded as an element of Gp. The projection n : tm M is an isomorphism 
from Tm onto M. Let M denote the closure of tm in the j^rojectiye variety G. Then M is 
a closed smooth two-dimensional subvariety of G. Define Moo '■= M \ TM. The projection 
TT : M — > M is an isomorphism from the complement tm of Moo in M, onto the complement 
M of Moo in M. On the other hand Moo is a smooth closed curve in M and the projection 
vr : Moo Moo is an unbranched two-fold covering. 

The mapping vr : M — ^ M is a so-called normalization of M, a regular mapping from an 
irreducible normal variety (every smooth variety is normal) onto M, which is a birational 
mapping and finite-to-one over every point of M, cf. |35s II. 5. 2]. Because normalizations are 
unique up to isomorphisms, one talks about the normalization of M. Our vr : M — > M is a 
simple, explicit one. 

The Riemann-Hurwitz formula (jl(J.25|) yields that the genus of Moo minus one is equal 
to 2 ■ (9 — 1) = 16, or that the genus of Moo is equal to 17. In this way we obtain a smooth 
completion M of M which is obtained by adding a smooth curve of genus 17 at infinity. 
Here the word "completion" is used in the algebraic sense. Proposition 110.81 below says that 
it can also be used in the sense that the flows of ^ and rj, with complex times, are complete 
on M in the sense that they define a transitive action on M/S of the additive group C^. 

Proposition 10.8 Suppose that M(c) = M, with c and M as in Subsection 110. IL Let M 
be the smooth completion of M described above, the normalization of the projective clusure 
of M , which is obtained by adding to M a smooth curve of genus 17 at infinity. Then the 
rational vector fields ^ and rj on M are everywhere finite and linearly independent. Their 
respective flows e*^ and e^^ with complex times t and s define a transitive action of the 
additive group on M, and for each p & M the mapping (t, s) ^— > e*^ oe^'^{p) defines an 
isomorphism from the complex torus C^/A onto M. Here 

A := {(s, t) I e'^oe''%p) = p} 

denotes the period lattice. It does not depend on the choice of p and has a Z-basis consisting 
of four elements o/C^ ~ which are linearly independent over R. 

Proof We first investigate the vector fields ()9.6|) — ()9.13p near infinity when z ^ 0, where 
we use projective coordinates with z = 1. With e as the last coordinate, this means that we 
identify the affine coordinates ((m, i), (y, (), 1) with p = ((m, 1), (t>, (), e), where 

u = z'-^ u, V = z~^ V, ( = z'^ C, € = z~'^, or (10.26) 
u = e-^u,v = e-^ v,C = e"^ C,z = (10.27) 
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Note that our notation means that we have to put tilde's over all the coordinates in the 
formulas (|9.6|) — (j9.13|) for the vector fields ^ and rj in the affine coordinate system. 

We will write the point p, at which we consider the vector fields ^ and 77, as an analytic 
function of e and a base point po which varies in the curve at infinity. It follows that we have 
a convergent power series expansion p = ^j>o Pjy which the coefficients pj for j > 1 
depend analytically on the point po in the curve at infinity. We may also assume that the 
vector pi is not tangent to the curve at infinity, which means that it can be identified with 
the vector p = d/d6 in the paragraph preceding (jl(J.19j) . Recall also that 

Cu + Cu-v = w = 9ux Ju, (10.28) 

in which the nonzero factor 6 is equal to a square root of — — 1. Note that there the coordi- 
nates of the base point po are denoted by {{u, 1), (f , (), 0), instead of the {{uq, 1), (fo, Co), 0) 
which we will use here. 

With these notations, we have that 

= —z~'^ i'z = — (m X J m, Jv) = —e^^{ux J u, J v). 

The constant term in the expression following is equal to zero, because Vq = (0 Uq. The 
first order term in its Taylor expansion with respect to e is equal to 

{Ui X Juo, Jvq) + (mo X J-Ui, Jvo) + {uq X J Uq, J Vi) = {uq X J Uq, J {vi "Co'^o)) 

where we again have used that Vq = Co ^o- Using (|1U.28|) . we obtain that attains the finite 
value 

= {u X J u, J id) at infinity. (10.29) 
Using that u = eu, we subsequently obtain that 

^u= {^e)u + e^u = r^ {e{^e)u + u x J (f -(«)), 

where we have used the homogeneity of ^ of degree 3 and z = 1. The constant term in the 
expression following is equal to zero, because Vq = (qUq. Using (jl0.29p . we obtain that 
the first order term in its Taylor expansion with respect to e is equal to 

{u X J u, J w) u + Ui X J (v — (u) + u X J {vi — (lU — ( Ui) = {u X J u, J w) u — u X Jw, 

where we have dropped all the subscripts in the notation. The inner product of this 
expression with u is equal to zero. Using that {u, Ju) = 1, we obtain that the inner product 
with Ju is equal to zero as well. Finally the inner product with u x Ju is equal to 

— {u X Jw, u X J u) = {u X {u X Ju), J w) = —{u, Jw) 

because ux{ux J u) = {u, Ju) u— {u, u) J u, {u, J u) = —1 and {u, u) = 0. Now it follows 
from Lemma flO.51 with z = 1, 'z = 0, and {Ju, u) = —1 that 

{Ju, w) = ( {J u, u) + ( {J u, u) — {J u, v) = 0. 
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Because u, Ju and u x Ju form a basis of C^, the conclusion is that the e~^-term in the 
expansion of in powers of e is equal to zero as well, or that is finite. 

Using that ^ is tangent to the surface M, we have obtained sufficient evidence to conclude 
that C, is finite in the complement of at most finitely many points of the curve Moo- In 
combination with the rationality of ^ this implies that is finite on M. 

For the vector field rj we begin with 

rje = —e^-^ {u x J u, v) . 

The constant term in the expression following is equal to zero, because vq = (qUq. The 
first order term in its Taylor expansion with respect to e is equal to 

{Ui X J Uq, Vq) + {uqX J Ui, Vq) + {uq X J Uq, Vi) = {uq X J Uq, Vi - (qUq) 

where we again have used that Vq = (qUq. Using (jlU.28|) . we obtain that rje attains the finite 
value 

rje = {u X J u, w) at infinity. (10.30) 

Note that w = 9ux Ju for a nonzero factor 9, and that {u x Ju, u x Ju) = —1, cf. ()10.18|1 . 
Therefore rje = —6^0at every point on the curve at infinity where z 0. 
Using that u = eu, we subsequently obtain that 

rju = (rje) u + erfu = (e {rje) u + u x v) , 

where we have used the homogeneity of rj of degree 3 and z = 1. The constant term in the 
expression following is equal to zero, because vq = (qUq. Using p().3()|l . we obtain that 
the first order term in its Taylor expansion with respect to e is equal to 

{u X Ju, w) U + Ui X V + U X Vi = {u X Ju, w) u — U X w, 

where we have dropped all the subscripts in the notation. The inner product of this 
expression with u is equal to zero. Using that {u, Ju) = 1, we obtain that the inner 
product with Jm is equal to zero as well. Finally the inner product with u x Ju is equal to 
— {u X w, u X Ju) =0, because w = 6 u x J u. Again using that u, Ju and u x Ju form 
a basis of C^, we obtain that the e~^-term in the expansion of rju in powers of e is equal to 
zero as well, or that t]u is finite. In the same way as for ^, we conclude that the vector field 
1] is finite on M. 

Let S denote the set of points in M where ^ and rj are linearly dependent. Proposition 
110.21 implies that fl M = 0, which means that S is contained in the curve M^o at infinity. 
Because ^ and rj commute, the set S is invariant under the fiow of both vector fields, and it 
follows that at every point of S both vector fields must be tangent to the curve at infinity. 
Because rje ^ at every point of the curve at infinity where z ^ 0, we are left with the 
points at infinity where z = 0. 

If z = 0, then we have (u, u) = and {u, Ju) = and it would follow that -u = if 
Ui = for some i. Therefore, assuming that m 7^ 0, we have for every i that Ui 7^ 0. In 
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the projective coordinate chart where Mj = 1 we have that e = u~ . It turns out that then 
r^e = at e = 0, which means that rj is tangent to the curve at infinity when z = 0. For this 
reason we turn to the computation of ^e, which is equal to the i-th coordinate of 

— = —e~^ u X (z J V — ( J u). 

Note that Vq = Xuq and Co = -^^o for the same factor A, and Zq = 0, which imphes that 
Co = as well. Therefore the constant term in the expression after is equal to zero and 
the first order term in its Taylor expansion with respect to e is equal to 

(^A — () u X J u, 

where we have dropped all the subscripts in the notation. Because X = Vi when Ui = 1, 
we conclude from ()1().23|1 that the factor SX — ( is not equal to zero, where we also use that 
the equations in Lemma 110.51 with z = imply that (w, Jw)=0 when w = ViU + ViU — v. 
Because m x Jm 7^ 0, it follows that for at least one choice of i we obtain that 7^ 0, which 
proves that ^ is not tangent to the curve at infinity when z = 0. 

The case that z = and m = is treated by interchanging the role of the vectors {u, z) 
and [v, (). Collecting all results, we have proved that S* = 0, or that C, and rj are linearly 
independent at every point of M. 

Using the branched covering over Uc in Subsection 111.11 one obtains that the complex 
level surface M is connected (in contrast to the real one), and therefore M is connected as 
well. The remaining conclusions of the proposition now follow by applying the argument of 
Arnol'd and Avez [3 Appendix 26] as at the end of Section [7| □ 

Remark 10.9 In the complex time coordinates on M ~ C^/A, the vector fields ^ and rj 
are constant (and linearly independent). Proposition 110. HI implies that the rotational motion 
of Chaplygin's sphere with horizontal moment is algebraically integrable according to the 
definition of Adler and van Moerbeke [21 p. 297]. In view of Subsection 19.21 this result 
remains true for arbitrary non- vertical moment. 

A very different proof of the algebraic integrability can be given by means of Chaplygin's 
integration of the system in terms of hyperelliptic integrals as described in Subsections 111.21 
and ITOl See Subsection im 

Remark 10.10 The surface M is invariant under the antipodal mapping x ^ —x. In 
projective coordinates near infinity, where we take one of the affine coordinates equal to 
1, this mapping is given by e 1^ — e, keeping the affine coordinates fixed. The set where 
e = is the projective space at infinity, which belongs to the fixed point set of the antipodal 
mapping. The coordinates for CP*/±1 near infinity are obtained by replacing e by e^. The 
antipodal mapping interchanges the two sheets along Moo, and it follows that M/±l is a 
smooth variety. Its curve at infinity, (M/±l) \ (M/±l), is isomorphic to M^. 

The antipodal mapping extends to an involution in M without fixed points, which leaves 
the vector fields ^ and 77 invariant. It follows that the projection from M to M/±l is a 
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twofold unbranched covering, which intertwines ^ and t] with two vector fields on M/±l, 
which we also denote by ^ and rj, which at every point are regular and linearly independent. 
Therefore the complex times of the fiows of ^ and t] lead to an identification of M/±l with 
a complex torus, on which the vector fields ^ and 7] are constant. 

Because the antipodal mapping belongs to the group S in ()10.3H) . we obtain an eight- 
fold unbranched covering vr : M/±l M/S such that the projection from M onto M/S 
in Proposition 110.111 below is equal to the composition of the twofold covering from M 
onto M/±l, followed by vr. In this way the curve Mqo/S of genus two is isomorphic to 
Moo/(S/±l). 



10.3 A Discrete Symmetry Group 

Let S denote the group of the 16 transformations 5* in of the form 

S{{u, z), {v, 0) = {{e,Ru, €2 z) , {e^Rv, ei 0) , (10.31) 

in which = ±1 and R G S0(3) is a diagonal matrix, with ±1 on the diagonal, two of 
them equal to —1 if i? 7^ 1. A straighforward computation show that every 5* G S leaves 
the functions in ()9.14|) invariant, and therefore leaves M invariant as well, for any choice 
of the constants q. Moreover, every S* G S alos leaves both vector fields ^ and rj invariant. 
Each linear transformations S has a natural extension to a projective linear transformation 
of CP*, which leaves M invariant. It also has a natural extension to the bundle G mentioned 
after Proposition 111). 81 and this extension leaves the smooth variety M and the vector fields 
^ and rj on it invariant. 

Proposition 10.11 Suppose that M(c) = M, with c and M as in Subs ection MO. 11 If S E 

and S 1 then S has no fixed points in M . 

As a consequence, the quotient M/S is a smooth complex projective algebraic surface. 
The projection from M onto M/S intertwines the vector fields ^ and rj with vector fields on 
M/S which we denote by the same symbols. The vector fields ^ and rj on M/S are regular 
and linearly independent at every point, and therefore M/S is isomorphic to a complex torus 
as well. 

Under the projection from M onto M/S, the curve of genus 17 is mapped onto a 
smooth curve Mqo/S of genus equal to 2. 

Proof Let F denote the set of fixed points of 5*. Because the vector fields ^ and r] are 
invariant under F is invariant under the fiows of ^ and rj with complex times. Because 
these fiows define a transitive action of on M, it follows that F is either void or equal to 
M. Because it is easily verified that M is not contained in F, the conclusion is that S has 
no fixed points in M. 

The restriction to P := Mqo of the projection from M onto M/S defines a 16- fold un- 
branched covering map from P onto C := Mqo/S. The Riemann-Hurwitz formula 110.251 



65 



therefore yields that 16 (genus(C) — 1) = genus(r) — 1 = 17 — 1 = 16, which imphes that 
genus(C) — 1 = 1, or the genus of C is equal to 2. □ 



Remark 10.12 Because every curve of genus 2 is hyperelliptic, cf. Farkas and Kra 
|131 Prop. III.7.2], we conclude that by adding a hyperelliptic curve of genus 2 at infinity, 
the manifold M/S can be completed to a complex torus, on which ^ and r] are linearly 
independent and constant vector fields. 

The torus M/S is isomorphic to the Jacobi variety of the hyperelliptic curve C which 
appears in Chaplygin's integration by means of hyperelliptic integrals. See Remark 111. 71 

Remark 10.13 Let Mqo/S denote the hyperelliptic curve of genus 2 which is added to 
M/E at infinity in order to obtain the torus M/S as the completion of M/S, cf. Remark 
110.121 Let J™ = diag(J2-^3, J?,Ji^ J1J2) be the comatrix of J as defined in ()9.30|1 . The 
rational function 

Moo 9 ((m, z), {v, C)) ^ -(^, r°u)/{u, Ju) 

induces a twofold branched covering from Mqo/S onto CP^, which branches over the points 
X = Ji (corresponding to Ui = 0) for z = 1, 2, 3, A = 00 (corresponding to 2; = 0), and 
the two zeros of the polynomial p{X) given by ()9.26j) . For the role of p(A), see also iv) in 
Subsection llO.il or ()10.33|) where the values bi of the functions hi are given in terms of the 
values Ci of the functions fi by means of ()9.41|1 — ()9.44j) . It follows that Mqo/S is isomorphic 
to the hyperelliptic curve which is defined by the equation 

= p(A) n - A) 

i=l 

between the projective coordinates (A, /i) in CP^. 

Of the six fixed points of the hyperelliptic involution (A, fi) t— > (A, — yu), the four corre- 
sponding to A = Ji, J2, J3, cxo do not depend on the values q of the functions whereas 
the other two, the zeros of p(A), move freely with the q, even with the constants of motion T 
and j of Chaplygin's sphere. This means that the curves Mqo/S are non-isomorphic for the 
generic variation of the constants of motion, and describe a two-dimensional subvariety of 
the three-dimensional moduli space of curves of genus two. If we are also varj^the constants 
Ji freely, then there is no restriction on the isomorphism class of the curve Mqo/S. 

Remark 111.81 contains an explicit verification that the curve Mqo/S is isomorphic to the 
hyperelliptic curve C introduced in (jll.SOp . 

Question 10.14 As observed in Remark ll().12^ the torus M/S is isomorphic to the Jacobi 
variety of the hyperelliptic curve C. According to Remark 110. 13^ C is isomorphic to the 
curve Moo/S which is added at infinity to the affine algebraic surface M/S in order to 
obtain the toral completion M/S. It follows from Matsusaka j2Zl that M/S is isomorphic 
to the Jacobi variety of the curve Mqo/S, and that Mqo/S is canonically embedded in its 
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Jacobi variety M /Ti, if and only if the self-intersection number of the curve Mqo/S in M/S 
is equal to two. (I owe this reference to Ben Moonen.) Is it possible to verify directly that 
the self-intersection number of the curve M^/H in M/S is equal to two? 

Remark 10.15 In terms of the parametrization of M by means of the complex times of the 
flows of the vector fields ^ and rj, cf. Proposition 110. 8( the condition that 5* commutes with 
these flows implies that S* is a translation. For each S* G S we have that S*^ = 1, cf. p0.31|) . 
Therefore, if we provide ^ with a real basis with respect to which the period lattice 
A is equal to Z^, we obtain that S is equal to a translation over a vector v E (|Z) /Z^. 
Because there are 2^ = 16 such vectors v, it follows that S is equal to the group of all 
translations of order two in the torus M. In other words, the covering of M M/S of the 
complex torus M/S is obtained by replacing the period lattice Aq of M/S by A = 2Ao. 

The group S is closely related to the set oftheta characteristics, as discussed in Mumford 
EH p. 163]. 

I owe this remark, and the encouragement that Proposition I1U.8I and Proposition 110.111 
might be the right picture of the projective completion of M, to Frans Oort. 



10.4 Jordan Rizov's answer to Question 110.141 

The following answer to Question 110. 141 has been kindly provided to me by Jordan Rizov. 
Let us collect in i) - vi) below the abstract data we shall be working with. 

i) Following Chapter El we consider the affine variety M C defined by the equations 
()9.14|) . in which /i, 1 < « < 6 are constants. According to Proposition 110. l| M is a 
nonsingular two-dimensional complex variety. 

ii) Consider the "standard" embedding C CP® and let M be the projective closure 
of M in CP® with respect to the complex topology. Then M is also the projective 
closure of M in the Zariski topolgy, because M is defined as the zeroset of polynomials. 
According to Proposition 110. 7( Moo := M \ M is a smooth algebraic curve of genus 9 
and M is singular along Moo- 

iii) Consider the normalization vr : M — M of M as constructed after Proposition 110. 7( 
where M is a nonsingular two-dimensional projective variety. The preimage M^o of 
Moo in M is a nonsingular projective curve and 

71 : Moo ^ Moo 

is an unramified two-fold covering. Hence, by the Riemann-Hurwitz theorem, the genus 
of Moo is 17. 

iv) Proposition llO.SI savs that M is a two-dimensional complex torus, and hence an Abelian 
variety (it is projective). 
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v) As described in the beginning of Subsection I1U.31 there is a group S of order 16 
acting on such that its action extends to an action on M and on M. Furthermore 
Proposition 110. ill savs that S acts freely on M. Hence, by |3m Ch. II, §7, Thm. 1], 
the quotient map 

TTs : M ^ M/S 

is etale (= an unramified covering map). Moreover, because the action is free, E acts 
a a finite group of translations on the Abelian surface M and by Ch. II, §7, Thm. 
4] the quotient M/E is an Abelian variety. 

vi) Since vr^ is etale, the nonsingular curve Moo is mapped onto a nonsignular complete 
curve Moo/S of genus 2. 

Before going on with any computations, let us simplify the notations a little bit by putting 

r := Moo/E, 
S := M/S, 
Ks '■= the canonical class of S, 

i.e. the divisor class of a top degree differential form. The question posed at the end of 
Question ll().14l is whether one can compute directly the self-instersection of F on S. We will 
do this using the 

Adjunction formula Let T be a nonsingular curve of genus g-^ on a nonsingular surface 
S with canonical class Ks- Then the following relation holds 

2g^-2 = V-{V + Ks). 

Proof The proof and the construction of the intersection pairing on a nonsingular surface, 
in an "algebraic" way, can be found in ^3 Ch. 5, §1], where Proposition 1.5 is the adjunction 
formula. An "analytic" proof is given in Ch. 4, §1]. □ 

Therefore, in order to compute the self-intersection F ■ F one has to enquire a little bit 
about the canicial class Ks of S. As we already saw, S is an Abelian surface, and the next 
result gives all we need. 

Fact // A is an Abelian variety of dimension g, then 

^'a - Oa, 

or equivalently, the canonical class Ka of A is trivial. 

Proof See for instance [TU Ch. 1, Prop. 1.5] or [301 Ch. 1, (5)], especially (*) on page 4, 
and ,30- Ch. 2, §4, (4) on p. 42]. □ 

The computation Applying the adjunction formula to S and F with g^^^^^^ = 2 and Ks 
trivial, one gets 

F-F = 2x2-2 = 2, 
which yields that the self-intersection number of F in S" is equal to two. 
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10.5 The system in = A 

Let hi be the functions on given by — In this subsection we will investigate 
the level set L = L{h) defined by the equations hi = bi, /i2 = 0, h^ = 63, /14 = 64. We will 
assume that the constants bi, 63, 64 satisfy the following conditions: 

61 ^ 0, ((trace J) bi + b^f ^ bi 63 and (10.32) 
p(J.) =61 J.2_ ((trace J) 61 + 64) Ji + 63^0 fori = l, 2, 3. (10.33) 

Let U denote the set of ((m, z), {v, ()) G such that {u, z) and {v, Q are linearly 
independent, U is an open subset of C®. Let / denote the mapping from U to defined 
by the functions /j in (j9.14|) . Let V be the smooth hypersurface in \ {0} defined by the 
equation /i2 = (g, r) = 0, and let H : V ^ denote the mapping defined by Hi = hi, 
H2 = hs and H3 = (trace J) hi + h^, in which the functions hi are given by ()9.36p — ()9.39|) . 
Let K be the mapping from to C'^ defined by 

K (c) = (ci C3 - C2^ C4 C6 - C5^ Ci C6 + C3 C4 - 2C2 C5) . 

Then the equations (lOijl mean that A (C*^) cVandHoA = Kof. 

Let c G be such that 

(61, 63, (trace J) bi + 64) = K{c). (10.34) 

Write F = ( and G = ( . Then the assumptions (fTim . (ITTOH|l just 
V ^2 C3 y \ C5 cq J 

mean that the polynomial p : A det(G — XF) given by ()9.28|) is of second order and has 

two distinct zeros, none of these equal to one of the Jj's. In other words, c satisfies the 

assumptions in Subsection 110.11 Note that these conditions imply that F and G are linearly 

independent, which in turn imply that M(c) C U and that the Jacobi matrix TcK at the 

point c of is surjective. Let x G M(c), which means that f{x) = c. Proposition 110. ll 

implies that the Jacobi matrix f at the point x of / is surjective. Write y = A{x). Then 

H{y) = H{A{x)) = K{f{x)) = K{c). The chain rule implies that Ty if o T^. A = TeKoT^/, 

which is surjective, and therefore Ty H is surjective as well. Using also that A intertwines 

the vector fields ^ and 77 in with the vector fields ^ and rj in e(3) ~ C^, as we have seen 

in Subsection 19.31 we have proved: 

Proposition 10.16 L = L{b) is a smooth two-dimensional affine subvariety of C®. // 
()10.34p holds then A\m{c) defines a twofold unbranched covering from M(c) onto L{b). It 
intertwines the vector fields ^ and rj on M(c) with the vector fields ^ and t] in . The 
latter vector fields are tangent to L{b) and linearly independent at every point of L{b). The 
mapping A|jv/(c) induces a birational isomorphism from M(c)/±1 onto L{b), which we will 
also denote by A. 

The statements that L{b) is smooth and ^ and r] are linearly independent at every point 
of L{b) can also be checked directly, but we found the proof which uses the system in 
simpler. 



69 



That the inverse of the rational map A : M(c)/±1 L{h) is rational follows from the 
general fact that if / : X — > F is a rational map between irreducible varieties X and Y 
of the same dimension, f{X) is dense in Y and / is injective over the preimage of a dense 
subset of y, then / is a birational isomorphism. Indeed, the homomorphism /* from the 
field C(y) of rational functions on Y to the field C(X) of rational functions on X is injective 
because f{X) is dense in Y . Furthermore the degree of the field extension of /*C(y) by 
C(X) is equal to the number of the elements of the generic fiber of /, cf. II.5.2.Thm. 7] 
(in characteristic zero every field extension is separable). In our case this implies that /* is 
surjective. Clearly / has a rational inverse if and only /* : C{y) — > C(X) is an isomorphism. 

According to Remark 110.101 the manifold M(c)/±1 is a complex torus to which the 
vector fields ^ and rj on M(c)/±1 extend as constant vector fields. Therefore M(c)/±1 can 
be viewed as a toral completion of L{h) to which the vector fields ^ and r] on L{h) extend 
as constant vector fields. This completion is obtained by adding a curve to infinity which is 
isomorphic to the smooth curve M(c)oo of genus 9. 

Proposition 10.17 Let M(c) and L{b) denote the closure of M{c) and L{b) in CP® and 
CP^, respectively. Assume that p0.34|) holds. Then A extends by continuity to a finite 
morphism from M(c) onto L{b), which factorizes through a morphism from M(c)/±1 onto 
L{b), which we also denote by A. M(c)/±1 is the complex torus of B.emark 110.101 and 
A : M(c)/±1 — > L{b) is a normalization of L{b). 

The restriction of A t o (M (c)/±1) \ (M(c)/±1) = M(c) \ M(c) = M(c)oo maps M(c)oo 
onto the curve L{b)oo '■= L{b) \ L{b) of L{b) at infinity. It assigns to the one- dimensional 
linear subspace 

C((m, z),iv, 0, 0)gM(c)oo 
of the one- dimensional linear subspace 

C {zu, J u X u, 0) 

o/C'^, where {z u, J u x u, 0) has to be replaced by {(v, Jv x v, 0) when u = and z = 0. 
The image L{b)oo is a smooth elliptic curve in CP^. The mapping A : M(c)oo L{b)oo is 
a twofold branched covering of the curve M{c)oo of genus 9 over the elliptic curve L{b) oo, 
where the 16 branch points in M(c)oo coincide with the branch points mentioned in the text 
preceding Proposition II . 71 

Proof The mapping A : C® — extends to a homogeneous polynomial mapping A : 
of degree two by means of the formula 

A((m, z), {v, C), e) = {u X V, zv - (u, e^). 

Near infinity, where e = 0, we can, as in the proof of Proposition ll().8| write p = {{u, z), {v, ()) 
as a convergent power series p = J2k>o Vk, in which C (po, 0) varies over the curve M(c)oo 
and the coefficients pk with k >1 depend analytically on p^. Because {uq, zq) and (f o, Co) are 
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linearly dependent when C {po, 0) G M(c)oo, we have that A (po) = 0. If the i-th coordinate 
of 

A' := D A (po) Pi = (ui xvo + UoAvi,ZiVo + ZqVi - Ci^o - Co^i) (10.35) 

is nonzero, then a division of all the other coordinates of A{p) by A{p)i yields the coordinates 
of A{p) in the standard projective coordinates in which the i-th coordinates is kept equal to 
1. Because A{p)i = ae in which a has a nonzero limit as e — 0, we obtain that A{p{e), e) 
converges in these coordinates as e — > 0. Moreover, its last coordinate e^/ A {p)i converges 
to zero as e — ^ 0, which means that the limit point belongs to the projective space CP^ = 
\ at infinity. 

The proof of Proposition 110.81 yields that the r-component of the vector A' in p0.35|) is 
nonzero whenever {po, 0) G M(c)oo- When zq ^ 0, we can work in the projective coordinate 
system where ^ = 1, hence zq = \^ z\ = 0, in which case the second component of A' is 
equal to the vector r := v\ — C,\Uq — C,qV\. In the proof of Proposition 110.81 we obtained 
that r = 6uq X Juq for some nonzero factor 6. Using that Vq = (qUq, we find that the 
p-component of A' then is equal to 

Uq X (t>i — (qVi) = Uq X r = 6uq x {uq x Juq) = —Ouq, 

because {uq, Uq) = and {uq, J Uq) + 1 = 0. 

If 2;o = and Mq 7^ 0, then we can work in the projective coordinate system where one 
of the coordinates Mq.i of "^o is identically equal to 1. We have Co = and Vq = fo,jMo; 
and therefore r = {zi Vqj — Ci) Uq, where in the proof of Porposition 110.81 we obtained that 
^1 "^0,1 ~ Ci 7^ 0. If Zq and Uq = 0, then we interchange the roles of the vectors {uq, zq) and 

(^^0, Co)- 

This concludes the proof that A has a continuous extension A to M(c) which maps M(c)oo 
into the projective space CP^ ~ CP^ at infinity. Furthermore, on the dense subset of M(c)oo 
where z 7^ it matches the description in Proposition 110.1 71 which therefore is valid at every 
point of M(c)oo- We also obtain that for 2; = 1 the restriction of A to M(c)oo is equal to the 
composition of the projection onto u, followed by the embedding u C {u, J u x u), where 
u runs over the elliptic curve E defined by {u, u) = 0, (m, Jm) + 1 = 0. This shows that 
A (M(c)oo) is isomorphic to E and that the restriction of A to M(c)oo is a twofold branched 
covering with the branch points as mentioned in the text preceding Proposition 110.71 (It is 
easy to verify that the points on M(c)oo with 2; = are no branch points.) 

The continuity of A, together with A(M(c)) = L{h) implies that L{h) C A (^M(c) j C L{h). 
On the other hand, because M(c) is compact, the continuity of A also implies that A (c)^ 
is a compact, hence closed subset of CP^. We conclude that A ^M(c) j is equal to the closure 

Uh) of L{h) in CP^ 

The graph of A : M(c)/±1 — > L{h) is equal to the projective closure of the graph of 
A : M(c)/±1 L{b), where the latter graph is an affine algebraic variety. It follows that 
the graph of A : M(c)/±1 L{b) is an algebraic variety, cf. Lojasiewicz jJEl P- 383], which 
implies that A : M(c)/±1 L{b) is an algebraic morphism. Because it is everywhere finte 
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and a birational isomorphism from M(c) onto L{b), it is a normalization of L{b). Note that 
M(c)/±1 is normal, because it is smooth. □ 



Proposition 10.18 We have L{b)oo = L{0)oo for every choice of the constants hi in the 
equations hi = hi which define L{b). L{0)oo is a smooth elliptic curve in the projective space 
CP^ = CP6 \ ~ CP5 at znfimty. 



Proof Let S be the additional projective coordinate for CP^ such that 6 = corresponds to 
the projective space at infinity. (We had 6 = in the proof of Proposition 110.1 71 ) Then the 
equations hi = hi for L{b) correspond to the homogenized equations hi — bi6'^ = 0. It follows 
that L{b) oo is contained in the subvariety L(0)oo of CP^ ~ CP^, which is determined by 
the equations 6 = and hi = 0. 

It follows from ()9.36p . ()9.38|) and ()9.39|) that the equations hi = 0, h^ = 0, h^ = are 

(l3 



three independent hnear equations for gi^, 



^ which have the solutions 



Qi 



(10.36) 



Here the index i is counted modulo 3 (cyclic notation). This determines the qt up to their 
signs in terms of r. Also note that r = implies that g = 0. because (g, r) = (0, 0) is 
excluded for the projective space at infinity, we have always that q ^ 0. 

Let ip be equal to minus the product of qi ri + g2 + qs r^, —qi ri + g2 ^2 + qs rs, qi ri — 
0.2 r2 + 93 and qi ri + g2 ^^2 — % r^. Then, for given r, there exists a solution q of (jl0.36|) 
and /i2 = if and only if there exists a solution q of (jl0.36|) and ip = 0. Two of the eight 
sign choices for the coordinates g^, then lead to two solution g and — g of (g, r) = 0. 

On the other hand 



; 4 4 o 2 2 2 2 

W = 2^ qi Ti -2 qi qi+i r^+i . 

Substituting (|10.36p in (|10.37p . we obtain after a straightforward calculation that 



(10.37) 



-1 2 



^ {Ji - Ji+lf Ti^ Ti+i' 
ieZ/3Z 



iGZ/3Z 



and therefore the equation ip = is equivalent to 



(Ji - J.+if ri\i+i^ = 0. 

jez/3Z 



;i0.38) 



Let F be the curve in CP^ defined by (jl0.38|) . The mapping r h-* (ri^, r2^, r^"^) defines 
a 16-fold branched covering of F over a nondegenerate quadric in CP^, which is isomorphic 
to CP^. The branching occurs when one of the coordinates of r is equal to zero, in which 
case another coordinate of r has to be equal to zero as well. Therefore the branching occurs 
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at the three coordinate axes, where the sheets of the covering are connected to each other. 
At these points, for instance r = 63, gi = 1/ (Ji — J2), q2 = —Qi, = 0, a. straightforward 
check shows that the derivatives of the functions hi defined by (|9.36|) . (j9.38|) and (j9.39p are 
hnearly independent. Therefore, although F is singular (has ordinary double points) at the 
coordinate axes, the curve M(0)oo is smooth at the corresponding points (g, r). 

For each r G -F we have two opposite g's which satisfy ()10.36|) and (g, r) = 0. The 
equations (jlU.36p have the solution g = if and only if 

ri+i' = r,_^\ t e Z/3Z, (10.39) 

and it is easily verified that these equations imply p0.38|) . Therefore the projection (g, r) 1-^ 
r defines a twofold covering from L(0)oo onto F, which is branched at the four points in F 
defined by (jl(J.39|) . These are smooth points of F. If ?/ is a local analytic coordinate of F 
near such a point, and we substitute y = z'^, then we obtain that the corresponding points 
(g, r) G L(0)oo can be written as q = zu (-z^)^, r = r (z^), in which u{y) and r{y) are analytic 
functions of y and m(0) 7^ 0. It follows that L{0)oo is smooth at (0, r{0)), and that this point 
is a simple branch point for the covering (g, r) i-^ r. We conclude that L{0)oq is smooth and 
connected, and therefore irreducible. Because the curve L{b)oo is a component of L(0)oo, 
it follows that L{h)oo = L{Q)ao- We know already from Proposition 110.171 that L{h)oo is an 
elliptic curve, but the above description can be also be used to verify directly that the curve 
L(0)oo is elliptic. □ 



Remark 10.19 At first sight the fact that the curve L{b)oo in L{b) at infinity does not 

depend on the values of bi, 63, 64 is quite disturbing.^ According to Remark 110.101 the 

quotient of M(c)oo by the group E/±l is isomorphic to Moo/E. According to Remark [10.131 
its isomorphism class varies in a two-dimensional subvariety of the three-dimensional moduli 
space of curves of genus two. As a consequence the curves M(c)oo in general will not be 
isomorphic either as we vary the constants of motion. The question is where these moduli 
appear in the completion of L{b), if the curve at infinity of the projective closure L{b) is the 
same for all b. 

The answer is that L{b)oo = L{0)oq is an ordinary double curve of L{b) at all points except 
the branch points of the twofold covering A : M(c)oo -^(O)oo- At these branch points, the 
variety L{b) has worse singularities. (We conjecture that, as in Mumford's appendix to [21, 
these are ordinary pinch points, where L{b) has local analytic equations = y z"^.) Accord- 
ing to Proposition 110.171 the branch points in L(0)oo of A : M(c)oo — > L{0)oo correspond to 
points (g, r) G CP^ for which p0.24|) holds with u replaced by g. Here c runs over the two 
solutions of the equation 

61 — ((trace J) bi + 64) c + 63 = 0, 

cf. fll0.33j) . and therefore the branch points move as a function of the moduli. 

The situation is very much similar to the description of F^ in Mumford's appendix to 
in the text starting with "Thus C is an ordinary double curve of F^ ..." and ending with 
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"... , hence C has genus 9", on p. 330 and 331. One difference is that our normahzation 
A : M(c)/±1 L{b) is a quite simple, concrete one, whereas Mumford's normahzation 
TT : Fc — > -Fc is abstract. 

Another difference is that our normahzation is equal to an 8-fold unbranched covering 
of the Jacobi variety of a hyperelliptic curve of genus 2, a characterization which does not 
appear in The symmetry group is the group S/±l, with E as in ()10.31|) . which on the 
(g, r)-space acts by means of the transformations 

S{q, r) = {eRq, Rr), 

in which e = ±1 and i? is a diagonal rotation as in (|l(J.3ip . 

If we take L{b) as defined by hi = hi with hi ^ 0, then S*hi = hi when i 7^ 2, but 
^2 = {q, r) satisfies S*h2 = — /i2- If ^2 7^ 0, we therefore can only divide out the subgroup 
of four elements S for which e = 1. According to the Riemann-Hurwitz formula ()10.25p . the 
quotient of the curve M(c) 00 of genus 9 by this group of four elements has genus equal to 3. 
The possibility of arriving at a curve of genus 2 may therefore be related to the fact that we 
restricted ourselves to the hypersurface {q, r) = 0. 

Question 10.20 What happens with the sytem on the surface hi = hi when 62 7^ 0? Is it 
still algebraically integrable? 

10.6 Chaplygin 

The themes of Section Uni do not occur in Chaplygin jH]. Because theta functions are defined 
in terms of complex coordinates, one might argue that the sentence "From (29) we see 
that u and v can be expressed in terms of theta functions of the two arguments a and 
r" in Chaplygin P[ after (30)] yields implicit evidence that Chaplygin did think of complex 
variables, as does the sentence "Solving equation (41) gives two real values for the quantity /" 
in Chaplygin [HI after (41)]. However, the inequalities between (27) and (28) in Chaplygin j^l 
§3] indicate that Chaplygin mainly focussed on the real system, whereas he also emphasizes 
that (41) has real solutions. Completion of the complexified system and tori (real or complex) 
definitely do not occur at all in Chaplygin 0. 

11 Hyperelliptic Integrals 

In this section we assume that the moments of inertia /j are different from each other, that 
p 7^ 0, that the constants of motion (j, T) are at a nonsingular level. We also assume in this 
section that the the moment j of the momentum around the point of contact is nonzero and 
horizontal, which means that 7^ and J3 = in ()3.9|) . As shown in Subsection 19.21 the 
rotational motion with arbitrary nonvertical j can be reduced to this case. 
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In order to obtain a smooth level surface, we will assume that the kinetic energy T is not 
equal to any of the critical levels 

TcriM:=||jir/2(/^ + p), ^ = 1, 2, 3, (11.1) 

of the function Tj on SO (3), cf. ()4.9j) with = 0. In order to obtain that the complex level 
surface is smooth, we will need furthermore that 

2Tp^\\jf. (11.2) 

If M has real points, then ()11.2p is a consequence of the assumption that T ^ Tcrit,i, because 
T < Terit,i = ||j||V2(/i+p) < ||j||V2p, cf. (UnD, and therefore 2T p < In the 

case that we allow arbitrary complex values for the parameters /j, p and we have to 

add ()11.2|) to the list of conditions. In other words, we make the same assumptions as in 
Subsection 110. II 

11.1 The Projection onto the First Vector 

Our next goal is to simplify the vector field ^ by means of a suitable substitution of variables 
in the w-space. We recall that u has the concrete interpretation that —ru is equal to the 
point of contact on the surface of the sphere in body coordinates, cf. fj2.7|) . 

We have the two complex surfaces L and Uc in which L is the surface in x C"^ 
determined by the equations ()3.9p and ()3.14p and f/c is the quadric 

Uc := {ueC^l (n, u) = 1} (11.3) 

in C^. The projection {u, f ) t— > u is a branched covering from L onto Uc, branching over 
the set of M e Uc for which there exists a solution t> G C'^ of {u, v) = 0, {v, v) = 
and ()3.14p . where the derivatives of {u, v), {v, v), and f{u, v) with respect to v are linearly 
dependent. If X{u) = 0, then f is a solution of ()3.14|) if and only if Y{u, v) = 0, in which 

case 

df{u, v)/dv = 2Y{u, v) dY{u, v)/dv - X{u) dZ{v)/dv = 0. 

Therefore the zeroset of X, which does not contain any real points in view of the remark 
preceding Lemma f7. 11 is contained in the branch locus. 

When u G Uc and X{u) ^ 0, then m is a branch point if and only if there exists a solution 

V of ()3.9|) and ()3.14|) such that the vectors u, v and 

w:= X{u)uj = Y{u, v) Ju + X{u) Jv (11.4) 

are linearly dependent. Here we have used the formula ()3.1H) for u. Assuming that j is 
not vertical, it follows from ()3.8|) that u and v are linearly independent, and we obtain 
that w = au + Pv for suitable constants a and p. With the abbreviations X = X{u), 

Y = Y {u, v) , we have 

a = {u, w) = Y {u, J u) + X {u, J v) = Y j p. 
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cf. fj!rTH|l and (nrTT)|) . and 
cf. fl!nH|l . dSm, and (HrTTj) 



With the notation 



(t;, w) = Y {v, Ju) +X {v, Jv) = 2TX, 



T := 



2T/iijir 



(11.5) 



the equation w = au + [3v takes the form X {J — t) v = Y {1 j p — J) u. Because = 
I + p, cf. (I2.18|l . we have that (1 — p J) = / and it follows that v is equal to a nonzero 
multiple of Ku, in which 



In other words, we have proved that away from X{u) = the branch locus is contained in 
the intersection of the quadric Uc with the quadratic cone defined by ()11.7|) . Because this 
intersection is irreducible, the conclusion is that the branch locus away from X{u) = is 
equal to the intersection of Uc with the quadratic cone defined by (|11.7|) . 

Remark 11.1 The matrix K in ()11.6|) is a diagonal matrix with eigenvalues equal to 



in which the Tcrit,i are the critical levels of Tj, cf. flll.5|) and pLlj) . The assumption that we 
are on a regular level set is equivalent to the condition that T 7^ Tcrit.i for every i = 1, 2, 3. 
If the level set contains real points, then we have that Tcrit,3 < T < Tcrit,2 < Tait,! or 
^crit,3 < Tcnt,2 < T < TcHt, 1- In the first case K has one negative and two positive eigenvalues 
and in the second case K has two negative and one positive eigenvalue. Therefore, in both 
cases the cone defined by pi.7|) has nonzero real points. The real cone intersects the real 
unit sphere U in two closed curves, diametrically opposite to each other. 

Because each of the two connected components of the real level surface of (j, T) is a 
torus, it can only be mapped onto the annular region between these two curves. (This can 
also be verified by means of explicit calculations, cf. the last part of Remark 111.21 ) The 
images of the quasiperiodic solution curves on the level set, cf. Corollarv 18.41 while running 
around the sphere in this band, run from one of the bounding curves to the other, every 
time with a contact of order two at the bounding curve. A similar behaviour occurs when 
j is neither horizontal nor vertical, but in that case the bounding curves are determined by 
more complicated equations. 

Remark 11.2 If = then, for given u ^ Uc, the solutions v of the equations (j3.9p and 
()3.14|) can be explicitly computed in the following way. Using the equations {v, v) = 
and ()11.5|) . we can write the kinetic energy equation ()3.14|) in the form 




(11.6) 



(u, Ku) = 0. 



(11.7) 



= hi (1 - r (/, + p)) = hi (1 - T/Tcrit,.) , 



(11.8) 



X{u) {v, {J -t) v)+Y{u, vf = Q, 



(11.9) 
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which is homogenous (of degree two) in the variable v. The condition that j is horizontal 
corresponds to the equation 

(m, v) = 0. (11.10) 

In solving the homogenous equations ()11.9|) and ()11.10|) for v with f 3 7^ we may put = 1. 
If U2 7^ 0, then we can solve V2 in terms of vi from ()11.10|1 . This leads to a quadratic equation 
for Vi of the form 

a-siu) fi^ + 2b2{u) vi + a2{u) = 0, 



in which 



a3[u) 
h2{u) 
aiiu) 



X{u) (Ci U2^ + C2 Ui^) + (Ci - C2f Ui^ M2^ (H-H) 

[X{u) C2 + (ci - C2) (c3 - C2) ^2^] Ml M3 and (11.12) 

X{u) (C2 Ms' + C3 M2') + (C2 - C3)' M2' (11.13) 



Here we have used the abbreviation 

Q:=l/(/, + p)-r. (11.14) 
With these notations, we obtain, using repeatedly that (m, u) = 1, that 

vi = {-h2{u) ± A(n)i/2) /a3(n), (11.15) 
in which the discriminant turns out to be given by 

/\2{u) \= b2{uY - a^iu) aiiu) = -CiC2CsU2^ X{u) {u, K u)/ p, (11.16) 
with K as in (1111)1) . 

The solution vector v is determined by taking V2 ■= — (mi f 1 + ^3) /u2- In order to obtain 
a vector of length we have to replace v by cv, in which c = ± ||. We obtain four 

solutions, consisting of two oppositie pairs. 

As expected, the discriminant is equal to a multiple oi X[u) {u, K u). Note that for real u 
we have real solutions v if and only if ci C2 C3 (u, m) < 0. It follows from ()11.8j) . pi.5|l and 
()11.14j) that Ci Bi = li/ {li + p) > 0, and therefore the determinant of ci C2 C3 -ft^ is positive, 
which implies in view of Remark 1 1 1 . 1 1 1 hat ci C2 C3 has two negative eigenvalues. It follows 
that the part on the unit sphere where ci C2 C3 {u, K u) < is the connected annular region, 
bounded by the two smooth curves determined by the equations {u, u) = 1, {u, K u) = 0. 

For u in the domain Ci C2 C3 {u, Ku) < 0, the interior of the annulus , we obtain four 
solutions V, which correspond to four possibilities for the velocity vector One pair of 
these velocity vectors correspond to one of the connected components of the level set, and 
their opposites correspond to the other connected component, cf. the discussion at the end 
of Section EJ These connected components are mapped to each other by means of a rotation 
over TT which maps j to its opposite, cf . the discussion after Lemma 14.11 

The intersection of the quadric {u, u) = 1 with the cone {u, Ku) = is equal to the 
intersection of the quadric {u, u) = 1 with the quadric {u, {1 + P K) u) = 1, in which /? is 
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any nonzero constant. If we choose /3 = r, then we obtain that the branch locus is equal to 
the intersection with Uc of the quadric defined by the equation 

{u, {1 - J ly^ u) = I, in which 7 := r/(l - rp). (11-17) 

The equation X{u) = is equivalent to {u, p {1+ p)^^u) = 1, cf. ()3.15p and ()2.18|) . The three 
quadrics {u, u) = 1, X{u) = and ()11.17p therefore belong to the one-parameter family of 
quadrics {u, (1 — XI)^^u) = 1, in which the parameter A takes the values A = 0, A = — 1/p 
and A = 7, respectively. The substitution of variables u = J^/^ x leads to 

(u, (1 -Xiy^u) = {x, /(I - A/)"^a;) = {x, {r^ - X)'^ x), 

and our family of quadrics is turned into the one-parameter family (pencil) of confocal 
quadrics 

E^ = l, (11-18) 

in which 

Qi := hence Ui = Xijai^^'^ . (11.19) 

11.2 Jacobi's Elliptic Coordinates 

Given x G C'^, the equation (jll.lSp corresponds to a polynomial equation of degree three 
for A, which has three solutions (Ai, A2, A3), which are called Jacobi's elliptic coordinates, 
cf. Jacobi [211 26. Vorlesung]. We briefly recall some of Jacobi's observations. 

Let Qi, 1 < i < n, he n different numbers, and let s G Z>o. Then the function z 1— >■ 
z'^ / Ylk ~ (^k) has simple poles at the points z = a^, with residue equal to a// H/cifc^j i^i ~ '^n 
It follows that 

z^ ^ a/ 1 

^ ^ Uk (^ - ^ Uk\k^^ («i -am) Z- ai 

is equal to a polynomial of degree s — n when s > n and equal to zero when s < n. Therefore, 
if s < n we have that 

z'^ aj^ 1 

E 



rifc - ak) Uk\k^i («i -a„,) z-ai' 

If we compare the expansions in powers oi 1/z for 2; — > cxd in the left and right hand side, 
we obtain that 

V = f ^ ^^^^"^ s = n-l, 

tr Uk\k^^ («« - «-) 10 *hen < s < n - 1, 

If Aj, 1 < 2 < ra, are arbitrary numbers and we write 

x.' = y.= ^' (a. -A.) ^^^21) 
\[k\k^i {a^-ak) 

78 



then an expansion of the numerators in the right hand side of 



in powers of Oj yields in combination with pi.2()j) that 



1=1 



E 



A, 



l<l<n. (11.22) 



In other words, if the Ui are defined by pi.2Hl . then each of the Aj's is a solution of 
the equation yi/ (a^ — A) = 1 for A. Note that ()11.21|1 defines a polynomial map- 
ping (Ai, ■ ■ ■ , A„) ^— > ■ ■ ■ , Un) from C" to C". It is a branched covering because for 
every permutation vr of {1, n} the vector {Xt,{i), A^(„)) has the same image as 
(Ai, A„). If the Xi satisfy ()11.21|) . then each Aj is a solution of the equation ()11.18|) . 
However, ()11.21|) does not define a mapping (Ai, ■ ■ ■ , A„) t— > (a;i, ■ ■ ■ , x„) , because the map- 
ping (a;i, ■ ■ ■ , Xn) ^ {yi, ■ ■ ■ , yn) is a 2"-fold branched covering — for any choice of signs ej 
the vector (ei xi, . . . , e„ Xn) has the same image as (xi, . . . , a;„). The branching occurs at 
the set where one of the coordinates Xi is equal to zero, which corresponds to the condition 
that one of the coordinates of A is equal to Oj. 

In our case we have n = 3 and ()11.19|) . The sign changes in the Xi correspond to sign 
changes in the Ui. For the equations of motion in the elliptic coordinates these will not cause 
too much trouble, because if i? is a diagonal matrix with ±l's on the diagonal, then the 
transformation (m, v) i— > [Ru, Rv) will leave the vector field invariant when deti? = 1 and 
turns the vector field into its opposite (corresponding with a time reversal) when det R = —1. 

We will keep Ai = 0, which implies that the equation (u, m) = 1 is fulfilled, and regard the 
remaining two elliptic coordinates (A2, A3) as coordinates on Uc- Apart from the problem 
that the velocity field for u is multi-valued, it will have singularities at all the branch loci, 
corresponding to the condition that A2 or A3 attains any of the five values — 1/p, 7, I//1, 
I//2 or I//3. 

The relation ()11.2H) induces a relation between tangent vectors X and A in the x-space 
and the A-space, respectively. If we take the logarithm of the left and right hand side and 
differentiate,then we obtain that 

2^ = 5^-^. (11.23) 

k 

Squaring this and inserting the formula 1)11.211) for Xi^, we obtain that 

* - 2^ — ff 7 ^ — ^^k^h + ^Yi '( V U-'--24j 
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The sum over i of the first sum in the right hand side, over the k and / with k ^ I, vanishes 
in view of (jl 1.201) with s < n — 1. On the other hand 



k) 



because as a function of the left and the right hand side both vanish at infinity and have 
the same poles and residues. This leads to Jacobi's conclusion, cf. [211 26. Vorlesung], that 

" A 2 

^E^'' = E n (^'-^'^O ^7^, in which (11.25) 

S{\k) := n (11.26) 
I 

From now on we use that Ai = 0, Ai = 0, and = 3. Then the quotient of the first sum 
in the right hand side of pi.24|) by vanishes in view of ()11.20j) with s < n — 3, where we 
have used that no terms appear with A; = 1 or Z = 1. On the other hand, for every /c 7^ 1 we 
have that 



3 



rimlm^i («i - «™) a* - ^k ]\i {ai - Xk 



^—1 i i.m\m^i 

because as a function of A^ the left and the right hand side both vanish at infinity and have 
the same poles and residues. With the notation ()11.26|) we therefore obtain in our case n = 3 
that 

4|:_.„,-A3)(^-^). (11.27, 

For the next equation we will use the index z in a cyclic manner, by taking i G Z/3Z. 
Then (111.2311 implies that 

diXi aj+iXj_|_iA / tti aj_|_i A (ttj — Oj+i) AfcAfc 



Xi Xi+i J ^ \ai — Xk Oj+i — Afc 7 (ttj — Afc) (aj+i — Afe) 



On the other hand it follows from ()11.21|) . Ai = and n = 3 that 

Xi'^Xi+i'^ {ai - A2) {tti - A3) (oi+i - A2) (flj+i - A3) 



didi+i {ai — aj+i) {ui — aj+2) (oi+i — O'i+2) {di+i — ctj) 

It follows that 

a^j^Xj-i-i^ / cii Xi ttij^i Xij^i 



\ Xi Xi^i 



{tti - A2) {tti - A3) {tti+i - A2) {tti+i - A3) Xk Afc 



(flj+2 — flj) (flj+1 — 0,1+2) \h ~ {Oi+1 — Xk) 
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which is equal to 1/ (ai+2 — (oi+i — ai+2) times 

(Oi - A3) (Oj+i - A3) 2 A 2 I o\ A \ A I ('^i ~ -^2) (Oj+l - A2) 2 a 2 
-— rT^2 A2 + ZA2 A2 A3 A3 + — — TT'^S A3 . 

(Oj - A2j (Oi+i - A2) [ai - A3) (Oi+i - A3) 

Now we have 

1 tti — Oj+i T-r , , 

(ai+2 - Ci) (Ci+i - ai+2) P .^^^^ 

{ai - A3) (fli+i - A3) _ (A3) ai+2 - A2 _ 5" (A3) ( A3 - A2 



(flj — A2) (fli+i — A2) 5* (A2) aj+2 — A3 S (A2) V ai+2 "~ A3 
the cychc sum over i oi a^ — Oj+i is equal to zero, and finally 



V 



v-^ _i_ ^ — (flj - Oj+i) (oj - A3) (oj+i - A3) — 



It follows that 



igz/^z '^^"^+^ ^ ^^+^ ^ V'5(A2) 5 (A3) ^ 



11.3 The Motion in Elliptic Coordinates 

In order to obtain the time derivative A = A of A corresponding to the velocity vector X = i, 
which in turn corresponds to the time derivative u := m x a; of m, we start with the observation 
that the horizontality of j implies that 

= (f, u) = {Ip,ui^, u) = {lu, u) = {u, lu), (11.29) 

cf. firmUl . dSHI), and (jTTH|l . It follows that 

ii X I u = {u X uj) X I u = {u, lu) cu, (11.30) 

which equation will be used in order to express uj in terms of u. 
From (|3.7p we obtain that 

2T = {loj.u) + p{u,u). (11.31) 

Furthermore, 

{I {il X I u) , il X I u) = ^ Ii {Ui+i Ii+2 Ui+2 - Ui+2 h+l Ui+if 
= ^ {li-l h+l^ Ui+i^ + /i_2 h-l^ Mj-l^) uj - 2/1 I2 h ili+l Ui+2 Ui+2 

«ez/3Z 

iez/3Z ^ * * ^ 
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where in the third equation we have used that 

= (it, = ^ ^2 ^.2 _^ 2Ui+i Ui+i Ui+2 Ui+2- 

iGZ/3Z 

It follows that the first term in the equation for the kinetic energy can be written in the form 

{lu, to) = hhh {I'^u, u)/{u, lu). (11.32) 
With the substitutions ()11.19p . we have 

3 

{I-^u, u) = Y,^l (11-33) 



i=l 



and 



(u, u 



) = (Jx, i;) = ^ x1/ai 



11.34) 



1=1 



which by means of (jll.25|) and (|11.27p . respectively, can be expressed in terms of the velocities 
of the elliptic coordinates. 

In order to express the denominator 



(u, I u) = {x, P x) = Xi'^/c 



;il.35) 



i=l 



in ()11.32p in terms of the elliptic coordinates, we observe that 

3 

1 = {U, U) = ^ Xi^/tti 
i=l 



11.36) 



implies that the third degree polynomial equation, which is obtained from (jll.l8p by mul- 
tiplication with S{X) = nLi i^i ~ '^)' A = Ai = as a solution. The second order 
equation for the two remaining solutions A2, A3 takes the form 



^ (xi^ - ai) 



i=l 



A - ^ Xi^ (aj+1 + a.;+2) + ^ (^h cih+i = 0, 

ieZ/3Z /16Z/3Z 



in which the constant term can be simplified to 



ieZ/3Z 

It follows that 



-di+l — (11+2 



ai 



E 



2 Ctj+i 0.4+2 \ ^ Xi 

Xi = fli 02 03 > — . 

Oj -^-^ Oj^ 



A2 + A3 = ^ 



11.37) 



i=l 
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and 

A2 A3 = ai 02 as (11.38) 

^ or 

1=1 

Combining (im^ . (im^l (ITT:^ for X = x, and we conclude that 



3 2 



/1/r \ A2 — A3 / A2A2 A3A3 \ 
Combining (111.311) and (111.391) with (111.3411 and (lll.27|) for X = i, we arrive at 

- ^ ^ - IS) - ^ - ( A - A)) • '-^«> 

In order to obtain the second equation for the two unknowns A2, A3, we start with the 
equation 

= {I u + p u — p {u , u) u, I u — p {u , u) u) = {I UJ, I u) + p {u, I Uj). 

Here we have used (13. 9|) and ()3.7p in the first equation, ()2.16|) in the second equation, and 
(jll.29p in the third one. The first term in the right hand side is equal to 

{I oj, luj) = {u, lu)^"^ ■ {I {u X lu), I {ii X lu)) 



jez/3Z 

22/ • ■ \ ^ 

{u,Iu) (010203) V 



where we have used (jll.3(jp in the first equation and (jll.l9p in the second one. Combining 
this with (ril.35|) . (111.381) . and (111.281) for X = x, we obtain that 



A/T T \ A2 — A3 /A2^A2 A3^A|\ 



(11.41) 



Also using (jll.4(Jj) . we therefore arrive at 



Ag — A3 / A2^ A2 As^ A| \ A2 — A3 / A2 A2 A3 A| 



For the unknowns 
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the equations ()11.4U|) and ()11.42|) take the form 

8T = ^2-^3 and 4||jf - 8pT = ^ 6 - ^ ^3, 
respectively. Substracting A2 times the second equation from the first one, we obtain that 



where we have used ()11.17|) in the second equation. 

At this stage we recall the change of the time parametrization defined by 

cf. Corollarv 18.41 and (|3.15|) . In order to express X{u) in terms of the elliptic coordinates, 
we recall that 



3 2 

pX(n) = l-^^^ with X= -1/p. 
i=i ~ ^ 



Now 

3 2 



^(^) - E ^ j = - A) (As - A) , (11.44) 

because as a polynomial in A the left and the right hand side have the same zeros Ai = 0, 
A2 and A3 and have the same leading coefficient. It follows that 

3 

X{u) = p-' (A2 + 1/p) (A3 + 1/p) / n + Vp) • (11-45) 



i=l 



Substituting A| = (^) /X{u) and (|11.45;i in (jll.4;-i|l . we arrive at 



('-^V/P(A.)=(^V. (11.46) 



V dr / \ A2 — A 

in which the polynomial P is defined by 

3 

P(A) := (-1/p - A) (7 - A) n - A) (11.47) 
and the constant c is determined by 

c':=(4||j|p-8pT)/n(pa. + l). (11.48) 

i=l 
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In a similar way we obtain the equation 



for dAs/ dr. 

At this point it becomes appropriate to introduce the hyperelliptic curve C defined by 
the polynomial P as the one point completion (compactification) of the affine curve 

{{X, z) eC^ \ = P{X)}. (11.50) 

The genus of the hyperelliptic curve C is equal to g if the degree of the polynomial P is 
equal to 2(yf + 1 or 2g + 2, cf. Shafarevich i35, Ch. Ill, §5], or Farkas and Kra 1.13^ III. 7. 4]. 
Because in our case the degree of P is equal to five, we have g = 2. 
Consider the motion on C x C which is determined by the choice 

1^^ = -^ and 11h = ^^ (11.51) 

Z2 dr A2 — A3 Z3 dr A2 — A3 

of the square roots in ()11.46|) and ()11.49p . respectively, where (A2, Z2) and (A3, Z3) denote 
the points in C which, under the two-fold branched covering C 3 (A, z) 1-^ A, project to A2 
and A3, respectively. 

The equations 1)11.511) imply that 

A2 dA2 _ A3 dAs _ Q J_ dA2 L 1:^ — c (1152) 

Z2 dr Z3 dr ' Z2 dr z^ dr ' 

or 

/•(A2(t),Z2(t)) ^ /•{X2{t),Z2{t)) 

/ — dA = constant, / z^^ dA = cr + constant. (11.53) 

J{Xa{T),Z-,{T)) Z J{X3{T),Z:i(T)) 

Here the integration is over a curve in C, running from (A3(r), 2:3 (r)) to (A2(r), 2^2 (t)), 
and depending smoothly on r. The primitives of the differential forms z~^ dA and ^ dA 
are the hyperelliptic integrals of the hyperelliptic curve C, and for this reason one says that 
flll.53j) implies that the problem is solved by quadratures in terms of the hyperelliptic integrals 
corresponding to the hyperelliptic curve C. The quantities A2 — A3 and A2 A3 are given by 
Jacobi's theta functions of the integrals in the left hand sides of ()11.53j) . cf. Shafarevich (HHl 
p. 419], and therefore one also talks about solving by quadratures in terms of theta functions. 

11 A The Jacobi Variety of the HypereUiptic Curve 

The equations in (|11.53|) lead to a beautiful interpretation of the system in terms of the 
Jacobi variety Jac(C) of the hyperelliptic curve C. The differential forms 

■= _ dA, i = 0, 1, (11.54) 
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extend to holomorphic differential forms on C and actually form a basis of the two-dimensional 
complex vector space (C) of all holomorphic differential forms of degree one on C, cf . Sha- 
farevich ^[ Ch. Ill, §5] or Farkas and Kra |13| III. 7. 5]. 

Remark 11.3 It is clear that the differential forms 

(3' ■= ^ dAs - — dAs, 1 = 0,1 (11.55) 

of degree one on C x C are holomorphic, and closed because of the separation of variables. 
In this notation, the equations in ()11.52j) mean that the vector field 

on C X C satisfies 

k/3[ = 0, k/3', = c, (11.57) 

if i^ (3 denotes the inner product of the differential form (3 with the vector field ^. This is a 
complex version of Proposition 18.31 where we had the differential forms j3 and 7 instead of 
P'q, and (3 = and ig 7 = —1. 

For each smooth curve 7 in C, we have the complex linear form 

n\C) [ p (11.58) 



on 1H}{C) — this defines an element of the dual space 1H}{C)* of T-C^{C). The integral 
depends only on the homology class [7] G IIi(C, Z) of 7. If we restrict ourselves to closed 
loops 7, then this leads to a homomorphism J from the Hi(C, Z) to 1H}{C)* . The image 

A(C) := j (Hi(C, Z)) c H\Cy (11.59) 

is an additive subgroup of Ti}{C)* , which is called the period lattice of C . For any Riemann 
surface (complete algebraic curves over C) of genus g we have that the complex dimension of 
TH}{C) is equal to g and therefore the real dimension of TH}{C)* is equal to 2g. Furthermore, 
the period lattice A(C) has a Z-basis which at the same time is an R-basis of 7Y^(C)*, cf. 
Farkas and Kra [1^ III. 2. 8]. Therefore the quotient space 

Jac(C) := 'H\Cy/k{C), (11.60) 

is compact, a torus of real dimension equal to 2g. Definition ()11.60p is the analytic definition 
of the Jacobi variety of the curve C, cf. ^2 p. 87] or jSH p. 143], whereas the algebraic 
definition is formulated in terms of divisors, cf. [221 P- 155] or cite[p. 3.28]TthII. 
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If we fix p, g G C, then the difference of the homology classes of two curves in C which 
run from p to g is an element of Hi(M, Z), and it follows that, for each (3 G T-C^{C), the 
element ^ 

f :=(/")+ A{C) G Jac(C) (11.61) 

J p J 'y 

does not depend on the choice of the curve 7 from p to q. This defines a smooth mapping 
/ : (P, q) ^ Ip from C X C to Jac(C). Note that = - for every (p, q) E C x C, 
and =0. If L denotes the involution (A, z) 1-^ (A, —z) in C, then l*P = —(3 for every 
(3 G Ti}{C), and it follows that / maps (p, q) and to the same element of Jac(C). 

Therefore the mapping j can also be viewed as a mapping to Jac(C) from the quotient 
C x^/ C of C X C by the involution l' : (p, q) {/-{l), '•(p)) in C x C. The diagonal 
{{p, q) E C X C \ p = q} in C X C , which is isomorphic to C, is mapped by the quotient 
map to a curve D in C x^/ C which is isomorphic to C/l ~ Pi(C). Because the diagonal is 
mapped to the origin, we have that / (D) = {0} as well. 

The mapping A is closely related to the Abel-Jacobi map A : C x C ^ Jac(C), which is 
defined as follows. Choose a base point po G C. Then A{p, q) is equal to the linear form on 
li}{C) modulo A(C) which is defined by 

rv PQ 

Because A{p, q) = A{q, p), A can be viewed as a mapping to Jac(C) from the symmetric 
power C^^) of C, which is defined as the quotient of C x C by the involution (p, q) 1— >• (g, p). 
It is a classical theorem of Abel and Jacobi that for genus two curves C the Abel-Jacobi map 
is surjective, cf. Farkas and Kra JSl III. 6. 6]. More precisely, it maps a genus zero curve Dq 
m to a point jo in Jac(C) and is a diffeomorphism from C*-^^ \ Dq onto Jac(C) \ {jo}, 
cf. Farkas and Kra 13, III. 11.8 and III. 11. 11]. This is an example of a blowing down, also 
called a sigma-process, cf. Shafarevich [23 Ch. II, §4 and Ch. IV, §3]. Now, if l (po) = po, 
then 

pq pp pq pi(p) pq pl{p) pq 

A{p, <!)= + = - =/-/ =/■ 

Jpo Jpo Jpo J (^{po) Jpo Jpo J i-ip) 

The condition that l (po) = Po means that po corresponds to one of the five zeros of P or 
to the point on C at infinity. It follows that the mapping J also maps a genus zero curve 
to a point and is a diffeomorphism from the complement of the curve to the complement in 
Jac(C) of the point. Because J{D) = {0}, the curve must be equal to D and / defines a 
diffeomorphism from C x,,/ C \ D onto Jac(C) \ {0}. 

In these terms the equations in ()11.53|) imply that the tangent map of / maps the vector 
field ^ in C X C to a constant vector field on Jac(C). More precisely, if we identify 1H}{C) 
with by means of the basis /3o, [3i of ()11.54|) . then pi.53p yields that ^ corresponds to the 
constant vector field with coordinates (c, 0). We therefore have verified that the rotational 
motion of Chaplygin's sphere with horizontal moment is algebraically integrable according 
to the definition of Adler and van Moerbeke 0, and in view of Subsection 19.21 this result 
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remains true if we only assume that the moment is not vertical. This verification is very 
different from the one via proposition I1U.81 

Because every holomorphic vector field on Jac(C) lifts to a bounded holomorphic vector 
field on TC^{C)*, and every bounded holomorphic function on is equal to a constant, we 
have that every holomorphic vector field on Jac(C) is constant. Therefore our conclusion 
that C, is mapped to a constant vector field on Jac(C) is equivalent to the conclusion that C, 
is mapped to a holomorphic vector field on Jac(C). According to (jll.51|) . the vector field ^ 
on C X C is rational with poles along the diagonal. The blowing down of the diagonal by 
the map / to a point (the origin) in Jac(C) apparently has the effect that it regularizes the 
vector field ^. 

Remark 11.4 If the rotational motion after time parametrization as in Corollarv 18.41 is 
truly quasi-periodic in the sense that the orbit on the (j, T)-level set is dense, then each 
continuous vector field which commutes with ^ is equal to a^ + for^, in which a = a{j, T) and 
b = b{j, T) are constants which may depend on j and T. Because the levels (j, T) , for which 
the ^-orbits are dense in the level set, form a dense subset of the set of all {j, T), it follows 
that if a continuous vector field is defined on all the regular level surfaces and also depends 
continuously on [j, T), then it is of the form + brj with a = a{j, T) and b = b{j, T) 
depending continuously on (j, T). If 60 and bi are arbitrary constants, then the vector field 
C, which in the elliptic coordinates corresponds to 



compare ()11.57|) . or on Jac(C) corresponds to the constant vector field with coordinates 
{bo, bi), commutes with ^ and depends continuously on {j, T). The conclusion is that there 
is a bijective linear correspondence between the vector fields which commute with ^ and have 
the aforementioned continuity properties, and the constant vector fields on Jac(C). 

In order to determine the coefficients 60 and 61 for which ( = t], where t] is as in Propo- 
sition IH31 we start with the time derivative x which corresponds to 



where in the second identity we have used ()3.1H) . ()2.16j) . and u x u = 0. It follows that 



Here we have used ()3.9|1 in the first equation, {u, v) = — js = in the second one and 
Ui = li^^^ Xi = aC^I'^ Xi in the last one, cf. In view of ()11.27j) . we obtain the relation 



(11.62) 



it := u X [I + p) uj = u X V 




3 




(11.63) 



for the corresponding time derivative A = A in elliptic coordinates. We have that 



dr ~ 



X(m)^/2 with X{u) given by (|11.45jl . satisfies (|n.62jl if and only if 
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or, equivalently, 

dA2 1^ _ bi - A3 bp dA3^^^ _ bi - A2 60 
dr A2 — A3 ' dr A2 — A3 

Squaring each of the expressions and using (jll.SOp . ()11.47|) . we can bring the equation (fTT 
into the form 

4||jf _ (A2 + I/P) (A2-7)(&i-A3&o)'-(A3 + l/p) (A3-7)(&i-A2&o)^ 
P' U^ (a. + 1/p) (A2 + 1/p) (A2 + 1/p) (A2 - A3) 

The numerator in the right hand side has to be equal to zero when A2 = which imphes 

that 

bi = -bo/p. (11.64) 
Substitution of (jll.64j) in the previous formula yields that 

4||jr _ ^^2 (A2 - 7) (A2 + 1/p) - (A2 - 7) (As + 1/p) ^ ^^2 ^ ^) 



p' n-=i + 1/p) - A3 

or 

bo' = A\\jr/p' (1/p + 7) n K + I/P), (11-65) 

1=1 

a formula similar to (jll.48p . The formulas (jll.65|) and (|11.64j) determine the constant vector 
field (60, bi) on Jac(C) corresponding to r], up to its sign. 

Until now we have not paid much attention to the fact that the substitutions which we 
have used are not bijective but, except for the mapping J : C x^/ C ^ Jac(C), are branched 
coverings. Recall that in order to obtain a single-valued vector field ^, we passed from the 
manifold L of solutions {u, v) of the equations ()3.9|) and ()3.14|) to its two-fold covering M, 
defined by ^HM . 

The projection [u, v) ^ u exhibits L as a fourfold branched covering over the quadric 

Uc := {u G I (n, u) = 1} 

in C'^, with branch locus at X{u) {u, Ku) =0, cf. Subsection II 1.11 and (jll.Tj) . At the generic 
points of Uc we have four possibilities for of the form ±,^', ±^", where ^' and ^" become 
equal at the branch locus. 

The branched covering L —>■ Uc is less than fourfold over the set D^o of the limit points 
of the u E Uc for which at least one of the solutions v of the equations ()3.9|) . ()3.14|) runs 
off to infinity. This is the set of m G Uc for which there exists a nonzero solution v of the 
homogeneous equations 

{u,v) = 0, {v,v) = 0, {u, Jv)^ + X{u) {v, Jv) =0. (11.66) 

This exhibits D^o as the resultant set of the three polynomials in v which appear in pi.66|) . 
and therefore Doo is a closed algebraic curve in Uc- 
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The set Dqo contains the set D of the points which correspond to points on the diagonal 
A2 = A3 in the elhptic coordinates, because according to (|11.51|) the vector field in the 
(^^2, A3)-space is infinite along the diagonal. At the points of D both vectors v and or 
w become infinite, which means that all solutions v of the equations ()3.9|) . ()3.14p run off 
to infinity when u approaches a point in D. This implies that the image of L under the 
projection (u, t>) 1— > u is contained in Uc \ D. According to pi.37p and ()11.38p . the point 
u E Uc corresponds to A2 = A3 in the elliptic coordinates, if and only if 

A{u) := {{r^u, u) - tmce - 4:{I u, u) / det I = 0. (11.67) 

Lemma 11.5 We have that = D, which implies that the projection L Uc is a fourfold 
branched covering from L onto Uc \ D- 

Proof (Sketch.) For (m, u) = 1 and Ui^ + 7^ the nonzero multiples of the f G 
such that 

= -Ml -^3 ± ^2 = -^3^2 =F = U2 + Ui^ 

parametrize the set of nonzero solutions v of the equations (m, f ) = and (f , v) = 0. 
Inserting this v in the last equation in (|11.66|) . where we use that J = (/ = p)~\ cf. (j2.18|) . 
we obtain the polynomial equation 

[{h - h) h Ul^ + {h - h) h U2^] [Ul^ + U2^] + [h - h) h {U^ - U2^) 

= 2 ± {h - h) h ui U2 Us. 
Squaring both sides we obtain a polynomial equation of the form 

Aiu) {m^ + u2'Y = 0, 

where we have used repeatedly that {u, u) = 1. A straightforward calculation shows that 
does not contain the points u such that {u, u) = 1 and Ui"^ + M2^ = 0, and therefore the 
equation for Doo is equivalent to the equation ()11.66|) for D. □ 

The mapping u ^ y, yi = Xi^ = ai Ui^ exhibits Uc as an eightfold branched covering of 
the plane 

P:= {yGC=^|X:- = l}, 

with branch locus at the three coordinate planes, at yi y2 y^ = 0. We arrive at a 64-fold 
branched covering M ^ P \ D, where we denoted the image of D C Uc in P under the 
mapping from Uc to P with the same symbol D. 

In order to reduce the order of the covering, we will use the group S introduced in ()1().31|1 . 
The group S leaves the fibers of the 64- fold branched covering M P \ D invariant, and 
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therefore the latter is equal to the composition of the projection M M/S and a uniquely 
determined mapping 

TT,,/, : M/S ^ P \ A (11.68) 
which is a four-fold branched covering. 

In the other direction we have the mapping (A2, A3) 1— > y defined by Jacobi's elliptic 
coordinates, this is a twofold branched covering with branch locus equal to the image of the 
diagonal A2 = A3, cf. Subsection 111.21 If we replace the (A2, A3)-space by the symmetric 
power C*^^-*, then we actually obtain a birational isomorphism with the inverse of ()11.2H) 
given by (jll.37|) . (|11.38|) . (Don Zagier pointed this out to me.) 

Let Caff denote the affine (finite) part of the hyperelliptic curve C . Then the projection 
(A, z) H-i> A exhibits Caff as a twofold branched covering over C. This leads to the following 
commuting diagram of branched coverings: 

Caff X Caff Caff X Caff 

i 4:1 [ 4:1 

C X C ^ C(2) 

Finally we have the mapping j : C x^i C —>■ Jac(C) which blows down the diagonal D to 
the origin in Jac(C) and defines an isomorphism from (C x^/ C) \D onto Jac(C) \ {0}. The 
image of Caff x^/ Caff in Jac(C) is equal to the complement in Jac(C) of the image of C under 
the mapping '■ Q ^ from C to Jac(C). (On the real axis, we have that 00 = —00.) 
This mapping is an embedding of C into Jac(C), cf. Farkas and Kra TT, III. 6. 4]. We denote 
the image of C in Jac(C) under this mapping by Note that /~ = 0, hence G It 
follows that J defines an isomorphism from (Caff x^/ Caff) \ D onto Jac(C) \ J^. 

The inverse of /, followed by the projection from (Caff x^/ Caff) \ -D to C'-^-* ~ P, defines 
a fourfold branched covering 

7rjae(c) : Jac(C) \ [ ^P\D. (11.69) 

In combination with the fourfold branched covering pi.68|) . this leads to the following 
identification of M/S with Jac(C) \ J^. 



Proposition 11.6 Let ^ be the velocity vector field on M/E and denote by ^•^^'^('^^ the con- 
stant vector field on Jac(C) defined by fjll.52j) . Then there is a unique isomorphism ip from 
M/Tj onto Jac(C) \ J^, such that tTj^^^^^ o ip = tt^^^^^, and 



{(.m) — T^(m) ^jac{C) (^^ip{m) ^) 

for every m G M/S. In other words, the diagram 

M/S ^ Jac(C) \ X 

P\D 



c 
00 



M/E ^1 ^ " ,Iac(C) 
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commutes and ip intertwines the vector field ^ on M/E with the constant vector field 

on Jac(C). The mapping ip intertwines the vector field rj defined by (j9.1Up — (j9.13|) with a 

constant vector field rj-^^^C^) on Jac(C), which up to a sign choice is determined by (jll.64p . 

dung). 

Proof For each p E P \ D there are four (not necessarily distinct) points m{i) G M/S, 
1 < « < 4 and jik) e Jac(C) \/f , 1 < k < A, such that tt^^^^ {m{t)) =p = tt^^^^^, {j{k)). The 
relation between the velocity field in the {u, f )-space and in Jacobi's elliptic coordinates, as 
described in Subsection 111 .31 yields that the set V{p) of the T„j(j) vr^^^^j, {^rn{i)) is equal to 

the set of the Tj(k) '^,u^(c) "^^^ P' of all p E P \ D such that V{p) consists of 

four distinct elements, of the form ±^', ±^", is equal to the complement of a closed curve in 
P\D. The set M' := (^t^j^j^^^ {P') and J' := (^t^j^^i^c)^ {P') is equal to the complement 

of a closed curve in M/S and Jac(C) \ J^, respectively. The mapping which assigns to 
m G M/S the pair 

is holomorphic from M/S to the tangent bundle T P of P. Its restriction 11^^, to M' is equal 
to an analytic diffeomorphism from M' onto a smooth two-dimensional submanifold X'p of 
T P. Similarly the mapping which assigns to j G Jac(C) \ the pair 

is holomorphic from Jac(C) \ J*" to TP, and its restriction 11^, to J' is equal to an analytic 
diffeomorphism from J' onto the same X'p. It follows that 

is an analytic diffeomorphism from M' onto J' . 

The branching properties imply that ip' has a continuous, hence analytic extension ip : 
M/E Jac(C) \ j^. Similarly {ip')'^ '■ J' —>■ M' has a continuous, hence analytic extension 

X '■ Jac(C) \ M/Tj. Because x ° is equal to the identity on M' and ip and x are 

continuous, we obtain that x o ?/; is equal to the identity on M/S. Similarly we obtain that 
Ip o X is equal to the identity on Jac(C) \ J^. □ 

Proposition II 1 . 61 implies that there exists a completion M/S of M/S, obtained by adding 
a curve at infinity which is isomorphic to the hyperelliptic curve C, such that the mapping ip 
in Proposition 111.61 extends to an isomorphism from M/S onto Jac(C). The vector fields ^ 
and rj extend to algebraic vector fields on M/S, which commute and are linearly independent 

at every point of M/S. The word "completion" is meant in the algebraic sense, but it can 
also be used in the sense that the fiows of and rj, with complex times, are complete on 
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M/S in the sense that they define a transitive action on M/S of the additive group C^. A 
completion with this property is unique up to isomorphism. 

Remark 11.7 In Proposition 110. ill we had obtained such a completion by adding the 
genus 2 hyperelliptic curve M^o/T, at infinity to M/S = M/S. There the completion is the 
complex torus M/S, in which M is a normalization of the projective closure M of M. The 
isomorphism of Proposition 111 . fj| leads to an isomorphism between Jac(C) and M/S and an 
isomorphism between C and the curve Mqo/S. 

The complex torus M is a 16-fold covering of the torus M/S ~ Jac(C), which can also 
be obtained by replacing the period lattice A of Jac(C) by 2A. See Remark [10.151 

Remark 11.8 The six fixed points of the hyperelliptic involution (A, z) ^ (A, —z) of 
C correspond to A = l//i, I//2, I//3, — 1/p, 7, 00, where 7 is given by pi.l7|) . pi.5|) . Of 
these, only 7 varies when we vary the constants of motion T and j. However, we are working 
on the assumption that ja = 0. The transformation of Subsection 19.21 to this case involved 
that we had to change the values of both T and p according to ()9.27|1 , and it follows that the 
isomorphism class of the hyperelliptic curve C varies in a two-dimensional subvariety of the 
three-dimensional moduli space of curves of genus two, as we vary the constants of motion 
T and j. If we also vary the moments of inertia li freely, then there is no restriction on the 
isomorphism class of the curve C. 

The fractional linear transformation A t-^ A/(l — p A) maps Jj = 1/ (/j + p) to 00 to 
— 1/p, the zero r = 2T/||j|p of p(A) to 7, and the zero 1/p of p{X) to 00. In this description 
of the zeros of p{X) we have used that js = 0. This leads to an explicit verification that the 
curve Moo/S in Remark 110.131 is isomorphic to C. 

Question 11.9 Can the constant vector fields on Jac(C) be determined more easily than 
in Subsection 111.31 by means of calculations at some special points, for instance at points 
corresponding to one of the fixed points of the hyperelliptic involution t? 

Question 11.10 (Richard Cushman) Could the isomorphism in Proposition 111.61 be ob- 
tained in a similar way as in Mumford [32, p. 3.57, 3.58] for Neumann's system? This 
question is suggested by the form ()3.14|) of the kinetic energy equation. 

11.5 The Translational Motion 

In order to obtain the motion of the point of contact p{t), we have to integrate the right hand 
side of ()3.H1 . Because this vector is horizontal, and the moment of momentum j around the 
point of contact is assumed to be horizontal as well, it is sufficient to determine the inner 
product with j and j x 63. For the latter one we have that 

{{A 00) X 63, 3 X 63) = {A 00, j) = {00, A-^j) = {00, v) = 2T, 
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where we have used fl3.8|) and ()3.12|) . It follows therefore that the component of p{t) orthog- 
onal to j grows linearly in time: 

{P{t), J X eg) = (p(0), J X eg) + 2rTt. (11.70) 

In order to determine the inner product of ^ with j in terms of Jacobi's elliptic co- 
ordinates, we begin with recalling the formula ()8.3j) . in which -^{p{t), j) is expressed in 
terms of the determinant of the vectors u, u, and lu. Under a transformation S" G S as 
in ()1U.3H) . the right hand side of ()8.3|) gets multiplied by €2- Therefore only the square of 
()8.3|) is a single- valued function on M/S, whereas fl8.3p is single- valued on the unbranched 
(= unramified) double covering M/Sq of M/S, where Sq denotes the group of 5* in p().3H) 
such that 62 = 0. In view of the isomorphism of M/S with Jac(C)/ in Proposition II 1 . 6( 

only j)^ can be a single valued function of Jacobi's elliptic coordinates. 
If C denotes the matrix {u, uj, luj), then 

({u, u) {u, oj) {u, I Oj) \ 
{uj,u) {uj,uj) {uo.Iuo) . 
{luO, u) {luj, LO) {I UJ, I UO) j 

Using that (m, u) = 1 and {loj, u) = 0, cf. ()11.29|1 . we therefore obtain that 

(detC)^ = {{00, 00) - {u, cj)2) -{lio, lio)- {io, Iuj)\ 
in which we can substitute 

{to, to) — {u, uj)"^ = {u X UJ, u X to) = {ii, it), 

cf. dSH). Substituting (fTT^ in combination with ()n.27|) with X = x, (jn.41|l . and ()n.4()jl . 
we obtain that 

^ ^ \2'\3' \[S{X2) S{X,))[S{X2) S{X2)J 

( A2 A2 _ A3 A3 \ I 
SiX2) SiXs)) j 

(A2 - XsY (-As^ - X2' - 2A2 A3) XI Aj _ (A2 - A3)' Aj Aj 

X2' A3' 5 (A2) 5 (A3) A2' A3' s (A2) s (A3) ■ 

Now we can write pi.43j) in the form 

(A2-A3)' A^ _^ (||j||'-2pT) (7-A2) 
A3'^(A2) p(-l/p-A3) ' 



94 



and we have a similar expression for A3, obtained by interchanging the indices 2 and 3. 
Substituting these, we obtain from (|8.3|) that 

,r ^ X(u) (det Cr ^ _ (pf -2pTf (.-A.) (o^-A3)^ 

\It'" p2 (I/P+A3) (I/P+A2) 

In combination with (jll.45|) . this leads to 

'd(p(r), j) 



dr 



-d' (7 - A2) (7 - A3) , (11.71) 



with 

nil + i/p) 



Remark 11.11 It follows from (111 .4411 with A = 7 and (jll.l9jl that 



- (7 - A2) (7 - A3) = ^ f 1 - E ^\ = ^ (1 - («, (1 - 7 1)-' u)) . (11.73) 

7\^ai-7y 7' 

The right hand side in pi.73|l is equal to zero if and only if n G Uc is in the branch locus 
for the branched covering (u, v) ^ u from L to f/c, cf. ()11.17|) . 

In fact, we could have concluded this at an earlier stage. The investigation of the branch 
locus started with the observation that this is the set of points where u, uj and I uj are 
linearly dependent, cf. (jll.4p . and according to (j8.3|) . this is equal to the condition that 
d(p(t), j)/dt = 0. 

In terms of the description of the real points u{t) on the unit sphere U in Remark lll.il 
the instances when u{t) reaches the boundary curves of the annulus in U coincide with the 
instances that (p(t), j) has a turning point. In view of the quasiperiodic motion of the 
rotational motion, this happens infinitely often. Keeping in mind that the j x es-component 
of the velocity of p{t) is equal to a positive constant, we obtain that the point of contact p{t) 
performs a swaying motion in the direction of j x 63. According to Corollarv l8.7l the function 
r I— >■ (p(t), j) is quasiperiodic if the rotational motion is not periodic and the irrational ratio 
^i(e) / ^2{^) mentioned after ()5.3p is sufficiently slowly approximated by rational numbers. 

As a function on the Jacobi variety Jac(C) of the hyperelliptic curve C, the function 

/ := - (7 - A2) (7 - A3) 

is rational, with poles along J^, zeros at J^^ and undetermined at the two points of 

n / = {0, / }. (11.74) 



'00 ^(7,0) J 00 

Note that J^'"'' is a point of order two in Jac(C), in the sense that 2 J^^'^'* = 0. 
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f has double zeros along J^,^ because when A is near 7, then on the hyperelliptic curve 

= P(A) = n (a. - A) i-l/p - A) (7 - A) 

i=l 

the variable 2; is a local parameter, in terms of which 7 — A is of order y'^. Also / has double 
poles along J*" , because if 6' is a local parameter for C at infinity, then A = 6^"^ u, z = v, 
where u and v are units, cf. Shafarevich [H^l IIL5.5]. Therefore, following a choice of a 
square root of / along curves in Jac(C), we find that there exists a rational function g such 

that / = g'^, either on Jac(C), or on a double covering Jac(C), unbranched, of Jac(C). Ben 
Moonen told me that g cannot be single- valued on Jac(C). Actually, if in the choice of the 
basis of A(C) corresponding to the closed curves Ai, A2, Bi, B2 as in Mumford |32l P- 3.76], 
with his point Oi in the role of (7, 0), then g turns into its opposite if we follow Ai around, 
and remains unchanged if we follow any of the other curves. Therefore g is a single valued 
rational function on the unbranched double covering 

jS^) = n\cy/Aoic), 

where Ao(C) is the sublattice of index two in A(C) which is generated by the basis elements 
corresponding to 2Ai, A2, Bi, B2. 

It follows that g is a rational function on Jac(C) with simple poles along the preimage of 
under the double covering : Jac(C) Jac(C), simple zeros along the preimage of J^^ q^, 
and undetermined at the preimage of pi.74|) . These properties characterize the function g 

on Jac(C) up to a constant factor. The function g can be identified with a quotient of theta 
functions as in Mumford (33 p. 3.80, 81], [SD Ch. II]. 

Question 11.12 We know that 

^ilhA = r X(«(r))V2 (A(r) u{r), 63 x j) (11.75) 
dr 

is a quasiperiodic function of r, because the rotational motion is a quasiperiodic function of 
r. We now have the additional information that the derivative of (p(t), j) is given hj dg 

along a straight line path in Jac(C), where g is the quotient of theta functions described 
above. Can this additional information be used in order to decide whether the problems 
with the integration (with respect to r) of the quasiperiodic function of r in the right hand 
side of ()11.75|) . which are mentioned in Subsection 15. 2| really occur? 



11.6 Chaplygin 

The last part of Chaplygin j^l §5], starting with "Now we discuss the curves traced out on 
the surface of the sphere . . . " , corresponds to our Remark 111.11 
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The proof that the elhptic coordinates ()11.21|) for u lead to a motion on the Jacobi variety 
of the hyperelliptic curve (|11.5(J|) with constant velocity is contained in Chaplygin "W, §3 up 
to (30)]. The notations in Chaplygin P §3] correspond to ours according to the following 
list, which is a continuation of the list in Subsection 13.51 

Chaplygin P §3 up to (30)] our notation 

(14) dSIZ I) 

(15) (jSSD, ^rrm and inmii 

(16) and (17) (fTTmUl 

(18) (fTT3T|) and (ITT3nj) 

(19) formulas following pi.40|) 
a^, 6^ c\d'^, g, k, a in (20) ai, as, ag (cf. firTTTHl V 1/p, 8T, 4 

X, z in (21) xi, X2, xs in pi.l9|l 

(22) (unsD, (dsn), (uni, dnsi, iin 

(23) formulas preceding pi.42p 

(24) (ITT:^ and (ITTTTj) 

(25) inrm 

(26) (Itt:^ 

formulas between (26) and (27) (fTTip . (ITm . (fm^ll 

(27), J, -(79Vj (da, 8T-4||j||Vp, 7 

(28) (imH|) 

(29) (fTT-53|l 

(30) (fTTITjl 

(33) (HHOD 

(34) (irrfTD 

In our Subsection 111.21 we have followed the miraculous calculations of Chaplygin ^ §3 up 
to (30)] quite closely, adding some more explanations in the hope to make these easier to 
read. In ()11.45|) we identified (up to a constant factor) the factor 

v/(A2 + 1/p) (A3 + 1/p) 

introduced in Chaplygin [HI (28)] with the integrating factor X('u)^/^ of Lemma mi the same 
as the factor y/X at the end of §2], and of Corollarv 18.41 In Subsection 111. 41 we added a 
discussion of the relation between the phase space of the rotational motion and the Jacobi 
variety of the hyperelliptic curve, about which Chaplygin did not say anything in jU]. 

Chaplygin did not tell how he came to the idea of using the elliptic coordinates (jll.21|) . 
The last part of Chaplygin §5] indicates that he had calculated the branch locus of the 
projection [u, v) ^ u from the [j, T)-level surface onto the u-sphere. Therefore he might 
have observed that in the complex domain it is equal to the union of two quadrics which, 
together with the sphere, belong to a one-parameter family of confocal quadrics, and this 
might have prompted him to use the elliptic coordinates pi. 211) . He might have refrained 
from mentioning this in his article, because in his time the use of elliptic coordinates in the 
presence of families of confocal quadrics was standard. 
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Although Chaplygin described the motion of u{t) in the annulus on the sphere in P| end 
of §5], he did not observe that u{t) reaches the boundary curves of the annulus precisely 
when d{p(t), j) / dt = 0, cf. Remark 111. Ill 

12 A Geometric Interpretation 

The kinetic energy equation in the form — X Z = 0, cf. ()3.14|) . is equal to the discriminant 
equation for the quadratic equation X{u) + 2Y{u, v) \ + Z{v) = in the variable A. Let 
a be an auxiliary parameter. Using ()3.9|) . the equation for A can be written in the form 

(Aw-f, {J + a){\u-v)) = {p-^ + a)\^-2aj3\ + 2T + a\\jf. (12.1) 

The discriminant of the right hand side is equal to zero if and only if 

(IIjIP-Js^) a'+{2T+\\jr/p) a + 2T/p = 0. (12.2) 

If (|12.2p holds, then the right hand side of (jl2.1|) is equal to 

(p^^ + a) {A - a js/ (p^^ + a) }^ 

and the equation (|12.ip is equivalent to 

{v, [p{I + p)-^ + a] v) = p-^ + a, (12.3) 

with V = 6 u — Tjv and 

e = \/ {X-aj^l [p-^ + a)] , T] = ll {X-aj^l [p-^ + a)] . 

It follows that 7] = {9 — 1) (p^^ + a) /aj^, and the conclusion is that the straight line / 
passing through m, with the direction vector equal to aj^u — (p^^ + a) v, is tangent to the 
quadric Q defined by the equation p2.3|) . Note that I = A~^{1), where I is equal to the 
straight line passing through —63, with direction vector equal to a J3 63 — (p~^ + a) j. Here 
the rotation A G SO (3) varies in the level surface in SO (3), corresponding to the equations 
()3.9|) and ()3.14|) . Also note that —Q = Q, and therefore the same properties hold with / 
replaced by — /. 

The equation ()12.2p has two solutions ai and 02, leading to two pairs of straight lines ±li 
and ib/2 and quadrics Qi and Q2 to which A~~^ i^h) and A~^ i^h) are tangent, respectively. 
The inner product of the two direction vectors is equal to 

(ai u - (p^^ + ai) V, a2 js u - (p^^ + 02) v) 
= {ai/p + a2/p + aia2) (||jf -Js') + IIjIIVp'- 

On the other hand it follows from (jl2.2|) that 

{\\jf-J3') (ai + «2)= -2T-\\jf/p, (||j|r-j3') «i«2 = 2r/p. 
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and we conclude that the direction vectors are perpendicular. Equivalently, the four straight 
lines ±/i, ±/2 form a rectangle. 

The above arguments work under the assumption that the moment j is neither vertical nor 
horizontal. If j 7^ is horizontal, js = 0, then the right hand side of ()12.H) is a square if 
a = — 2Tp/||j|p, and a constant if a = —l/p- In the first case the line passing through u 
with direction vector v is tangent to the quadric defined by ()12.3p . with a = — 2Tp/||j|p. 
In the second case the line passing through —v with the direction vector u is tangent to the 
quadric {v, [(/ + p)-' ~ p-'] v) =2T - Wjf/p. 

Question 12.1 What are all the straight lines / in the plane spanned by 63 and j, and 
quadrics Q, such that for each A G SO (3) in the level surface corresponding to the equations 
()3.9|) and ()3.14|) we have that A~^{1) is tangent to Ql This question may be related to 
Question 111.101 



12.1 Chaplygin 

Section [121 corresponds to Chaplygin's §5]. Chaplygin multiphed the figures by the radius 
r of the sphere, in order to have the corner point — r 63 of the rectangle ±/i, ±/2 attached to 
the point of contact of the sphere with the plane. 
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